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NASA/BLM APT: PHASE II FINAL REPORT 
VOLUME II “ TECHNOLOGY DEMONSTRATION 


1. INTRODUCTION . 

The NASA/BLM APT Final Report is divided into three 
volumes: Volume I - Executive Summary, Volume II - Technology 
Demonstration, and Volume III - Technology Transfer. Volume I 
is a summary of major aspects of the Phase II program, and 
Volume III is a sunj'tary of the workshops and project reviews 
performed to aid in transferring to BLM the technology demonstrated 
during Phase II operations. This volume. Volume II, is a technical 
discussion of activities and results of the technology demonstration 
aspects of Phase II. To avoid disrupting the flow of the project 
narrative, appendices have been used wherever possible to allow 
efficient inclusion of necessary backup data in support of the 
project activities. 


1 » i Overview o f Objectives for the Technology Demonstration . 

The primary objective of Phase II was to demonstrate a 
quantitatively-based integration of remote sensing data with 
existing ground data collection techniques for two applications: 
vegetation type mapping and productivity estimation. The 
vegetation types of interest are listed in the BLM-developed 
framework classification in Appendix 1-A. The productivity 
estimation was directed at the vegetation strata and associated 
estimates listed below: 


Strata 


Estimate 


Range 

Woodland 

Forest 


Pounds/acre (kilograms/hectare) of cattle 
forage 

Cubic feet/acre (cubic meters/hectare) of 
pinyon and juniper 

Board feet/acre of Ponderosa pine 
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A secondary objective was to demonstrate procedures 
for optimizing sampling while considering information requirements, 
precision specifications on the estimates and budget constraints. 
The approach, procedures and results presented in this report are 
a product of a combination of past experience and new knowledge 
and insight discovered and applied during project performance. 
Subsequently, it is anticipated this project will serve as 
another step in refining further the integration of advanced data 
collection techniques with existing methods for resource inventory 
and assessment. 

1 . 2 Overview of the Approach to the Technology Demonstration . 

Given the objectives above, an overall project plan 
was developed that designated the input data sources, the data 
to be measured or est-.mated from each source and the analysis 
procedures to be used in combining the data sets to produce the 
desired outputs. Appendix 1-b is a flow diagram that illustrates 
the many tasks and the interrelationships between these tasks for 
the Phase II technology demonstration. 

The inputs defined for Phase II were as follows: 

1. Landsat Digital Data - computer classification of the 
Landsat data would ptoduce a sis^ctral stratification of the test 
site to drive sampling and estimation efforts. Spectral categories 
would be described in terms of ground vegetation based on quanti- 
tative analysis of inputs 2 and 3 below and output in vegetation 
map form. 
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2. Large Scale Aerial Photography (LSP) - 35mm color 
transparency stereo pairs at nominal scales of 1:750 to 1:1200 
would be interpreted to estimate percent vegetation cover by 
species and to estimate tree volumes. The LSP was also used to 
select and locate ground plots. 

3. Ground Data - ground measurements on selected LSP plots 
would provide the basis for estimating productivity by strata 
(range, woodland and forest) when related to and combined with 
the LSP interpretation. 

4. Digital Terrain Data - National Cartographic Informa- 
tion Center (NCIC) terrain data tapes would provide elevation, 
slope and aspect information to describe Landsat spectral 
categories in topographic terms (thereby augmenting the vegeta- 
tion descriptions) . 

The following general data analysis procedures were 
to be applied during Phase II as follows: 

1. Regression Analysis - ground measurements and photo 
estimates would be combined to determine correlations and 
develop estimates of productivity. 

2. Analysis of Variance - vegetation composition by 
species of Landsat spectral categories would be derived from 
photo interpretation on the LSP. 

3. Contingency Analysis - Landsat spectral categories 
would be described in topographic terms (elevation, slope and 
aspect) when combined with digital terrain data. 
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The overall project approach can be more easily 
understood by examining the uiagrams in Appendix 1-B in view of 
the expected contributions of data inputs and the analysis 
procedures applied to the inputs. 

1.3 Technology Demonstration Participants . 

Individuals from a number of groups and agencies 
were involved in the accomplishment of the Phase II technology 
demonstration. The roles of the various participants is 
described in this section. 

BLM/Denver Federal Center, Scientific Systems Develop- 
ment (SSD) Group was the primary source of program requirements 
and acted as the liaison between the Arizona Strip District 
personnel, NASA-Johnson Space Center, and ESL Incorporated, the 
prime contractor. They also provided all the large-scale aerial 
photography used during the project. The major thrust for 
all phases of the NASA/BLM APT came from SSD and this assured 
that Phase II performance and output was in keeping with the 
long-term goals and requirement of BLM. 

BLM/Arizona Strip District personnel provided the 
detailed knowledge the study area which was vital in their 

assistance with the Landsat classification task. They developed 
the vegetation classification framework (Appendix 1-A) used in 
the vegetation description. And, importantly, they provided 
assistance and support during the field data collection efforts. 

NASA-Johnson Space Center provided both the overall 
project management and the majority of funding. It is to NASA's 
credit that the project management philosophy for this APT 
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1j2 

(including Phase II) was to allow BLM personnel to dictate the 
technical performance requirements. This approach ensured 
overall program relevance. 

ESL Incorporated was the prime contractor to NASA 
on Phase II, and was responsible for completing the major 
elements of the project: Landsat classification, sample selec- 
tion, photo interpretation, data reduction and analysis and 
output product generation. Assisting in the performance of 
these elements were three subcontractors : Resource Inventory 
Services (RIS) , San Jose, California; Remote Sensing Research 
Program (RSRP), U.C. Berkeley; and Rangeland Resources Incorporated 
(RRI) , Mancos, Colorado. RIS provided consultant services in 
quality assurance with particular attention to sampling procedures. 
RSRP developed initial sample design recommendations and performed 
sample allocations utilizing a survey planning model. RRI 
collected the ground data used in the productivity estimation. 

The remainder of this report will describe each of 
the major project elements listed above in detail. 


2 . 


METHODOLOGY . 


2 . 1 Landsat Data Processing . 

The Landsat data processing for this project was 
divided into six major categories: 

(1) Preprocessing steps 

(2) Training of the classifier 

(3) Maximum likelihood classification 

(4) Class description 

(5) Digitizing 

(6) Digital terrain data. 

The following sections will discuss in detail the procedures 
followed in each of these categories. 

2. 1.1 Preprocessing Steps . 

Landsat digital data requires a certain amount of 
manipulation befojre it can be used effectively and efficiently 
in a resource assessment/inventory project. To be used more 
effectively, the data is placed into the proper format, corrected 
for radiometric and geometric distortions, and related to a 
preselected ground coordinate system. To be used more efficiently, 
the project area is extracted from the appropriate Landsat scenes 
to eliminate unnecessary manipulation of data that is outside 
project boundaries. 

2. 1.1.1 Tape Reformatting and Study Area Extraction . 

The Arizona Test Site was located entirely within 
Landsat scene ID #20947-17074, taken on 26 August 1977; so the 
computer compatible tape (CCT) for that scene was entered into 
ESL's IDIMS using the function ERTSENTR. The study site was 
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located within the data by bringing up the Landsat image on 
the IDIMS color display and comparing it with existing maps of 
the area. The appropriate subsection of the scene covering 
the site was identified and extracted for further work and 
thus reduced subsequent processing costs by eliminating unnecessary 
data. Tape reformatting in preparation for subsequent processing 
was then performed using IDIMS function PLENTER. This function 
converted the digital data from the EROS format to the IDIMS format, 
removed the synthetic pixels, mosaicked the four image strips 
together, and associated the SIAT (Special Image Annotation Tape) 
file with the data The SIAT file is required for the geometric 
rectification of the data while the other elements from PLENTER are 
necessary prior to the radiometric correction steps. The 
reformatted scene subsection was then ready for the remaining 
preprocessing steps. 


2.1.1, 2 Radiometric Correction. 


The radiometr,' corrections to the data remove the 
"consistent striping" caused by slight differences in sensitivity 
of the scanner detectors on-board the satellite and also 
eliminate "intermittent striping" or bad data lines. "Consistent 
striping" was detected by viewing the study area at full image 
resolution on the IDIMS color display. No bad data lines were 
detected at this point; however, one was noted during later 
steps (see Section 2.1.3). The IDIMS function HISTEQ was used to 
generate frequency histograms for each sensor (six) and for 
each band (four), 24 histograms in all, to quantify the extent 
of banding by intensity value by sensor. Utilizing this infor- 
mation, the Landsat data was passed through HISTEQ again and 
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the appropriate adjustments made to the data to generate uniform 
data for each sensor for each band. The result of this processing 
was a radiometrically corrected Landsat image for the Arizona 
Test Site, an intermediate product in the sequence of steps 
towards a multispectral classification of Landsat digital data. 

2. 1.1. 3 Geometric Rectification . 

Two types of geometric rectification were required 
to orient the Landsat data properly with the ground: discontinuous 
error correction and continuous error correction. The discon- 
tinuous error correction accounted for anomalies in the Landsat 
data caused by satellite and multispectral scanner conditions 
noted at the time the image was first taken and required use of 
the SIAT file mentioned in Section 2. 1.1.1. The continuous error 
correction registered the Landsat data to the selected ground 
coordinate system, generated square pixels (in raw data they are 
rectangular, nominally covering 79 meters long by 56 meters 
wide on the ground) and required use of a ground control point 
network developed specifically for the area. 

2. 1.1. 3.1 Discontinuous Error Corrections . 

The correction of errors using the SIAT file actually 
occurred when precision geometric correction (PGC) was performed 
on the classified Landsat data (Section 2.1.3). However, those 
error sources, as listed below, are briefly described in this 
section since they are part of the overall geometric rectification 
required to produce useful Landsat classification products: 

(1) Mirror scan velocity profile 

(2) Panoramic distortions 

(3) Scan skew 

(4) Perspective geometry 

(5) Attitude (of spacecraft) . 
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The effects of these error sources are discussed below and infor- 
mation derived from the SIAT file is used to correct them during 
PGC. 

Mirror Scan Velocity Profile . There is a scanning 
mirror in the Landsat Multispectral Scanner (MSS) that 
deflects the line-of-sight to produce the scan. 

Its velocity is only approximately linear. The 
deviation causes a fairly repeatable error in the 
along-scan direction with a displacement of pixels 
that varies approximately +4 00m. 

Panoramic Distortion . The region of the projection 
plane image is proportional to the tangent of the 
scan angle rather than the scan angle itself, which 
produces an along-scan distortion. 

Scan-Skew . The satellite moves 216m during the 
active period of one scan. Thus the swath scanned 
is not perpendicular to the ground track, so the 
X and y axeis of the scene are not perpendicular. 

The error is across-scan, and not related to the 
earth's rotational skew. 

Perspective Geometry . MSS data is collected in a 
perspective sense; that is, all imaged pixels are 
projected to the tangent plane via lines that meet 
at a point (the detector) . Because the MSS is a 
line scanner rather than a framing system, the 
effect is in the along-scan direction only. 

Attitude . The sensor vertical axis is nominally 
normal to the earth's surface, and one horizontal 
axis is along the velocity vector. Maximum pitch 
and roll values are 7 mrad. Maximum yaw is 10 mrad. 
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Maximum rate is 0.26 mrad/sec or 7 mrad/scene for 
each angle. Pitch induces an across-scan distor- 
tion. Yaw rotates the scene and is equivalent to 
an along-Bcan error and an across-scan error. 

Roll is an across-scan effect. Attitude changes 
slowly and may be treated as constant for one 
swath, but induces a spatially varying distortion 
over the scene. 

Note ; Certain discontinuous geometric distortions within 
the Landsat data had already been corrected during tape 
reformatting using known models and parameters. They included: 


The 24 detectors are 


sampled sequentially in such a manner that each 
line out of each group of six lines should be 
offset by 1/12 pixel from the one above it. The 
distortion is distributed in the along scan 
direction. 


Earth Rotation . The earth rotates 13 km in the time 
one scene is recorded (at the equator) . Thus, there 
is a gradual westward shift of the swath viewed. The 
error is along-scan, and acts on six lines of the 
image at a time, because there are six detectors in 
each spectral band. 


Synthetic Pixels . For data processing purposes, 
NASA GSFC inserts synthetic pixels within each 
line to make all lines the same length, which is 
a multiple of 24. This introduces discontinuities 
in the along-scan direction which must be removed 
in the geometric correction process. 


2*'5 


2.1. 1.3.1 —Continued. 


Spectral Band Offset . The four groups (four 
spectral bands) of detectors are offset 112 meters. 
However, this offset is corrected by the insertion 
of filler pixels at the start end of the scan lines 
on the CCT's produced by NASA. 

2. 1.1. 3. 2 Continuous Error Correction . 

The errors that were corrected by using a ground 
control point network included: 

(1) Spacecraft velocity 

(2) Altitude 

(3) Map projection. 

The effects of these errors are discussed briefly below. It is 
important to note that the control point network used for this 
geometric rectification element was also used to overlay the 
digitized boundaries of both the study area and the allotments 
to be sampled after classification (Section 2.1.5). 

Spacecraft Velocity , Although the satellite is 
in a very nearly circular orbit, there are slight 
velocity variations. These variations cause an 
across-scan distortion which tends to introduce 
an overlap of scan swaths if the velocity is too 
low, and a gap if it is too high. 
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Altitude . Deviation of the satellite from the 
normal altitude produces scale change distortion. 

The altitude will vary between 900 and 950 km 
during the mission lifetime. The major effect will 
be a change in scale from scene to scene. The 
second order effect is due to the altitude variation 
within a scene. 

Map Projection . To maximize its utility, the 
Landsat data should be cast into the same projection 
as a particular map or set of maps. Doing this is 
equivalent to inducing a distortion in the data 
and can be achieved by the same processes that 
remove distortions. 

The control points to be used for continuous error 
correction were located first within the Landsat data by selecting 
those distinct features viewed on the IDIMS color display that 
could also be located on USGS 7.5 minute qua'^r;: sgle maps of the 
area. The Landsat sample and line value for each selected 
point (a pixel) was placed in a computer file and stored on disk. 
This file, named SLARIZ .ARIZ .GIS , eventually contained 
95 points with corresponding sample/line values. As each point 
was selected cn the screen, its location was also annotated on 
the 7.5 minute quads for digitizing later. An IMMS software 
subsystem, GES (Geographic Entry System) was used to generate 
the five control point transform^^ cions necessary to relate the 
Landsat data to the ground. Table 2-1 lists the transformations 
(and their applications) and the following paragraphs describe 
the procedures followed to generate them. 
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Table 2-1 

. NASA/BLM APT 

Control 

Point Transformations. 

Transforma- 
tion Name 

Source/ 

Destination 

Order 
of Fit 

Applications 

1. ARIZTl 

GES-to-Landsat 

Third 

-to evaluate control point 
network distribution and 
residuals 




-to determine locations of 
1978 large scale photo- 
graphy transects on 
uncorrected Landsat data 

2. ARIZT2 

GES-to-50 Meter 
Grid 

Second 

-to pass the data retained 
in GES files to the 50 
meter grid in the form of 
a strata mask file 
(digitized boundaries, etc.) 

3. ARIZT3 

50 Meter Grid- 
to-Landsat 

Third 

-to document the transforma- 
tion to a ground coordinate 
system 

4. ARIZT4 

50 Meter Grid- 
to-UTM 

First 

-to determine UTM coordinates 
for areas or points selected 
for sampling (on corrected, 
PGC, image) 

5. PGC 

Lar.dsat-to- 
50 Meter Grid 

Third 

-to register uncorrected 
Landsat image to 50 meter 
grid using SIAT file data 
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GES utilizes a geoblock concept for establishing the 
area on the earth's surface (specified in latitude and longitude) 
where all geographically oriented data for a given project is 
located. The Shivwits Planning Area and all control points weri 
referenced to a geoblock whose origin was set at 35® 30 '00" N 
and 114®30'00" W and covered 105 minutes north and 105 minutes 
east. This geoblock and its associated files contain all the 
digitizing work for this project. The 7.5 minute quads with the 
control points were registered to the Shivwits geoblock based on 
the latitude/longitude coverages of each map. The 95 control 
points were then digitized off these maps and automatically 
stored in GES files. A third order least squares transformation 
(ARIZTI) was generated between the GSS files and the Landsat 
file SLARIZ using the TRNSFORM program. Through an iterative 
process of rejecting points that had excessively large residuals 
and rerunning the transformation, a final group of 71 of the 
original 95 control points were accepted and had residual means 
of .7 pixels in the sample direction e,nd .49 pixels in the line 
direction. 

Next, the program ALLCOORD was used to read the GES 
files containing the digitized control points and create a file 
of control point values relative to a 50 meter grid 
(ARIZ50M.ARIZ.GIS) corresponding to UTM Grid Zone 12 with the 
50 meter grid origin (bias) set at 4,100,000 N, 220,000 E. The 
UTM grid origin was obtained by extending the zone 12 values 
west to include that portion of the study area that occurs 
west of 114® longitude. The program TRNSFORM was used to develop 
a second order least squares transformation (ARIZT2) between 
the GES internal units and the 50 meter grid values in ARIZ50M. 
The same 71 control points used in transformation ARIZTI were 
used in ARIZT2 and resulted in residual means of .002 meters in 
the X direction and .006 meters in the y direction. The residual 
means for this transformation are so small because this is a 
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simple translation from one rectangular coordinate system 
(GES internal units) to the rectangular UTM grid. 

A third order least squares transformation (ARIZT3) 
was generated next between the 50 meter grid locations (file: 
ARIZ50M) and the Landsat sample/line values (file: SLARIZ) for 
the 71 control points. This transformation was essentially a 
duplicate of the transformation generated duriiig the precision 
geometric correction (PGC) of the classified Landsat data 
(Section 2.1.3); the only difference was that ARIZT3 did not 
take account of any of the SIAT file data. Since PGC does not 
keep a permanent record of this transformation during image 
registration, ARIZT3 was generated simply for documentation. 

The final transformation (ARIZT4) was generated by 
first using ALLCOORD to create a text file of the control points 
referenced to UTM Zone 12 (ARIZUTM. ARIZ .GIS) . The 50 meter grid 
file, ARIZ50M, was specified as the source and the UTM file, 
ARIZUTM, was specified for the destination of the transformed 
point locations (1st order transformation) . Residual means 
for ARIZT4 were .57 meters in the x direction and 6.07 meters 
in the y. The residuals are so small because the transformation 
is a simple translation from the subsectioned UTM grid (the 
registered image) to the UTM coordinate system (zone 12). 

All transformations discussed above (except PGC) have 
been stored in the file ARIZTMTX. ARIZ .GIS and documented with 
computer line printer outputs, CRT hardcopies and session 
histories for delivery to BLM. 

As a point of reference, the radiometric corrections 
were applied prior to scene classification; however, the 
geometric transformations were not applied until after classifi- 
cation (see Section 2.1.3). 
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Training of the Classifier . 


Unsupervised and supervised training procedures were 
utilized to generate the spectral statistics to be applied in 
the maximum likelihood classification algorithm used to classify 
the Landsat data. The supervised training was done to improve 
the results of the preliminary classification which was based 
on unsupervised (clustering) procedures. The overall training 
procedures generated 83 spectral clusters which were used to 
classify the raw Landsat data. 

2.1.2. 1 Unsupervised Training . 

The initial unsupervised training was based on 
clustering each of the 10 x 30 pixel "windows" which were 
located and extracted for each of the 120 large scale aerial 
photography transects flown in 1978 by the BLM. These 
transects consisted of ten "nested" photo plots flown in an 
east-west direction with two 35mm cameras to produce ten 
separate stereo pairs at 1:2000 scale from one camera, plus one 
additional photo for each pair at 1:16,000 scale. The 10 x 30 
window approximated the area covered on the ground by the ten 
frames taken at the smaller of the two scales. These windows 
were first stratified into four major environmental types (low 
desert, high desert, forest and woodland), and then the combined 
flight lines falling in a given environmental stratum were 
clustered using the IDIMS function ISOCLS. (The parameter values 
used in the ISOCLS runs were DLMIN=3, SEP=0, STDMAX=3.2, 

IST0P=15, LNCAT=1, NMIN=30, MAXCLS=30, KRN=10 and CHNTHS=3.0). 

These spectral clusters were then given initial 
descriptions through photo interpretation of the transect 
photography. First, line printer maps of the clustered data 
were generated for each 10 x 30 "window" as shown in Figure 2-1 
BLM personnel from Arizona photo interpreted the vegetation 
type for each cluster within a "window" by projecting the transect 
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Figure 2-1 . Annotated Line Printer Hap from ISOCLS 
of a 1978 LSP Flight Line*. 
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photography for that window onto a viewing screen and determin- 
ing the individual cluster descriptions through group discussion. 
Figure 2-2 is an example of the photo interpretation form used 
by the BLM in generating the preliminary cluster descriptions. 

As was expected, this procedure resulted in a large number of 
clusters and also a certain amount of duplication of clusters. 

The total group of spectral clusters was edited to 
reduce both the number as well as the duplicity of clusters 
to establish the initial statistics package for maximum 
likelihood classification. A cluster was deleted when both 
its spectral characteristics and vegetation description was 
similar to another cluster. IDIMS functions COMPARE and DIVERGE 
were used to determine spectral similarities. Figure 2-3 
is an example of a COMPARE plot of spectral clusters. The 
cluster number is given in the center of each cluster represen- 
tation; and the location, shape and size of each representation is 
a graphical illustration of the means and standard deviations 
of the spectral responses for the clusters as calculated for 
each of the Landsat bands noted (IDIMS band 2 == Landsat band 5 
and IDIMS band 4 = Landsat band 7) . The plot is used to 
determine spectral overlap between clusters. Figure 2-4 is an 
example of a DIVERGE output for spectral cluster distance 
analysis. The matrix of numbers represents a measure of the 
distance between cluster pairs (intersections of rows and 
columns are the pairings) where the larger the number is 
between any two clusters, the more separable the clusters are. 

This analysis considers the statistics for all four Landsat bands 
simultaneously as opposed to COMPARE which is two bands at a 
time. In this project any DIVERGE value greater than 1300 was 
deemed acceptable; cluster pairs with values below 1300 were 
examined on COMPARE plots for overlaps in statistics and checked 
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for similarities in PI descriptions as well prior to deciding 
which clusters to delete. This editing process resulted in an 
initial statistics package with separate training statistics 
for each environmental strata and was used to perform test 
classifications on four specifically selected subsections 
(Intensive Study Areas) of the study site. 

The Intensive Study Areas (ISA's) were identified and 
located so as to have one ISA for each of the four environmental 
strata: low desert, high desert, woodland and forest. Each 
ISA was 512 X 512 pixels in size and their locations with 
respect to the Arizona Test Site are illustrated in 
Figure 2-5 . Since a separate set of statistics was developed 
for each strata, each ISA was classified individually with the 
appropriate training statistics using the IDIMC maximum likeli- 
hood classification algorithm CLASFY. Classification results 
for the ISA's were viewed on the IDIMS color display and from 
this evaluation it was determined that both water and Ponderosa 
pine types required additional training statistics to produce 
a more accurate classification. 

2. 1.2. 2 Supervised Training . 

The required additional training statistics were 
generated using supervised training procedures. In contrast 
to the unsupervised approach above where spectral clusters 
were described after clusteri'ng, these procedures were based 
on ground conditions being identified before clustering. Homo- 
geneous areas of water and of Ponderosa pine were identified 
and located within the unclassified Landsat image and then 
extracted for clustering by the IDIMS function TSSELECT. The 
spectral cluster statistics generated from this operation were 
then combined with the existing environmental strata statistics 
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to create a single training set to be used for the full Arizona Test 
Site image classification. First, however, the now combined 
set was edited (same procedures as described previously) 
to eliminate duplicate clusters that had occurred between 
environmental strata training sets. 

The edited statistics package along with the Arizona Test 
Site Landsat image subsection were input to ESL's Advanced 
Scientific Array Processor (ASAP) on the In-house IBM370 
computer for running the maximum likelihood classification 
algorithm on the data set. This classification was controlled 
by a five percent "threshold” to readily identify any areas that 
did not fit within the training (spectral cluster statistics) . 
Essentially, the threshold value specified that assignment to a 
spectral class for any pixel required a .95 or greater likelihood 
of being correct in order to be accepted. All class assignments 
with likelihood values below .95 were rejected, and those pixels 
were assigned a value of zero in the classified image (refer to 
ESL Technical Manual TM705, Rs‘v. A, Vol. I, pp. 4-55 to 4-62). 

It is important to note that not all rejected class assignments 
were necessarily incorrect. A small number of low likelihood 
but correctly classified pixels were undoubtedly rejected along 
with the greater number of incorrectly classified ones. 

Evaluation of the thresholded classification on the 
IDIMS color display, with the assistance of BLM Arizona person- 
nel, identified the requirement for additional training 
statistics to be generated for sedimented water, shadow, 
barren land and low desert vegetation. Additionally, a "bad" 
data line from the raw Landsat input image was detected in 
band 6. The additional statistics were generated through 
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supervised training procedures wherein BLM Arizona personnel 
specified areas of known ground conditions corresponding to 
those types missing from the preliminary classification. 

TSSELECT was used to outline and cluster those areas, and the 
resulting statistics were combined with the previous training 
set. After final editing of this combined set (a third cycle 
of editing) , a final spectral cluster training set made up of 
83 classes was available for final classification of the 
Arizona Test Site image. 

2.1.3 Maximum Likelihood Classification and Precision 

Geometric Correction . 

Before this step in the Landsat processing could be 
performed, the "bad" data line identified above had to be 
removed. The IDIMS function FIXLINE was used to replace the 
line with a new line formed by averaging the gray levels from 
the lines above and below the "bad" one on a pixel-by-pixel 
basis across the entire length of the image. The resultant 
image was acceptable for the final maximum likelihood classifi- 
cation run. 

The 83 spectral clusters (training statistics) and 
the "clean" Arizona Test Site Landsat image were submitted 
to the IBM370-based ASAP f_r classification. The classified 
output image was checked on the IDIMS color display and approved 
for further processing. To this point, all processing had 
been performed on uncorrected Landsat data, i.e., the data had 
not yet been registered to the ground even though all of the 
necessary coordinate transformations to do so had been built. 
Precision Geometric Correction (PGC) of the approved classifi- 
cation accomplished this registration and resulted in an image 
ready for sampling. 
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As discussed in Section 2. 1.1. 3 (Geometric Rectifi- 
cation) , a control point network and associated coordinate 
transformations were developed to prepare for the precision 
geometric correction procedures. IDIMS function PGC was used 
to combine the necessary information from the SUiT file, the 
appropriate control point files (image and map) related to 
the 50 meter UTM-based grid and the uncorrected Landsat 
classified image to produce a registered Landsat classification 
image. The PGC produced image was rotated to a north-south 
orientation, resampled (nearest-neighbor) to a 50 meter square 
pixel resolution and related to the UTM zone 12 ground coordinate 
system. The next step in the processing was to describe and 
combine the 83 spectral classes into preliminary summary 
vegetation classes. 


2.1.4 


Class Description . 


The Landsat classification, to be used as a strati- 
fication for directing the multistage sampling efforts, required 
further analysis and processing. First, dissimilar vegetation 
types that appeared to have the same spectral response (i.e., 
placed in the same sp'^ctral cluster) had to be separated. 

Then, after thesj spectral "confusions" had been eliminated, the 
resulting classes had to be grouped into various summary 
categories of similar vegetation and the categories had to be 
described in terms of species composition on the ground. 

Grouping the spectral classes was necessary because the budget 
was not large enough to adequately sample the individual 
classes. Sample allocation and selection based on the similar 
vegetation category groupings, however, were supported with 
species composition descriptions for each based on analysis 
of variance performed on the 1978 LSP photo interpretation 
results as applied to the Landsat classification. The following 
sections describe the elimination of "confusion" classes, 
grouping of the classes and descriptions of those classes for 
sampling . 

2. 1. 4. 1 Reduction of "Confusions" Between Classes . 

The approved final Landsat classification for the 
Arizona Test Site consisted of 83 spectral classes. BLM/Arizona 
personnel evaluated these classes as to their locations and 
supplied to ESL a class-by-class list of "confusions." It was 
decided that an environmental stratification based on elevation 
would be used to eliminate these "confusions." BLM/Arizona 
personnel supplied ESL with a strata map for the area (see 
Figure 2-6 ) based on delineation of contour lines off 

1:250,000 scale quad sheets as follows: 
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Figure 2-6. Environmental Strata Kask Used 
to Reduce Spectral Confusion. 
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(1) Low Desert - below 4000' contour line 

(2) High Desert - below 6000' contour line but above 4000' 

(3) Mountain - above 6000' contour 

(a) Known locations of Ponderosa pine stands 

(b) Remaining sites with no Ponderosa pine stands. 

The original 83 spectral classes were grouped first 
into 26 summary classes as listed in Table 2-2 . The groupings 
were based on the original preliminary names given the clusters 
during the training phase (Section 2.L2.1) and also on 
similarities in mean spectral response for each band from the 
final training statistics. Within each summary class the 
"confusions" noted by BLM were checked on a spectral class-by- 
spectral class basis giving special attention to the location 
of each class with respect to the enviromnental strata (as 
viewed on the IDIMS color display) throughout the Arizona Test Site. 
The following general rules were applied: 

(1) If the spectral class was located in the Low Desert 
strata and was assigned a preliminary name indicative 
of a Low Desert type, no change was made. 

(2) If the spectral class was located in the Low Desert 
strata but was assigned a preliminary name indicative 
of either High Desert or Mountain types, a new 
spectral class identification number was assigned 

to only those pixels within the Low Desert strata 
and then that ID number was placed in the Low Desert 
grouping it most closely matched based on mean 
spectral response (all four Landsat bands) . 

(3) If the spectral class was located in the High Desert 
strata and was assigned a preliminary name indicative 
of a high desert vegetation type, no change was made. 
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Table 2-2 . Results of Assignment of 83 Spectral 
Clusters to Summary Categories. 



Strata 

Description 

Spectral Clusters 

1 

Low Desert 

Creosote-Bursage 
(rocky soil) 

1,2,7,10,11 

2 

Low Desert 

Creosote-Bursage 
(sandy soil) 

3,5,9 

3 

Low Desert 

Creosote-Pure 

4 

4 

Low Desert 

Upland Desert 
Shrub - Creosote 
Dominant 

6,12,13,14 

5 

Low Desert 

Blackbrush 

45 

6 

Low Desert 

Mixed Desert 
Shrub - Creosote 
and Cactus 

8,74 

7 

All 

Riparian 

Woodland 

25,26,27 

8 

High Desert 

Shrub-Grass 

28,34,38 

9 

High Desert 

Grassland-Shrub 

31,39 

10 

High Desert 

Snakeweed-Grass 

29,30 

11 

High Desert 

Sage-Mix Shrub 

32,44 

12 

High Desert 

Sage 

35,49,51,54,55, 

13 

High Desert 

Saltshrub 

36,37 

14 

High Desert 

Pinyon- Juniper- 
Shrub 

33,50 

15 

High Desert 

Pinyon-Juniper- 

Sage 

53,56,57,58,59,61, 

62,63 

16 

High Desert 

Pinyon- Juniper 

40,41,46,52 

17 

Mountain 

Mountain Shrub 

42 

18 

Mountain 

Mixed Chaparral 

60 

19 

Low Desert 

Agriculture 

15,16,17,18,19,20, 

21,22,23,24 

20 

Mountain 

Ponderosa Pine- 
Oak 

66 
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St rata 

Description 

Spectral Clusters 

21 

Mountain 

Ponderosa Pine- 
Mix 

65 

22 

Mountain 

Ponderosa Pine 

64 

23 

All 

Shadow 

47,83 

24 

All 

Water 

67,68,69,70,71,72, 

73,82 

25 

All 

Bare 

76,77,78,79,80,81 

26 

Low Desert 

Upland Desert 
Shrub-Blackbrush 

43,48,75 







2. 1.4.1 


Continued. 


(4) If the spectral class was located in the High Desert 
strata but assigned a preliminary name indicative of 
either Low Desert or Mountain types, a new spectral 
class identification number was assigned to only 
those pixels located within the High Desert strata 
and then that new ID number was placed in the High 
Desert grouping it most closely matched based on mean 
spectral response (all four Landsat bands) . 

(5) If the spectral class was located in the Mountain 
strata and was assigned a preliminary name indicative 
of a Mountain type, no change was made unless the 
class was a Ponderosa pine type and was located 
outside of a designated Ponderosa area. Reassignment 
of numbers and classes for Ponderosa types in non- 
Ponderosa areas followed similar guidelines as in (2) 
and (4) above. 

(6) If the spectral class was located in a Mountain 
strata but had been assigned a preliminary name 
indicative of other than a Mountain type, the same 
reassignment procedures were followed as in (2) and 
(4) but applied in the Mountain strata. 

The rationale behind these rules was that a spectral 
category could conceivably represent more than one vegetation 
type and these different vegetation types could be accounted 
for by assignment of new class ID's to vegetation types not 
represented in the final training statistics. As a result of 
the reassignment procedure, the original 83 spectral clusters 
were expanded to 117 classes and grouped into 27 summary categories 
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(orioinally 26). The 117 classes eliminated the "confusions" 
noted by BLM and the summary categories provided a framework 
for efficient sampling for productivity estimation. The 27th 
summary category was created by combining four of the reassigned 
classes in the High Desert strata that did not match well with 
any of the other 26 original categories but exhibited very similar 
mean spectral responses; and were all generally located within the 
same area indicating the potential of all being of a similar type. 
Table 2-3 illustrates the reassignment procedures and Table 2-4 
lists the overall results by summary class. 

2. 1.4. 2 Summary Class Description with 1978 LSP Data . 

Once the summary classes were created, quantitative 
descriptions of those classes were developed for sample alloca^ 
tion optimization by using the photo interpretation results 
from the 1978 LSP ( supplied by the BLfI ) and analysis 
of variance procedures. The individual LSP plots were paired 
with their corresponding Landsat spectral classes based on 
location; the photo interpretation associated with each plot 
(% cover) was entered into ESL’s Earth Resources Inventory 
System (ERIS) on IDIMS and a one-way analysis of variance was 
run on a combined data set of photo interpretation/Landsat 
summary class (groupings of spectral classes) to produce means 
and covariance matrices of % cover by species for each summary 
class. The following paragraphs describe these procedures in 
more detail. 

The map locations of each photo plot (approx. 1200 plots) 
were provided to ESL by BLM. The Landsat spectral class 
corresponding to each plot was determined by first digitizing 
the plot locations and then transforming the list of map 
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Table 2-3 . Reassignment of Spectral Classes Where 

"Confusions" Exist Between Environmental 
Strata. 


ORIGINAL ASSIGNMENT REASSIGNMENT TO; 


Preliminary 

Name 

Spectral 

Class 

Summary 

r 

Spectral 

Class 

Summary 

Category 

Cateqory’ 

Strata 

"Confusion” Noted 

P. Pine-Mix 

65 

21 

Mountain 

Found in High Desert 

85 

15 

P. Pine 

64 

21 

Mountain 

Found on Poverty Mtn.< 

84 

15 

P. Pine Mix 

65 

21 

fl 

II II If fl 

85 

15 

Sage 

51 

12 

lighDesert 

Found in Low Desert 

86 

2 

II 

35 

12 

II II 

II II If II 

87 

1 

H 

49 

12 

II fl 

II II fl fl 

88 

1 

II 

54 

12 

II II 

fl II II fl 

89 

1 

ft 

55 

12 

II II 

II II II II 

90 

1 

Sage/Mlx. Shrub 

32 

11 

lighDesert 

Found in Low Desert 

91 

1 

Up. Des. Shrub 

75 

26 

LowDesert 

Found in High Desert 

117 

27 

fl 

45 

If 

II II 

M II II If 

101 

27 

Agriculture 

15-24 

19 

LowDesert 

Found in High Desert 

113 

18 

Mix. Chapparal 

60 

18 

Mountain 

Found in Low Desert 

114 

19 

Saltshrub 

36 

13 

HlghDeser 

: Found in Low Desert 

92 

1 

tf 

37 

If 

II II 

II II II II 

93 

1 

Snaxeweed /Grass 

29 

10 

MlghDeser 

: Frand in Low Desert 

94 

2 

II It 

30 

II 

fl If 

II II II II 

95 

2 

Grassland/ :»hrub 

31 

9 

lighDeser 

Found in Low Desert 

96 

4 

ft ft 

39 

II 

II II 

II II II II 

97 

4 

Shrub-Grass 

28 

8 

lighDeser 

Found in Low Desert 

98 

1 

ff ft 

* 34 

8 

II II 

fl It If ff 

99 

2 

If ff 

38 

ff 

fl If 

II If ff If 

100 , 

1 

4ix . Deser tShrub 
(Creosote&Cactu 

8 

i) 

6 

LowDeser 

t Found in High Desert 

115 

27 


*Note: Mo Ponderosa Pine Is found on Poverty Mountain even though it is above 6000 feet. 
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Table 2-4 


Results of Assignment of All Spectral 
Clusters to Summary Categories. 


H^nQ^rejUII 

■' ^ 

Strata 

Descrip^ion 

Spectral Clusters 

1 

Low Desert 

Creosote-Bur sage 
(rocky soil) 

1,2,7,10,11»87,88, 

89,90,91,92,93,98, 

100 

2 

Low De&cit 

Creosote-Bur sage 
(sandy soil) 

3,5,9,86,94,95,99 

3 

Low Desert 

Creosote-Pure 

4 

4 

Low Desert 

Upland Desert 
Shrub - Creosote 
Dominant 

6,12,13,14,96,97 

5 

Low Desert 

Blackbrush 

45 

6 

Low Desert 

Mixed Desert 
Shrub - Creosote 
and Cactus 

8,74 

7 

All 

Riparian 

Woodland 

25,26,27 

8 

High Desert 

Shrub-Grass 

28,34,38 

9 

High Desert 

Grassland-Shrub 

31,39 

10 

High Desert 

Snakeweed-Grass 

29,30,106,107,108 

11 

High Desert 

Sage-Mix Shrub 

32,44,102,103,105, 

111 

12 

High Desert 

Sage 

35,49,51,54,55, 

110,112 

13 

High Desert 

Saltshrub 

36,37,104,109 

14 

High Desert 

Pinyon- Juniper- 
Sage 

33, 50 

15 

High Desert 

Pinyon- Juniper- 
Shrub 

53,56,57,58,59,61, 

62,63,84,85 

16 

High Desert 

Pinyon- Juniper 

40,41,46,52 

17 

Mountain 

Mountain Shrub 

42 

18 

Mountain 

Mixed Chaparral 

60,113 

19 

Low Desert 

Agriculture 

15,16,17,18,19,20, 

21,22,23,24,114 

20 

Mountain 

Ponderosa Pine- 
Oak 

66 
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liQjmjM 

Strata 

Description 

Spectral Clusters 

21 

Mountain 

Ponderosa Pine- 
Mix 

65 

22 

Mountain 

Ponderosa Pine 

64 

23 

All 

Shadow 

47,83 

24 

All 

Water 

67,68,69,70,71,72, 

73,82 

25 

All 

Bare 

76,77,78,79,80,81 

26 

Low Desert 

Upland Desert 
Shrub-Blackbrush 

43,48,75 

27 

High Desert 

Sage-Grass 

101,115,116,117 







2. 1.4. 2 — Continued. 


coordinates to a list of Landsat coordinates. The Landsat 
coordinates were input to the IDIMS function GRNDSPOT to 
determine the Landsat spectral class occurring at each photo 
plot location. The resulting list of spectral class by photo 
plot was ready to be paired with the photo interpretation (PI) 
data existing for each plot. 

The following data from the 1978 LSP PI was input to 
the IDIMS Earth Resources Inventory System (ERIS) for analysis 
by summary class: 

1. All Summary Classes 

a. % cover - trees 

b. % cover - shrubs 

c. % cover - grasses/forbs 

d. % cover - nonvegetation 

2. Pinyon-Juniper Summary Classes 

a. % crown closure 

b. number of trees 

c. total tree height 

3. Ponderosa Pine Summary Classes 

a. % crown closure 

b. number of trees 

c. total tree height. 

The list of spectral class by photo plot was restructured in 
ERIS into a file of photo plots by summary class (see Table 2-4 
for spectral class by summary class) and this was in turn used 
to restructure the PI data-by-photo plot listings into a merged 
file of PI data by photo plot by summary class. ERIS was then 
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used to perforin the analysis of variance (ANOVA) on the 
merged-f ile . Table 2-5 is an example of the results of 
these ANOVA runs. The column labeled "level" indicates 
summary class number. N is the number of photo plots repre- 
senting each summary class (note: the 1978 LSP could not be 
allocated with an intent of adequately sampling the summary 
classes since this Landsat processing phase occurred after 
the LSP acquisition). The "mean" and "st. dev." (standard 
deviation) are calculated from the total PI observations in 
each summary class, e.g., 256 observations for summary class 15. 
The mean value for each PI parameter listed above was used in 
the productivity estimation sample allocation procedures as well 
as covariance matrices calculated between PI parameters for 
each summary class. Table 2-6 shows an example of the 
covariance matrices calculated from the combined Landsat data 
and 1978 LSP PI data. 

2. 1.4. 3 Combining Allotment Boundaries and Summary Classes 

For Sampling . 

After the summary classes were described in terms 
of mean and covariance for % cover, number of trees and tree 
height, the allotments specified by BLM (the area to be 
sampled) had to be combined with the classification to 
determine the areal extent of each summary class. This combined 
data set, then, represented the Landsat-based stratification 
input to the sample allocation procedures (Section 2.2.2). 

After the boundaries of the allotments to be sampled were 
digitized from maps (Section 2.1.5), these boundaries were used 
to masJc the Arizona Test Site classification such that only 
the classified areas of interest remained. The allotments 
were grouped into range, woodland and foreet strata so that 
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Table 2-5. Example of the 1978 PI ANOVA Results 
For Percent Cover of Trees. 

NOTE TRCC GHYi 

ONCWAX ANOVA FOR TRCC SUBaCRlRTC IN GRCVI 


ANALYRia OP VARIANCE 


buc TO 

OF 

fa Nfsaa/DF f*ratio 

FACTOR 

23 

112931, 

4992,66 49.9092 

ERROR 

1192 

119216, 

100,960 

TOTAL 

1209 

231747, 



PHOTO 



SWIMARY CLASS 

PLOTS 



LEVEL 

N 

MEAN j 

IT, OiV, 

1 

If 

,0000001400 

•000000C400 

2 

79 

•126991E-01 

.112909 

3 

21 

,000000C400 

.000000E400 

4 

2 

,000000c400 

.OOOOOOC400 

9 

9 

•000000E400 

.192999C-04 

A 

20 

,990001 

2,49967 

• 

14 

,000000c400 

,169261C*04 

9 

12 

3,33333 

9,19919 

10 

69 

,304349 

1,94992 

11 

99 

,724491 

2,14903 

12 

199 

3,20649 

6,62912 

13 

14 

,0000ti0C400 

,16929tC*04 

14 

107 

7,91402 

10,9159 

19 

296 

17,1690 

14,2919 

16 

127 

24.0630 

19,1194 

17 

9 

21,2222 

13,4090 

11 

9 

11,7779 

16.9616 

20 

2 

47,0000 

4,24264 

21 

31 

27,1949 

14,7412 

22 

9 

44,aooo 

21,9002 

23 

9 

19,6000 

10,3344 

29 

4 

,O00OOOE4OO 

,ooooooc#oo 

26 

39 

4,17143 

9,19931 

27 

24 

2,09333 

4,40273 


POOLED ST, DCV, ■ t0,0429 


srsssw 
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Table 2-6. Example of Covariance Matrices 
From 1978 PI Data Analysis. 

>KCTR 

CNTCH INPUT PILE NANC ••ASPZGf 
LIST? 

CNTCA CAie tCLCCTXON CAXTCAXA 
**GACYlel 

CNTCA VAPXABLC NAMES ONE PEN LINE (ARRAY SUBSCRIPTS OPTXONA^J 


IS COLUMNS AND 81 CASES LOADED 
>COV TREE SHRUB CP NVEG 



TREE 

SHRUB 

GP 

NVCG 

TREE 

• 000 




SHRUB 

.000 

74,251 



GP 

,C00 

•16,472 

111,265 


NVEG 

• 000 

-J7,77S 

•274,791 

}12,571 

>COV PP 

PPNT PPHT 





PP 

PPNT 

PPHT 


PP 

,000 




PPN? 

• 000 

,000 



PPHT 

,000 

.000 

,000 


>COV PJ 

PJNT PJHT 




• 

PJ 

PJNT 

PJHT 


PJ 

,000 




PJNT 

•ooo 

,000 


• 

PJHT 

,000 

,000 

,000 
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sampling could be performed within each stratum separately. 
Figure 2-7 illustrates these allotment boundaries with 
respect to the project area for the three strata. Within each 
stratum, the total number of pixels assigned to each summary 
class was determined using the IDIMS function PIXCOUNT (see 
Table 2-7 ) . 


At this point in the Landsat processing, each summary 
class was described in terms of its mean % vegetative cover 
by type, the covariance between those types, and its total 
coverage (number of pixels) within each of the range, woodland 
and forest strata. However, two remaining Landsat processing 
tasks. Digitizing and Digital Terrain Data, will be covered 
prior to the discussion of the use of the data set above in 
data collection procedures. 

2.1.5 Digitizing . 

This aspect of thfi Landsat processing facilitated the 
use of map-based information in both the classification phase 
and the sampling/estimation phase. Environmental strata and 
administrative boundaries were digitized from a number of maps 
and registered within the same geoblock containing the Landsat 
data. It is important to note that the 71 control points used 
during the correction of continuous geometric errors (Section 
2.1.1. 3.2) were also used in overlaying the digitized data. The 
following sections describe the data digitized and the objectives 
and applications of the digitizing results. 
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Figure 2-7 (a). Range Allotment Boundaries. 
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Figure 2— 7 (b) . Woodland Allotment Boundaries. 
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Figure 2-7 (c). Map Depicting Forest Areas A and B 

Within the Arizona Test Site. 
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Table 2-7 (a). Pixel Counts for Range Stratum. 


Pixel 

PANGC 


COUNT 


TOTAL COUNTEn* 4300000 NOT COUNTED* 0 

CLASS COUNT PERCENT 


0 

0 

,00000 

1 

146990 

.03409 

3 

41670 

,00969 

3 

28902 

,00663 

4 

282 

,00007 

5 

10334 

,00240 

6 

29105 

,00677 

7 

439 

,00010 

9 

9689 

,00225 

9 

1973 

,00046 

10 

63301 

,01472 

n 

53893 

,01253 

12 

114149 

,02655 

13 

1413 

,90033 

14 

80770 

,01878 

19 

97882 

,02376 

1 A 

74074 

,01723 

17 

3908 

,00091 

18 

4689 

,00109 

19 

588 

,00014 

20 

263 

,00006 

21 

4753 

,00111 

32 

141 

,00003 

33 

3344 

,00078 

24 

262 

,00006 

39 

2877 

,00067 

36 

33377 

,00776 

27 

32746 

,00529 


2-41 



Table 2-7 (b). Pixel Counts for Woodland Stratum 


ORIGINAL PAGE 53 
OF POOR QilALITY 


P f X e L 
*OOOLAM^ 

TOTXL COUNTCOi 


COUNT 

455000 


NOT COUNTtni 


t) 

0 

,00000 

1 

6800 

,01495 

2 

4968 

,01092 

} 

264 

,00058 

4 

0 

,00000 

5 

650 

,00143 

6 

die 

,00194 

7 

26 

,00006 

9 

8 

,00002 

9 

146 

,00032 

10 

178 

,00039 

n 

651 

,00143 

12 

5u95 

,0l 120 

1 1 

2 

,00000 

1 4 

10703 

,02352 

15 

29395 

,06460 

15 

25275 

,05555 

17 

1296 

,00293 

19 

274 

,00060 

19 

0 

,00000 

20 

114 

,00025 

21 

1216 

,00267 

22 

32 

,00007 

21 

362 

,00080 

24 

16 

,00004 

25 

5 

,00001 

76 

4928 

,01061 

27 

3917 

,00861 


0 
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original 
OF POOR 

Table 2-7 (c). Pixel Counts for Forest Stratum. 


page ts 
quality 


PIXEL COUNT 

rOXBST 

TOTAL COUNTEO* IbOOOOO NOT COUNTED* 0 


Class count percent 


0 

0 

•00000 

1 

42AS1 

•01190 

2 

40209 

• 01117 

1 

196S0 

•00546 

4 

756 

•00021 

5 

56B6 

•00159 

6 

16953 

•00469 

7 

399 

*00011 

9 

393 

• OOUl 1 

9 

3213 

•00099 

10 

4690 

•OOllO 

11 

11493 

•00319 

12 

4007b 

,01113 

11 

154 

• 00004 

14 

54260 

•01507 

IS 

155996 

•04330 

16 

1 19941 

•03304 

17 

7979 

•UU222 

lA 

6013 

,00167 

19 

249 

,00007 


20 

1460 

,00041 

21 

19061 

,00502 

22 

7373 

,00205 

23 

4294 

,00119 

24 

4 

,00000 

25 

166 

,00005 

26 

32477 

,00902 

27 

19996 

,00525 
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2. 1.5.1 AcJministrative Boundaries. 


There were two general classes of administrative 
boundaries digitized; "ownership” and management. The project 
area and all lands within the project area administered by BLM, 
the state and those privately owned, were delineated on maps 
and represented "ownership." Within the project boundary, a 
number of allotments and pastures within those allotments were 
specified by BLM as areas for range and woodland productivity 
estimation and therefore represented management boundaries. 
Additionally, two intensive mapping areas, A and B, were identified 
by BLM for forest productivity estimation and their boundaries were 
also digitized. The project area was digitized from a series of 
uses 1:250,000 scale paper map sheets. Administrative boundaries 
were digitized from a series of 1:63,360 scale mylar copies produced 
from the general highway map set for Mojave County (1969) as 
originally published by the Arizona Highway Department. Figure 2-8 
illustrates the allotments and pastures digitized within the 
project boundaries (ownership patterns are not shown to minimize 
confusion; see bolov.) . Figure 2-9 illustrates the intensive 
mapping areas A and B. 

The digitized administrative boundaries (pastures 
and allotments) were converted to image masks and overlayed 
onto the Landsat classification results for use in controlling 
subsequent sampling and data summarization efforts. Even though 
the productivity estimates were to be for each allotment and 
pasture for all ownerships, no field plots were to be taken on 
private lands. As a result, the numerous scattered small parcels 
of private ownership within the project area were digitized 
to control only the placement of samples and not the data 
summarization. Section 2.2, Data Collection, describes the use 
of the administrative boundaries in sample allocation and 
selection. 





Figure 2-8. Allotment and Pasture Boundaries Digitized 

Within Arizona Test Site. 
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Figure 2-9. Map Depicting Intensive Mapping Areas A and B 
Within the Arizona Test Site. 
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2.1. 5. 2 Environmental Boundaries. 


As discussed in Section 2. 1.4.1, an environmental 
stratification based on elevation was used to eliminate con- 
fusions in the classification prior to sampling. This 
environmental stratification was developed by delineating 
specific contour lines and special areas on 1:250,000 scale USGS maps 
covering the project area. BLM-Arizona personnel produced the 
input map contining the 4000 foot and 6000 foot contour lines and 
also areas known to contain Ponderosa pine. The resulting 
digitization of this map was converted to an image mask and 
used to control the reassignment of spectral class numbers. 

2.1.6 Digital Terrain Data (DTD) . 

The detailed descriptions of spectral classes in 
terms of percent vegetation composition were to be augmented 
by topographic information. Specifically, the ranges of 
elevation, slope and aspect for each class would be determined 
and reported to further explain class differences based on 
topographic influences. Defense Mapping Agency (DMA) digital 
terrain data (DTD) tapes distributed by the National Cartographic 
Information Center (NCIC) were used to provide the elevation 
data while slope data and aspect data were interpolated from 
the elevation data. 

DTD tapes are organized by one degree blocks of latitude 
and longitude and correspond to the east and to the west halves 
of the 1:250,000 scale USGS map series. Portions of four of 
these blocks were required to cover the Arizona Test Site 
(see Figure 2-10). 
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2 . 1 . 6.1 


DTD Reformatting . 


The input for the IDIMS function TERRAIN is the 
DTD tape disci' ibuted by NCIC. The terriin data is recorded 
on 9-track 1600 bpi with IBM's variable block physical record 
format {block size = 32,756 bytes). The output of the TERRAIN 
function is a digital image whose intensity values are in 
integer format. Each line of this image represents a south to 
north quantized elevation profile with a resolution of one foot 
in elevation. Since "neat ' unds" of the digitizing is a 
latitude or longitude line (not orthogonal on a transverse 
mercator projection) , fill is required to create the rectangular 
image in IDIMS. As the DTD images are entered into IDIMS, a 
report is produced from the header section of each image listing 
information about the map sheet from which the data was recorded. 
Also included in this report are the total second values for the 
four corner points that were recorded prior to digitizing the 
actual elevation data and the resulting line and sample values 
in the Dtu image for the four points. 

2. 1.6. 2 DTD-to-Landsat Registration Procedures . 

As described previously, only the four corner points 
of a DTD image have values in latitude/longitude. To install 
new points, the 15 minute tic intersections within each 1® block 
for the Arizona Test Site were computed by adding increments 
of 900 seconds to the latitude and longitude values of the 
southeast corner of that block. These points, described 
by coordinates of total secc i ilr . were converted by the program 
ALLCOORD into UTM coordinates for UTM zone 12. A total of 25 
points for each 1® block were established in this fashion. 

The four corner coordinates for each 1® block in terms 
of IDIMS DTD image line and sample values were determined and 



2.1.6. 2 — Continued. 


paired with their corresponding UTM values from the step above. 
The piogram TRNSFORM was used to create a transformation from 
the UTM coordinates to the DTD images. A first order least 
squares fit was all that was needed since both the source (UTM) 
and destination (DTD image) coordinate systems were rectangular. 
The transformation equations were then used to locate the 
DTD image line and sample values for the 25 UTM points for 
each 1° block. 

The 25 UTM points for each block were next used in 
the ALLCOORD program to determine their GES internal coordinate 
system values for subsequent input to the GES-to-50 meter (T2) 
transformation equations. This resulted in the location of the 
25 UTM points and their coordinates relative to the subsectioned 
50 meter UTM grid, the geometrically corrected classification 
image space. These values were then used as a source in 
TRNSFORM with the 25 DTD image points as destination to create 
first order least squares transformation equations. The 
function REGISTER was used to warp the DTD images to the 
geometrically corrected classification image. The registration 
of the Las Vegas (East) 1“ block was accomplished by extending 
the UTM zone beyond the 144* longitude boundary. 
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2 . 2.1 


Overview of Sample Design . 

2. 2. 1.1 General Description . 

The sample design for each frame (woodland, forest and 
rangeland) included a two-stage sample, where the first stage 
primary sample units (PSU's) were defined by the 50 meter UTM 
grid and the second stage was comprised of a double sample of 
photo and ground plots (secondary sample units, or SSU's) from 
within each selected PSU (see Appendix 2-A) . The sampling within 
each frame was done independently of the others, although the PSU's 
for each frame were defined by a common rectangular grid overlay 
of 9x40 pixels (450 x 2000 meters)* and there was some overlap between 
frames. Other features of the sample design varied with the frame. 

2. 2. 1.2 Rangeland Design , 

The PSU grid overlay of 9 lines by 40 samples resulted 
in a population size of 3032 PSU's covering 215,276 hectares 
(531,957 acres). A compilation of the 27 Landsat summary classes 
for each PSU 'was produced using the IDIMS function SAMPLET. 

Then the S’^iumary classes (or in some instances combinations of 
classes) were used to stratify the PSU's. It was originally 
planned to use a plurality rule. However, it was decided that 
this would seriously bias the sampling because the scattered 
occurrences of a given class may be characteristically different 
than the instances where there are dense concentrations of that 
class. Also, some classes may have few if any PSU's under the 
plurality rule. Therefore, the strata were formed by the 
existence of any pixels of the given class in the PSU. This 
resulted in overlapping strata in that a particular PSU might be 
selected from any of several strata. Equal probability sampling 
of PSU's was used from within the stratum representing a class 
regardless of the number of pixels of that class in the PSU. 

* Note: This PSU size was used in actual sample selection-location 
efforts. A PSU size of 8x40 pixels was used in the Survey Planning 
Model that determined the number of samples necessary to achieve 
the accuracy and precision design goals. 
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2 . 2 . 1.2 


Continued . 


The number of PSU's to sample for each stratum was determined 
from the U.C. Berkeley Survey Planning Model (Section 2.2.2). 

The population and sample sizes by Landsat PSU strata are summarized 
in Table 2-8. 

The 108 sampled PSU's were used as locations for flight 
lines, each comprising 15 large scale photo plots. These are 
the sampled SSU's. Photo interpretation was performed for 
stratifying all 1620 plots into one of the categories corresponding 
to the Landsat summary classes (Section 2. 2. 3.1). The number of photo 
plots occurring in each of the photo strata was also input into the 
Survey Planning Model which determined the number of ground plots 
to sample. This information is summarized in Table 2-9. In 
addition to the photo stratification, interpretation was made for 
percent cover, height class and foliar density by species (Section 
2 . 2 . 3 . 2 ) , 


The photo plots were sampled for ground visits in the 
following manner. Within a photo stratum a plot was selected 
by simple random sampling. Then an additional plot was selected 
by simple random sampling of the remaining plots of the same 
flight line regardless of the photo stratum. These two steps were 
repeated until the required sample size was obtained. Subsampling 
of the ground plots was also performed using subplots located 
along transect lines (Section 2. 2. 4.1). 
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Table 2-B. 


Summary of Landsat Stratification for Rangeland PSU's 


Landsat Summary 
Class No. 

Total No. * 
of PSU's 

Sample Size 
Prom SPM 

Actual Sampl' 
Size Chosen 

1 

1249 

9 

15 

2 

748 

2 

4 

3 

1052 

2 

4 

6 

803 

2 

4 

8 

464 

2 

4 

9/12/27 

1786 

9 

15 

10/13 

971 

5 

8 

11 

1275 

4 

8 

14 

2036 

5 

8 

15 

1802 

6 

10 

16 

1470 

5 

8 

17/18 

524 

2 

4 

21 

244 

2 

4 

23 

375 

2 

4 

26 

1195 

4 

8 

All Classes 

3032 

61 

108 


* The total number of PSU's for each summary class was derived 
by tallying all PSU's that contained at least one pixel of 
that class; therefore, PSU's could be counted more than once 

** The sample size from the SPM was increased for each class to 
the value shown in the "actual" sample size column because 
budget existed to allow more samples. The SPM value is the 
sample size expected to achieve the sample design accuracy 
goals at least cost . 
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Table 2-9. 


Summary of Photo Stratification for Rangeland SSU's 


Photo Stratum 

No. of 

Photo Plots* 

Actual Ground 
Sample Size 

1 

171 

15 

2 

80 

10 

3 

26 

6 

4 

106 

6 

5 

14 

0 

6 

112 

10 

7 

2 

0 

8 

90 

7 

9/12/27 

230 

13 

10/13 

8 

3 

11 

189 

13 

14 

16S 

11 

15 

163 

14 

16 

141 

14 

17/18 

17 

4 

19 

1 

0 

20 

0 

0 

21 

4 

2 

22 

0 

0 

23 

16 

2 

24 

7 

0 

25 

41 

1 

26 

32 

4 

Uninterpretable 

2 

— 

Totals: 

1620 

135 


* 108 PSU's for the Rangeland with 15 SSU's (photo plots) per 
PSU totals 1620. This column reflects the sorting by a 
photo interpreter of the plots into categories related to 
the Landsat summary classes. 
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2.2. 1.3 Woodland Design . 


The same PSU grid was used for the woodland population 
resulting in 388 PSU's covering 26,449 hectares (65,358 acres). 
However, the Landsat summary classes were used in a different 
manner for PSU selection. No stratification on the classes was 
used — rather the sum of the pixel counts for classes 14, 15, 

16 and 17 was used as sampling weight! for unequal probabilities 
of selection for FSU's. Forty-one of the PSU’s covering 2,187 
hectares (5,405 acres) had no pixels in any of those classes 
and therefore could not be sampled. The sampling frame therefore 
had 347 PSU’s covering 24,263 hectares (59,953 acres). 

A total of 45 PSU’s were selected with replacement 
as flight lines. Each flight line had 15 photo plots for a total 
of 675. Again the photo stratification was performed and the 
total number of plots per stratum was compiled. The SPM was 
again used to determine the number of ground plots to sample 
(Section 2.2.2), but the only stratification used in selecting 
plots for ground visits was one which separated pinyon juniper 
types from all others. In addition to the photo stratification, 
interpretation was performed for percent cover, height class 
and foliar density class by species. Detailed photo interpretation 
was performed for a subset of trees on each plot (Section 2. 2. 3. 2. 2). 
Subsampling of the ground plots was performed in order to collect 
ground measurements on the same trees that were photo interpreted 
(Section 2 . 2 . 4 . 2 ) . 



2-55 


?. .2.1.4 Forest Design . 


The PSD grid overlayed on the forest population 
resulted in 2,063 PSU's covering 153,114 hectares (378,351 acres). 
The Landsat sunuiiary class counts were again used as a stratifica- 
tion for selection of PSU's for flight lines. Tha existence of 
any pixels in class 20, 21 or 22 (all Ponderosa pine types) 
qualified a PSU for inclusion into the correspondixig stratum. 

This reduced the sampling frame to 441 PSU's covering 35,214 
hectares (87,015 acres). The SPM was used to determine the required 
number of flight lines to be allocated to the three strata 
(Section 2.2.2). These results are summariicied in Table 2-10. 

Table 2-10. Summary of Landsat stratification 
for forest PSU's. 

Landsat Total PSU Sample Actual PSU 

Class No. No. of PSU's Size from SPM Sample Siz e * 


20 

101 

2 

4 

21 

‘’39 

15 

25 

22 

108 

10 

17 

Classes 

441 

27 

46 


* Increased to take advantage of available budget. 

The PSU's were selected with replacement and with equal probability 
of selection. Each of the 46 flight lines contained 15 plots for 
a total of 690 plots (originally there was a 47th flight line, but 
it was dropped because we did not receive an acceptable set of 
photography) . 

Stratification of the photos was performed for all 
plots. There were 56 plots in class 20, 62 in class 21 and 
47 in class 22. The remaining 525 photo plots were in non- 
Ponderosa pine types. Ground plots were selected by the same 
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2 . 2 . 1 . 4 


— Continued. 


method used for rangeland, 2 plots from class 20, 15 from class 21 
and 10 from class 22 for a total of 27 ground plots from 46 PSU's. 
Photo interpretation for percent cover, height class and foliar 
density by species was also performed. Detailed interpretation 
was made on a subset of trees on each plot (Section 2. 2. 3. 2. 2). 

The ground plots were subsampled to obtain measurements for the 
same trees as interpreted on the photos (Section 2. 2. 4. 2). 
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2 . 2.2 


Allocation and Selection of Primary Sample Units (PSU's) . 


2 . 2 . 2 . 1 U.C. Berkeley Survey Planning Model Allocation 

Optimization Efforts . 

The UCB Survey Planning Model (SPM) was used to 
determine sample sizes for estimating cattle-usable forage 
biomass in the Rangeland stratum and wood volume in the Forest 
and Woodland strata. The SPM utilizes a nonlinear programming 
procedure for finding the number of samples to be allocated 
between substrata and sample stages (PSU, SSU, ground) that 
will minimize total variable cost subject to meeting prior 
constraints on sampling error for the parameter (s) to be 
estimated. A full description of the SPM and its application 
was given in the NASA/BLM ASVT Phase I report (ESL Incorporated 
1979) . In addition, a summary of the SPM is given in 
Appendix 2-B. 

A general description of the stratified two stage 
with double sampling design has been given earlier (Section 
2.2.1). Specific assumptions used to compute sample size with 
the SPM are given in Table 2-11. These assumptions were based 

on initial design work done for the Arizona test site in 
Phase I by UCB (see ESL 1979) and modified in subsequent Phase II 
discussions between UCB, ESL, Resource Inventory Services, BLM and 
NASA. Note in particular that the 15 SSU per PSU size was set based on 
(1) results from Phase I work for Arizona showing SSU sample 
requirements generally did not exceed 15 and (2) on practical/ 
cost considerations concerning the number of stereo pairs to 
be obtained per flight line. The 2 km by 400 m PSU size was set 
to approximate the rectangular sample units reported in English 
measure during Phase I. 
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Table 2-11. Assumptions Used With the U.C. Berkeley Survey Planning 
Model Regarding the Two Stage With Double Sampling 
Design to be Used In the 1979 Arizona Test Site 
Inventory. 


Design Characteristic 


Assumption 


I. Fixed 


A) Sampling substrata 


B) Maximum number of sampling 
substrata 

C) PSU size 


D) PSU's assigned to sampling 
substrata based on 


E) SSU layout 


F) Large scale photo SSU 
selection probabilities 

G) Ground SSU selection 
probabilities 

H) Correlation between ground 
parameters of Interest (e.g,* 
forage biomass and wood 
volume) and substitute 
parameters (e.g., percent cover, 
height, number of individuals) 
actually simulated in the SPM 

I) Allowable sampling error 

1) Cattle - usable forage 
biomass 

2) Forest (ponderosa pine) 
volume 

3) Woodland (pinyon and 
juniper combined) volume 


Composed of Landsat summary classes, 
individually or in combination 

15 (a SPM limitation as of spring 
1979) 

40 cells horizontally by 8 cells 
vertically (2 km by 400 m) 

Plurality rule: i.e., PSU assigned 
to sampling substratum having the 
most number of pixels in that PSU 

15 systematically-spaced SSU's along 
the long axis of a PSU 

Equal within a PSU 


Equal among all SSU's photo-sampled; 
otherwise zero probability of selection 

High (effectively r=1.0) 


+ 20 % 0 805^ confidence level over the 
entire rangeland stratum 

+20% @ 80% confidence level over the 
entire forest stratum 

+^20% 0 80% confidence level over the 
entire woodland stratum 
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Table 2-ll. “-Continued. 


Design Characteristic 
J) Sample size 

1 ) PSU's/substratum (h) 

2) Photo SSU's/PSU (i) selected 

3) Ground SSU 's/substratum in 
PSU's selected 

II. Not Fixed 

A) PSU selection probability 


B) Large scale photo-to-ground 
correlation 

C) Data acquisition/measurement 
costs 


Assumption 

Min. allowable < size _< max. (total 
number available) 

1 < nfi < 

1 1 < 15 

3 £ n^ < 15*njlj 

Equal or variable (i.e,, pss) within 
sampling substrata 

Variable (see Table 2-15) 

To be defined (see Table 2-14 ) 


2-60 


2 . 2 . 2.1 


— Continued. 


Attention should also be drawn to the inventory 
sampling error goals. They were +20% at the 80% 
confidence level for cattle-usable forage biomass in the 
Rangeland stratum and similarly for Ponderosa pine wood 
volume and Pinyon-Juniper wood volume in the Forest and 
Woodland strata, respectively. In other v/ords, the true 
value for the parameter to be estimated (e.g., forage 
biomass) was to fall within +20% of the estimated toual, 

80 times out of 100 under the following assumptions: 

(1) sample plot values for the random variable representing 
biomass or volume tend to be normally distributed within 
sampling substrata, (2) measurements of biomass or volume 
are unbiased, (3) the estimator used to link sample 
stage data together are unbiased, and (4) the estimates of 
resource parameter means and covariances used in the SPM 
represent faithfully the relative magnitude and variability 
of the ground parameters of interest (e.g., biomass and 
volume) . 


Since this is a first-time inventory the fourth 
assumption cited above must also be treated carefully. 
With no reliable ground forage biomass or wood volume 
data for all sampling substrata, substitute resource 
parameter information had to be used instead in the SPM. 
Percent composition of trees, shrubs, and herbaceous 
materials obtained from interpretation of large scale 
aerial photography (LSP) obtained during 1978 were em- 
ployed as forage biomass “correlates" in SPM runs. If 
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2 . 2 . 2.1 


— Continued. 


these percent composition data were not strongly correlated 
with palatable forage biomass magnitude and variability, then 
sample allocation results would be less efficient (larger cost 
for fixed error) than if forage biomass was used directly. 
Similarly, LSP-derived tree height class, number, and percent 
composition datt were substituted for volume in the Forest and 
Woodland strata. This source of sample inefficiency can be 
removed in subsequent inventories by using the forage biomass 
and volume data resulting from the first inventory directly in 
the SPM. 


2.2.2. 1.1 Inputs . 

2. 2. 2. 1.1.1 Mean Vectors and Covariance Matrices from 1978 
PI Data for Landsat~~Summary Classes . 

The first data input required for Survey Planning 
Model computer runs consisted of a specification of sampling 
substrata and the mean resource parameter values and variance- 
covariance matrix associated with each substratum. 

Sampling substrata in the Arizona Test Site were 
composed of Landsat summary classes, either individually or in 
combination. The maximum number of sampling substrata within 
a given stratum (Rangeland, Forest, or Woodland) was limited to 
15. This limitation was imposed by the SPM capabilities as 
of March 1979, but was also consistent with the bi 'get available 
for within substratum LSP and ground plot sampling. As a 
consequence, the 27 Landsat summary classes in the Rangeland 
stratum were grouped into 15 sampling substrata. Table 2-12 
shows which classes were combined and which were eliminated from 
consideration. Classes having similar species composition and 
thought to be equally forage-productive were combined, while classes 
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Table 2-12. Original Landsat Classes Deleted or Combined 
in Rangeland Stratum 


Summary 
Class No. 

Deleted 

Combined 

Reason 

4 

/ 


Low (<2 percent) grass and forb 
cover and very small area 

5 

/ 


Very low grass/forb cover; this class 
represented blackbrush in rocky areas 

7 

/ 


Negligible forage and small area 

19 

/ 


No rangeland forage 

20 

/ 


Low forage cover and small area 

22 

/ 


Low forage cover and small area 

24 

/ 


No forage 

25 

/ 


Low forage cover 

17,18 


/ 

Similar type: mountain shrub/mixed 
chaparral 

9,12,27 


/ 

Similar type: combination grass/ 
shrub/sagebrush 

10,13 


/ 

Similar type: combination snakeweed- 
grass/salt shrub 
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2. 2. 2. 1.1.1 — Continued. 


were eliminated if the area represented was small and/or if 
little cattle-usable forage was available. All other Landsat 
summary classes not listed in Table 2-12 were treated as 
separate sampling substrata. 

In the SPH, sampling rates were determined by minimiz- 
ing cost subject to constraints on sampling variance for the 
random variables representing resource parameters of inventory 
interest. Nominally, the resource parameters on which basis 
sample intensity would be controlled in the 1979 Arizona Test 
Site inventory were (a) cattle-usable forage biomass in the 
Rangeland stratum and (b) wood volume for Ponderosa pine 
and Pinyon-Juniper in the Forest and Woodland strata, respectively. 
However, since this was a first-time inventory, data on the 
magnitude and spatial variability of these resource parameters 
did not exist by sampling substratum. Thus the substitute 
parameters listed in Table 2-13 were used instead. The 
assumption to be made was that these substitutes, when taken as 
a group, faithfully represented the relative spatial variability 
of the original parameters. Substratum-specific estimates of 
means and covariances for these substitute parameters were 
generated from photo interpretation of the large scale aerial 
photography acquired over the Arizona Shivwits? Unit by BLM in 
1978. Mean vector and covariance matrix estimates for substrata 
consisting of combinations of Landsat classes are given in 
Table 2-13. The mean and covariance data were used in the SPM 
to simulate substitute parameter values for all 50m x 50m 
Landsat class map cells. When cells were aggregated into 
primary sampling units (PSU's), between and within PSU 
covariances could be computed for each sampling substratum. 

These covariance data were in turn used as a measure of variability 
in the SPM nonlinear programming module to calculate sample 
sizes designed to meet sampling error goals. 
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Table 2-13. Substitute Parameters Simulated in the SPM and 

Corresponding Mean Vector and Covariance Matrices 
for Landsat Summary Classes That Were Combined 
Into Sampling Substrata. 


Substitute Parameters: 


rangeland 


Forest 


Woodland 


1 ) Percent cover 
for total trees 

2) Percent cover 
for total 
shrubs 

3) Percent cover 
for total grass 
and forbs 


substi- 
tute for 
forage 
biomass 


1) Avg. height 
class for 
ponderosa 

2) No. of stems 
for ponderosa 


3) Percent cover 
for ponderosa 


substi- 

tute 

'for 

wood 

volume 


1) Avg. height for 
pinyon & juniper 

2) No. of stems for 
pinyon & juniper 
combined 

3) Percent cover for 
pinyon & juniper 
combined 


Mean Vector and Upper Triangle of Covariance Matrices for Combined 
Landsat Summary Classes in the Rangeland Stratum: 



Classes 17,18 

Classes 9,12,27 

j 

Classes 10,13 


Means 

Means 

Means 

Trees 

16.5 

3.1 

0.3 

Shrubs 

28.6 

11.8 

6.4 

Grass/ 

Forbs 

4.9 

15.8 

14.5 

Covariance Matrix 

Covariance Matrix 

Covariance Matrix 

259.6 

-353.9 -38,8 

39.7 -12.4 -34.3 

3.4 1.0 -3.0 


1096.7 -197.2 

111.5 -46.7 

29.3 -26.7 


230.7 

336.4 

207.7 
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2. 2. 2. 1.1. 2 Expected Correlations Between Photo and Ground . 


Another input to the SPM required for the double 
sampling phase of the inventory design was the expected corre- 
lation between variables measured on the LSP and the parameters 
of interest (e.g., forage biomass) on corresponding ground 
plots. In the Ra igeland stratum, no conclusive evidence was 
available to establish within sampling strata correlations 
or aven regional correlations. Work done by Cohen (1979) in 
Nevada suggested that correlations between percent composition 
of big sagebrush (A. tridentata) taken from LSP (70mm, scales 
approximately 1:500-1:1000) versus ground forage biomass were 
in the .8 to .9 range. However, this vegetation type was 
somewhat different than the sagebrush-perennial grass type in 
Arizona. In addition, no data existed for other vegetation 
types. Consequently, Rangeland correlations ranging between 
.6 and .85 were used in the SPM. 

More conclusive evidence regarding correlation between 
LSP measured variables and wood volume existed for the Forest 
and Woodland strata. A number of studies have shown 
correlations between .8 and .9 for LSP height, count, crown area, 
and percent cover data in the Ponderosa pine type. Similar 
correlations between weight and aerial photo-derived measurements 
of crown area and crown diameter have been documented for 
Pinyon-Juniper types by Meeuwig et al. (1979 ). In order to 
determine the sensitivity of sample size results to correla- 
tion, photo-to-ground correlations for wood volume were set to 
.7, .8, ,85 and .9 in the SPM. A given correlation was 
assumed to hold for all sampling substrata simultaneously. 
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2. 2. 2. 1.1. 3 Expected Coefficients of Variation . 


Coefficients of variation (standard deviation divided 
by mean value) were computed for wood volume and number of 
trees by Resource Inventory Services (San Jose) for ESL 
Incorporated. These CV's were determined for varying lengths 
of ground sample transect using the 1978 LSP PI data and a 
technique for simulating tree populations and samples thereof 
(see Appendix 2-C) , A wood volume CV versus cost analysis 
performed by Resource Inventory Services established a 150 foot 
ground transect as the most efficient ground sample layout. 
Coefficient of variation data on number of trees per 150 foot 
transect was then provided to UCB for use in the sample alloca- 
tion work. These data were (a) Ponderosa Pine Forest volume: 

CV = .65 - .70 with expected 2.4 to 3.6 trees per transect; 
and (b) Pinyon-Juniper volume: CV = .40 - .45 with expected 
4 to 7 trees per transect. The coefficient of variation data 
was then used to compute ground and LSP measurement costs for 
use in the SPM as described in the next section. 

2 . 2 . 2 . 1 . 1 . 4 Expected Costs Itemized for Landsat, LSP 

&;quisition/'Interpretation, and Ground Data 
Collection . " 


The cost function minimized during calculation of 
the optimal sample size distribution in the SPM was 


'TVC 


L 

= E 

h=l 


'°hl"h ^ 'h2"h"h"h ^ =h3 


+ ‘=h4"S"h + ■=h5'^ 


for the two stage with double sampling design. In this equation 


Ctvc “ total variable cost, i.e,, the total of costs that 
increase as a function of sample size. Fixed 
overhead costs are not included. 

h = a sampling substratum index 

L = total number of sampling substrata to be included 
in a given design 
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2. 2. 2. 1.1.4 — Continued. 


C .1 = cost of selecting and locating Landsat first 
“ stage cell clusters (PSU's) in substratum h 

Ch 2 = cost of acquiring and measuring a secondary sampling 
unit (SSU) on LSP in substratum h 

c. - = cost of travel between PSU's in substratum h 
n j 

c. ^ = cost of obtaining ground data for an SSU in 
substratum h 

c-c = cost of travel between ground plots in a given 
PSU in substratum h 

n* = number of PSU's to be sampled in substratum h 

n/ = number of SSU's to be sampled with LSP within a 
" PSU selected for sampling in substratum h 

n. = number of SSU's to be ground sampled within a 
“ PSU selected for sampling in substratum h 

njt.n^ = total number of SSU's to be sampled with LSP 
in substratum h 

n*.n. = total number of SSU's to be ground sampled in 
“ substratum h. 

The cost coefficients ^h2' *^h3' *^h4 °h5 

must be specified prior to running the sample allocation module 
of the SPM. Table 2-14v presents the values of these coef- 
ficients for the Rangeland, Forest, and Woodland strata together 
with their derivation. The subscript h has been eliminated as 
costs were assumed to be similar within sampling substrata. 

Table 2-14 was based on the best available information as of 
March 1979 from ESL Incorporated, Resource Inventory Services, 
and U.C. Berkeley. An extensive footnote section at the end 
of Table 2-14 documents the assumptions made in determining 
the cost figures. 


Table 2-14. Derivation of Cost Coefficients Used in the Survey Plenning Model 
for the 1979 Arizona Test Site Inventory. 
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Subitems Included Subitem Costs Total Cost 

Coefficient Per Sample Unit Per Sample Unit Per Sample Unit 
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Subitems Included Subitem Cost' Total Cost 

Coefficient Per Sample Unit Per Sample Unit Per Sample Unit 
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(Travel Cost 
Between Ground 
Plots) 



All numbers not footnoted represent assumptions made by U.C. Berkeley based upon prior knowledge, 
past experience, or educated guess. 
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FIO Arbitrary assumption that 1-1/2 less trees will be observed by the photo interpreter per transect 
than on corresponding ground transect (see footnote Fll) due to overlapping the crowns. 
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2. 2. 2. 1.1. 4 — Continued. 


Note also that two different costs were computed for 
Rangeland LSP interpretation and Rangeland ground measurement. 
These upper and lower figures represented the range of costs 
expected to be encountered. They were used to determine the 
sensitivity of the SPM solution for sample size to change in 
cost. 


2. 2. 2. 1.2 Outputs (Optimized Sample Sizes Based on Inputs to 
the Survey Planning Model) . 


2. 2. 2. 1.2.1 Number of Sample Units to Select in the Rangeland 
Stratum . 

Given the fixed sample design assumptions cited in 
Table 2-11 and the inputs just described, the SPM was used 
to compute sample allocation to trampling substrata and sample 
stages within substrata that would meet sampling error goals. 
Sample sizes were obtained under varying assumptions about cost 
coefficients, photo to-ground correlation, and PSU selection 
probability as summarized in Table 2-15. 

For the Rangeland stratum, the topmost cost coefficient 
set was designated primary on the basis that the costs for 
photo and ground measurement were most conservative (i.e., 
highest) . Similarly, the most conservative correlation coef- 
ficient (.6) was selected as primary given the lack of supporting 
evidence for higher values within the Arizona Test Site. Table 
2-16 shows the results under the primary assumptions and given 
that PSU selection was equal probability within sampling 
substrata. The first column on the left gives the Landsat 
summary class number (s), the next column the sampling substratum 
name, and the third column lists the number of PSU's assigned 
to each substratum by the plurality rule. Moving to the right. 
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Table 2-15, Summary of Cost, Correlation, and Probability 
of PSD Selection Assumptions, Used for Two Stage 
With Double Sampling Design. 


Cost Coefficients ($) 

(Cy ^4* ^5^ Correlation PSU Selection Probability 


I. Rangeland 


(2.51, 6. 

14, 

55.19, 

120.75, 

A* -7. .75, 

Equal 

and 

Unequal 

30.19) 




.8, .83 




(2.51, 5. 

31, 

55.19, 

120.75, 

.83 

Equal 

and 

Unequal 

30.19) 








(2.51, 6. 

14, 

55.19, 

110.69, 

.83 

Equal 

and 

Unequal 

30,19) 









II. Forest 

(2.51, 3.64, 52.75, 25.00, .7, .8, jM, .9 Equal and Unequal^ 

20.13) 

I II. Woodland 

(2.51, 4.26, 52.75, 30.00, .7, .8, .9 Equal and Unequal^ 

20.13) 


Footnotes 

1 Values underlined represent primary assumptions. 

2 Unequal probabilities In the Rangeland stratum represented the number of 
cells in all 15 substrata in a given PSU divided by the total number of 
cells in the 15 classes over the whole Rangeland stratum. 

3 Unequal probabilities In the Forest stratum represented the number of cells 
In Landsat summary classes 20, 21, 22 combined divided by their total 

over the whole area. Similarly for Woodland using classes 14, 15, and 16. 
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Table 2-16. U.C. Berkeley SPM Sample Allocation Results Recommended for Implementation 
in 1979 Arizona Test Site Rangeland Forage Biomass Inventory. 3 
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Numbers in parentheses represent final number of sample units reconsnended. 

The SPM was used to determine the number of samples required for estimating productivity; however 
since the project also had the objective of quantitatively-based vegetation mapping, all 15 
SSU's per PSU were recommended for photo acquisition and interpretation, regardless of the SPM 
results given in the table. 



Footnotes (continued): 


( 



I 

I 


• 

tx> 

tH 

I 

CM 


OJ 


.*0 



(0 




3 




cr 




0/ 




• • 

05 



lA 


O) 


0) 

• 

c 



O LT> 

0) -o T- 



LO LT) 

r- 05 P 


•f 

1 +J 

O. t/5 3 


r“ 

me* 

E 3 0) T> 


•r* 

^ 0) "S' 

re x: 



O r— 

(/) 01 p *o 


re 

• • • 

xa 0) 


X) 

01 CD <0 

r— » P 


o 

■M Q. 

re o p 05 


L. 

re ** 

3 P 1— P 


C3. 

E O f“ 

P 3 C 



•t- «\j in 

O =5 Irt 3 


C 


re tn 0) o 


o 

l/) •• CM 

O. P u 


•r" 

o> u 

c c 


4-> 

O 

05 05 re 0) 


u 

E 

xz xz 


OJ 

O l> •• 

X -p </) p 



■D 0) ^ 

< 05 


'aJ 

0) in 

— 'P X3 


t/» 

QJ O) o 

>5 0 E 



i- c 

— » . 3 


■o 

»♦- -r- •> 

p to C 


et 

1 r— ^ 

iO 05 05 


3 

4-0.0—' 

r- p C t— 


O 

O E 0> 

r— C *r- re 


u 

ire « 1— 

05 05 r- 3 


o> 

to (/) CO t 

O ID P 



O' o o 

O U 


■D 

41 0) CO 

05 05 X re 


c 

im ir^ •' 

x: P 


re 

05 XI CM « 

>5p 05 



01 re o tn 

X3 c x; 


o 

■c X 

P 05 P 


•M 

o •> 

— * re 4) 


O 

ZD r“ r— O 

X .c X c 


f 

m 1— o CM 

P p re 


a. 

Q- C r- 

05 x: 




(/) O XI P 
*— 155 




r— P P 


<u 

4- 05 .C 

05 05 O 05 



. . . . . - 

U C C P 




o re 




O -D -O 05 

re* CP 




155 re 05 




o re 



■ ■ 

P X to 155 



re 

f“ re 



■P 

-o to 05 X 



re 

45 T- X 



I- 

to JZ ‘r— ^ 



•M 

re p- OlO. 



1/5 

05 in 


t/> 

re 

p p 



•P c o 

U . o 05 


f— 

re •— ',(■? 

C T3 XT 


<v 

P x: 

•r- 05 05 P 


u 

■P 4J 




1/5 •— • • 

to P -r- >t 


00 

X C 

re O r— X5 



05 O 

X «P 


>, 

C >, -r* 

p re "O 


JQ 

•»- +J +J 

s: 05 05 



r“ •— re 

a. o. to P • 



CL — • f— 

tn re p to 


X 

E •>“ 0) 

re .a p 

to X O Q. 



05 re Q. 05 



1/5 re P 

x: X c 05 p 


(/) 

X3 o 

p o P to 



o o u 

C •!- 

• ■ 

r-» 

V P 

C O P to to 

U) 

0) 

a. "u 

•f- •— re - >5- 

c 

U 

4-> OD c 

P O Z3 to 

o 


C 1— C 3 

•o re o tn >) 

•r- 

o 

<U 3 O O 

05 o 1— a. 



E P •- P 

•r- o re 

Q. 


C 4-> OI 

p 05 4- c 

E 

■ • 

or) >> U i 

•5“ core 

3 

0) 

•— -P 4) O 

O-O-r- 

l/» 

N 

to -1- r- P 

05 C r- P "O 

1/5 

*1— 

1/5 r“ 05 1 

CL re 05 C 


(/) 

re re I/I I— o 

to p X3 re 



p rep 

c x: E 

o> 

13 

=3 3 Z3 3 O 

05 O 05 3 05 

c 

00 

oo f— m (T x: 

N -r- •— C C 

•r“ 

Q. 

a. CL (X OI Q. 

•f— P r— •— 

r— 


to O 4- f— r- 

cx 



05 re Q. 

E 


. 

=5 1— P P E 



.a u -o 

tn 05 o o re 


■'w' 


CL to P 4D to 


CO 


2-7b 


2.2. 2. 1.2.1 — Continued. 


the fourth column gives the number of PSU's to be sampled per 
substratum, the fifth column the number of SSU's out of 15 to 
be sampled per PSU , and the sixth column the number of SSU's 
per PSU to be ground sampled. Numbers enclosed in brackets 
represent the final number of sample units recommended for 
sampling. This subjective increase in sample size was made to 
take full advantage of the available ground campling budget by 
increasing the number of ground plots in substrata with low 
sam.ple size. 

It should also be noted that a rounding rule of ^0.1 
has been used in all SPM results reported here. That is, since 
the SPM determines optimality in fractions of sample units, 
the results must be rounded to integer units. 0,1 was selected 
as a conservative rounding rule whereby all sample size rcsu.1 ts 
of some integer value plus at least 0.1 were rounded to the next 
highest integer. 

The sample sizes listed in Table 2-16 were those 
recommended for implementation in the Arizona inventory.* 

Choice of the equal probability of PSU selection option over 
unequal probability (see Table 2-17 for unequal probability 
selection results) was made because; (1) costs were not 
significantly different, and (2) selection probabilities on 
the average did not differ significantly enough to provide any 
advantage to maximizing the number of LSP plots falling on 
areas of interest by selecting PSU's with variable probability. 


* With the exception that all 15 SSU's per PSU were recommended 
for photo measurement, regardless of SPM results to be used 
in vegetation mapping as well. 
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Table 2-17. Auxiliary Table Showing SPM Results For Unequal Probability of PSU Selection 
Photo-to-6round Correlation of .60, All Other Assumptions as in Table 2-16. 
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2. 2. 2. 1.2.1 — Continued. 


Tables 2-18 and 2-19 give results for equal and 
unequal probability of PSU selection when the photo-to-ground 
correlation is .83. This correlation was consistent with that 
obtained by Cohen (1979) between LSP-derived (70nm\, approximate 
scale 1:500-1:1,000) big sagebrush percent cover and ground 
forage biomass in central and northern Nevada. The principal 
difference between these results and those shown in Tables 2-16 
and 2-17 is a somewhat lower PSU sampling rate (generally one 
or two PSU's) in the larger sampling substrata. In contrast, 
the number of photo SSU's per PSU is increased by one or two 
in the largest sampling substrata. The results in Tables 2-18 
and 2-19 are consistent with the mathematical formulation of 
the variance estimators, in that higher correlation should 
lead to a higher photo sampling rate, lower sampling variance, 
and hence a smaller PSU sample to reach sampling error goals. 

Ground sampling rates did not decline as expected, since these 
were already near minimal levels. 

2. 2. 2. 1.2. 2 Number of Sample Units to Select in Forest and 
Woodland Types . 

Survey Planning Model results for the Forest and 
Woodland strata are shown in Tables 2-20 and 2-21 respectively. 
Sample sizes shown there were computed for both equal and 
variable probability of PSU selection, given a large scale 
photo-to-ground volume correlation of 0.85. The probability 
for a given PSU was equal to the sum of cells in summary classes 
20, 21, 22 for the Forest stratum (14, 15, 16 for the Woodlari 
stratum) in that PSU, divided by the sum of those same classes 
throughout all PSU’s in the given stratum. The correlation of 
0.85 was recommended by ESL Incorporated as a realistically 
obtainable value and was thus used as the standard in the SPM. 
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Table 2-18. Auxiliary Table Showing SPM Results For Equal Probability of PSU Selection for Rangeland 
Photo-to-Ground Correlation of .83, All Other Assumptions as in Table 2-16. 
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Table 2-19. Auxiliary Table Showing SPM Results For Unequal Probability of PSU Selection for Rangeland 
Photo-to-Ground Correlation of .83, All Other Assumptions as in Table 2-16. 
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Table 2-20. U.C, Berkeley SPM Sample Allocation Results for Prime Cost and 

Correlation Assumptions in the 1979 Arizona Test Site Forest Inventory. 
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Table 2-21. U.C, Berkeley SPM Sample Allocation Results for Prime Cost and 

Correlation Assumptions in the 1979 Arizona Test Site Woodland Inventory. 
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2. 2. 2. 1.2. 2 — Continued. 


Forest types Ponderosa pine-oak (class 20 ) , Ponderosa 
pine-mixed shrub (class 21), and Ponderosa pine (class 22) 
were, when taken together, distributed in a fairly contiguous 
fashion. Normally under these circumstances, little difference 
would then be expected for sample allocation based on equal 
versus variable probability of PSU selection. However, examina- 
tion of Table 2-20 shows that many more PSU’s are required 
to meet the 20 percent allowable error goal when variable as 
opposed to equal probability of selection is employed. This 
phenomenon can be explained as follows. The mean values for 
the volume-substitute variables (a) percent cover of Ponderosa 
pine, (b) number of Ponderosa trees, and (c) average Ponderosa 
height, differ significantly between class 21 and the other two 
classes. Consequently, the correlation between the sum of 
cell counts in classes 20, 21, and 22 combined in a given 
PSU, and the corresponding PSU SPM-simulated totals for cover, 
number, and height will tend to be low. This, in turn, means 
that, on the average, more high covariance PSU's will occur 
when computing within and between substratum covariances. The 
net result is that higher sampling rates are required to achieve 
a given sampling error goal. A similar, though less pronounced, 
situation was responsible for the higher PSU sizes in summary 
classes 20 and 22. Thus the sampling rates for equal 
probability of PSU selection were recommended for implementation. 

Results for the Woodland stratum shown in Table 2-21 
illustrate a situation in which variable probability selection 
did decrease required sample size. In this case, there was 
much closer agreement between Landsat summary classes in percent 
Pinyon-Juniper cover, number of stems, and average height than 
in the case of Ponderosa pine. Another factor increasing the 
correlation between PSU cell counts (for classes 14, 15, 16) 
and PSU totals for percent cover, number of stems, and height 


2. 2. 2. 1.2. 2 — Continued. 


was the more dispersed nature of Pinyon-Juniper classes among 
other Landsat summary classes. That is, a higher proportion 
of PSU's had low Pinyon-Juniper cell counts than compared to 
a similar histogram of counts developed for Ponderosa pine 
in the Forest stratum. Hence, the contrast in total percent 
cover, stems, and height between low and high cell count 
PSU's was greater for Pinyon-Juniper. 

Thus the PSU sampling rate for Pinyon-Juniper using 
variable PSU selection probability was approximately half that 
required with equal probability. Normally, this sampling 
option and associated sample allocation would have been 
recommended for implementation. However, since this was a first 
time inventory, it was important to ensure an adequate photo 
and ground sample. As a result, the PSU sample sizes with equal 
probability of selection were recommended, but PSU's were to 
be actually selected with variable probability so as to minimize 
sample error. 
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2.2. 2. Selection of PSU's. 


The selection of PSU’s for sampling was made with equal 
probability of selection for range and forest areas and with 
probability proportional to size for the woodland areas. The 
PSU’s were generated by placing a 450m x 2000m (9 pixels by 
40 pixels) grid over the classified and geometrically corrected 
image using the IDIMS function SAMPLET. The output from this 
function was an ERIS file containing a count of the pixels by 
class for each PSU. From this file, lists of PSU’s by Landsat 
strata (see Table 2-8 ) were generated by weighting the 

classes making up a ntrata by 1 and weighting all other classes 
by 0. The list was then filtered to eliminate all the PSU’s with 
a weighted total of 0. These lists were then used for selection 
of PSU’s. 


The selection process, for range and forest, consisted 
of generating a list of random numbers between 1 and the number 
of PSU’s in the stratum being sampled (see Table 2-8 ). This 

list of random numbers determined which PSU's in the list were 
sampled. The procedure for the woodland samples consisted of 
calculating the cumulative sum of the weighted total for each 
PSU in the list, then generating a list of random numbers between 
1 and the total cumulative sum for all the PSU’s. The random 
numbers then determined the PSU's selected by finding tha first 
PSU in the list with a cumulative total greater than or equal 
to the random number in the list. The selected PSU’s are listed 
in Appendx 2-D. 

The selected PSU’s were then plotted on USGS 7.5’ 
and 15’ quads to assist in acquisition of the large scale 
photography (LSP) . This consisted of calculating the UTM 
coordinates of the center line of the PSU using the UTM grid 
origin and the PSU number. These UTM coordinates were then used 
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2 . 2 . 2. 2 


— Continued. 


to do the plotting and labeling of the 200 PSU’s. The PSU*s 
were also located on the 1:250,000 USGS Grand Canyon quad sheet. 
These maps were then sent to the BLM for LSP acquisition. 

2.2.3 Interpretation of Secondary Sample Units (LSP Plots 

Along PSU's) . 

2. 2. 3.1 Initial Interpretation of Secondary Sargple Units 

(SSU's) for Selection o f Plots for Ground Visit . 

Each photo plot of each PSU (15 plots per PSU) was 
sorted into one of the summary categories used in the original 
PSU allocation and selection procedures above (200 PSU's total 
resulting in 3000 photo plots) , Based on the primary vegetation 
cover types noted for each plot an interpreter tallied, by PSU, the 
total number of SSU's for each summary category. The PSU's were 
grouped according to their vegetation strata (forest, woodland 
and range) and then each summary category was totaled across 
PSU's for each strata. The results of the planning model runs 
specified the total number of ground plots to be selected for 
each summary category by strata (see Table 2-22 ) with 

selection by equal probability (random selection with replacement) . 

Since the original SSU's were 35mm color transparency 
stereo pairs, it was considered inappropriate to use them 
during ground data collection efforts. Instead, the transparencies 
for each SSU selected for ground visit were printed as three inch 
by five inch photo pairs and then labeled as to PSU and SSU number 
and left or right stereo position. These prints were supplied to 
the ground crews for field navigation along with maps showing 
PSU location and 1:30,000 scale resource photography showing 
plot location. Also, in the case of the forest and woodland plots, 
the transect lines to be visited and the trees to be measured 
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2. 2. 3.1 — Continued . 

on each transect were annotated on the photos to ensure 
the ground data coincided with the detailed photo inter- 
pretation for estimating tree volume (see Section 2. 2. 3.2). 

The photo prin'-.s were delivered to the field 
crews, and the transparencies from the selected ground plots 
were returned to the photo interpreter for detailed 
interpretation. It is important to point out that the 
interpreter was not informed as to which plots were 
being measured in the field to avoid potentially biasing 
his efforts on those plots with respect to nonvisited 
plots. 


2. 2. 3. 2 Detailed Photo Interpretation to Produce Photo 

Estimates . 

Detailed photo interpretation of each stereo pair 
was accomplished utilizing a Richards Model GFL-3040 light 
table, a portable 4X stereo viewer, and a 7X loop magnifier. 
There were a total of 199 flight lines selected for photo 
interpretation over range, woodland, and forest areas. 

Each flight line contained a total of 15 stereo pairs. 

Thus, a total of 2985 stereo pairs were interpreted. Ad- 
ditionally, each stereo pair had an associated small 
scale (wide angle) transparency. This was utilized to 
determine differing cover types within a 75 meter radius 
of the center of the transparency. 
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2 . 2 . 3. 2 


Coni', inued. 


One person conducted all the photo interpreta- 
tion activity.. This was done to eliminate the multiple 
biases that may occur using multiple interpreters. Thus, 
any bias in the interpretation of species type, percent 
ground cover, height,, or crown diameter iis constant and can 
more easily be corrected. 

Species composition information was recorded 
lor each vegetation type (range, woodland and forest) . 

This data consisted of species identification (see Table 
2-23 for the species codes) and cover to the nearest 1%. 
Additionally, each species was assigned to a height class 
(see Table 2-24) and evaluated for foliar density (see 
Table 2-25) . This data was recorded on the forms illus- 
trated in Figures 2-11 and 2-12) . Flight line and plot 
numbers, scale, photointerpretation template number and 
homogeneity index code (see Table 2-26) were also recorded. 

The homogeneity index was obtained by inter- 
preting the area within a template (see Figure 2-1 3C) 
representing a nominal ground distance of 75 meters in 
radius on the small scale (wide angle) transparency 
associated with each plot. The two digit homogeneity index 
consisted of the first digit obtained fr im interpreting 
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Table 2-23. Species Codes. 


Species Species 


Code 

or Ground Cover 

Code 

or Ground Cover 

01 

Apache Plume 

25 

Serviceberry-CeanothuB 

02 

Arroweed 

26 

Shadscale 

03 

Aspen 

27 

Snakeweed 

04 

Big Rabbi thrush 

28 

Turbine 11a Oak 

05 

Big Sagebrush 

29 

Turpentine Bush-Jimmy Weed 

06 

Blackbrush 

30 

Willow-Salix and Chilopsis 

07 

Bursage 

31 

Winterf at 

08 

Cholla 

32 

Yucca 

09 

Clif frose 

33 

Other Desert Shrub 

10 

Cottonwood 

34 

Other Tree 

11 

Creosote 

35 

Other Cactus 

12 

Fir 

36 

Perennial Grasses 

13 

Four-Wing Saltbush 

37 

Annual Grasses and Forbs 

14 

Gambel's Oak 

38 

Bare Ground-Rocky Soil 

15 

Joshua 

39 

Bare Ground-Sandy Soil 

16 

Juniper 

40 

Water 

17 

Little Rabbitbrush 

41 

Riparian Shrub 

18 

Locust 

42 

Shadow 

19 

Manisanita 

43 

Sage-like Low Desert Shrub 

20 

Mixed Chapparal/Mixed Mtn. Shrub 

44 

Deadbrush 

21 

Mesijuite-Acacia 

45 

Dead Tree 

22 

Pinyon Pine 

46 

Agriculture 

23 

Ponderosa Pine 

98 

Uninterpre table 

24 

Saltcedar 

99 

(Dark Exposure) 

Uninterpretable 

(No Stereo Pair or Overlap) 
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Table 2-24. Height of Vegetation. 


CODE 


HEIGHT (ft.) 


1 

2 

3 

4 

5 

6 


0-1 

1-3 

3-6 

6-15 

15-25 

25+ 


Table 2-25. Toliar Density of Vegetation 


CODE 


DENSITY {%) 


1 

2 

3 

4 


0-25 

26-50 

51-75 

76-100 
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ure 2-12. PHOTO INTERPRETATION FORM FOR WOODLAND AND FOREST AREAS 






Table 2-26. Homogeneity Index Codes. 


Code 


Descriptio n 


1 One cover type within the plot 

2 Two cover types within the plot 

3 Three or more cover types within the plot 

“0 No other cover types within 75 meters of 

the plot center 

-1 Only one other cover type within 75 meters 

of the plot center 

-2 Two other cover types within 75 meters of 

the plot center 

-3 Three or more other cover types within 75 

meters of the plot center 
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Nominal Scale Scale Range 
1:600 1:530-1:649 


Nominal Scale 
1:600 


Scale Range 
1:550-1:649 



A. Example Woodland and Forest B. Example Range I»lot Template. 

Transect Template. 


Nominal Scale Scale Range 
1:600 1:550-1:649 



C. Example Homogeneity Index Templates. 


Figure 2-13. Photo Interpretation Templates. 
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2.2. 3.2 — Continued 


the large scale transparency (used for species composition) and 
a second digit from interpretation of the small scale transparency 
using the 75 meter homogeniety index template. 

The photo interpreter was familiar with each of the various 
species types encountered in the rangeland, woodland, and forest 
areas. This familiarity was obtained while visiting the test area 
with BLM personnel prior to the beginning of all photo interpreta- 
tion. While in the field, ground photos and descriptions of 
each species type were recorded. Upon return from fhe field, 
this data was compiled into a photo interpretation aid and used 
by the interpreter to assist in species identification (see 
Appendix 2-E) . 

2. 2. 3. 2.1 Range. 

A total of 108 rangeland flight lines (15 stereo pairs each) 
were interpreted. Each stereo pair v;as recorded by flight line 
and plot number on the form found in Figure 2-1?. Each pair of 
transparencies were aligned to stereo and a scaled template over- 
layed. The template represented a nominal ground area of 42 feet 
by 48.75 feet and was positioned with its center corresponding to 
the center of the stereo overlap of the transparencies (see Figure 
2-13B). The scale of the stereo pair determined which template was 
used for interpretation (see Table 2-27) . Only the ground area 
falling within the boundary of the template was 
interpreted. All types of ground cover found within a template 
boundary were recorded individually to the nearest 1%. All ground 
cover types contained within a template totaled 100% cover. 

Any portion of the ground area falling within a template, 
but not within the stereo overlap, was interpreted for ground cover 
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Table 2-27 


Photo scales and templates used 
in the photo-interpretation. 


TEMPLATE 1 

TEMPLATE BASED 
ON NOMINAL SCALE OP 

RANGE OF PHOTO 
SCALES APPLIED 

1 

1:500 

450-549 

2 

1:600 

550-649 

3 

1:700 

650-749 

4 

1:800 

750-899 

5 

1:1000 

900-1099 

6 

1:1200 

1100-1349 

7 

1:1500 

1350-1749 

8 

1:2000 

1750-2499 

9 

1:3000 

25004- 
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2 . 2 . 3 . 2*1 


--Continued, 


type. The percentage of the total template area within this 
criteria, for any given stereo pair, was recorded as percent nonstereo 
coverage to the nearest 10 percent. For each type of vegetation 
found within a stereo-pair template, an estimate of height and 
foliar density was recorded. These were recorded on the PI form 
with the codes listed in Tables 2-24 and 2-25 . These measurements 
were intuitive and not based on photogrammetric procedures. 

2. 2. 3. 2. 2 Woodland . 

A total of 45 woodland flight lines (15 stereo pairs 
each) were interpreted. Species composition data for each stereo 
pair was recorded by flight line and plot number on the form 
found in Figure 2-11 . Each photo pair was adjusted to stereo and 
the species and ground cover types contained within the entire area f 
stereo overlap were r^i'oorded individually. Each cover type was 
recorded to the nearest 1% of the total stereo overlap. Heights 
and foliar densities for each species were also recorded using 
the same method as described for range (Sec. 2. 2. 3. 2.1). 

A scaled transect template (Fig.2-13A, p. 2-99) was 
applied to each stereo pair based on the scales listed in Table 2-27. 
The template consisted of two parallel horizontal lines which tri- 
sected the area of a single 35mm slide. The length of the hori- 
zontal transects represented a nominal ground distance of 75 feet. 

All transect templates were overlayed on the left frame of a 
stereo pair ond were registered to the right edge of the left 
frame. The right frame was then adjusted to stereo and the tem- 
plate adjustment was complete (see Figure 2-14 ) . 

Additional information for each stereo pair was also 
recorded on this form. This information included flight line and 
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2 . 2 . 3 . 2. 2 


— Continued 


plot numbers, scale (as supplied by the aerial photographer), tem- 
plate number applied, total number of trees, and percent non- 
coverage. Percent non-coverage refers to the total amount of 
transect line length which does not fall within the overlap of a 
stereo pair. It is recorded to the nearest 10 percent as transect 
line non-coverage. 

Any Pinyon Pine or Juniper tree whose crown intersected 
any part of either transect line was interpreted. Each tree was 
numbered in sequtv^ce starting with the left-most tree on the 
upper transect line and proceeding right across the line. When 
all trees on the upper transect line were numbered, the numbering 
sequence was continued with the left-most tree on the lower 
transect line and proceeded to the right through the last tree to 
intersect the lower line. This numbering sequence is illustrated 
in Figure 2-14. 


Each numbered tree on the transect lines was analyzed 
and specific information recorded or; the tree measurement form 
(see Figure 2-15). Species was recorded by code number (see Table 
2-23) . The long and short dimensions of the tree crown were 
measured to the nearest 1/100 of an inch and recorded. Finally, a 
height estimate was recorded to the nearest 5 foot increment for 
each tree. Tree heights were not photogrammetrically measured and 
were thus recorded as estimates based on relative size of surround- 
ing trees and a general Jcnowledge on the part of the interpreter 
of the range of tree heights to be expected for a given species. 

2.2. 3.2. 3 Forest . 

A total of 46 forest area flight lines (15 stereo pairs 
each) were interpreted. Species composition data for each stereo 
pair was recorded by flight line and plot number on the form shown 
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in Figure 2-11. Similar to the woodland are^««, no area 
template was applied to each photo pair, and therefore, a 
0 was recorded in the Template Number column on this form. 

Each photo pair was adjusted to stereo and the percent of 
overlap was analyzed. The nominal stereo overlap was 60 per- 
cent forward. Any overlap less than this was recorded in 
the % non-coverage column of the form. The number recorded 
in this column represents the percent of stereo overlap 
missing from the nominal 60 percent forward overlap. This 
figure was recorded to the nearest 10 percent. 

All ground cover types and species found within the 
entire area of stereo overlap were identified and recorded 
by code number (see Table 2-23) . Each voer type was recorded 
to the nearest 1% of total cover within the area of stereo 
overlap. All cover types totaled 100 percent cover. Heights 
and foliar densities for each cover type were also recorded 
by code number (see Tables 2-24 and 2-25) using the method 
described for Range (Section 2. 2. 3. 2.1). 

The pvoper transect overlay was then applied to 
each stereo pa;Lr as per the procedures described for Woodland 
(Section 2. 2. 3. 2. 2). The transect applied was based on the 
scale of the stereo pair as supplied by the aerial photographer. 
The various templates are listed in Table 2-27. Figure 2-14 
depicts the method employed to overlay the template on the 
stereo pair. 

Any Ponderosa pine trees that intersected the 
transect lines were assigned a tree number. The number 
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sequence was the same as was used in the analysis of the woodland 
areas (Section 2. 2. 3. 2. 2). Each tree number was recorded on the 
tree measurement form (see Figure 2-15) . This form also contained 
information concerning flight line and plot numbers, scale, the 
template number applied, total number of trees, and the percent of 
template non-stereo coverage. Each numbered Ponderosa pine on a 
transect line was analyzed and specific information was recorded on 
the tree measurement form. The long and short dimensions of each 
tree crown were measured to the nearest 1/100 of an inch. A 
height estimate was also recorded and was based on the size of 
surrounding trees and a general knowledge of the interpreter about 
this species type. This number was not a» photogrammetric measure- 
ment, but was an estimate and was recorded to the nearest 5 foot 
increment. 

A second tree measurement form was used in the forest 
areas to record all Pinyon Pine and Juniper trees which touched 
or intersected the transect lines. Information similar to that 
recorded for the Ponderosa Pine was recorded for these two tree 
types using the methods previously described. A second tree 
measurement form for these two species types was used for ease of 
data handling due to the high density of trees and tree types 
found on some forest area flight lines and plots. This data was 
not used in the estimation procedures described in Section 2.3. 

The information is provided for use by the BLM for estimation of 
Pinyon and Juniper volume in the forest areas if that becomes 
an item of interest at a later date. 

2.2.4 Data Collection on Ground Plots . 

The data collection procedures for the ground plots were 
designed to be as close to BLM procedures as the multistage sample 
design would allow. The procedures used on the range plots were 
modified from SVIM while the forest and woodland procedures were 
based on the BLM Extensive Forest Inventory Field Handbook. 
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2.2.4. 1 Range Plots 


Data collection on the 135 range plots was designed to 

be as close to SVIM procedures as practical. Plot locations were 

defined as the cf?nter of stereo overlap of a selected pair or 

LSP. This location was pin-pricked on LSP prints and also 

located on BLM 1:30^000 resource photography to assist in location 

of the plot. This location was used as the center of a 42 foot 

by 48.75 foot plot. This plot consisted of a 200 point grid on 

5 lines with a spacing of 15 inches along a line and a distance 

of 10.5 feet between lines. These 200 points were used to 

determine the species composition of the plot. A systematic set 

2 

of 20 subplots (.1 M ) were used to take SVIM shrub and weight 
characterization measurements. This data was recorded on 
modified SVIM forms (see Appendix 2-F) . Additional information 
was recorded for each plot as follows: the name of the crew mem- 
ber who made the weight estimates and the ground distance for 
two points on the LSF for use in making photo scale calculations. 

2. 2. 4. 2 Forest and Woodland Plots . 

Plot location for the 16 woodland plots and 27 forest 
plots was determined in advance by ESL drawing two transects 
(75* nominal) in the stereo overlap of each selected stereo pair. 
These transects were annotated on the photography and the 
individual trees to measure were pin-pricked and numbered. 

The information recorded for each tree was diameter (ground and 
stump height for woodland and breast height for forest ) , height, 
crown diameter, number of stems and growth and age information 
on a subset of the trees. The data coding forms and ground 
data collection procedures are detailed in Appendix 2-F. 
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Data Reductions. 


All data reductions tasks were performed using the Earth 
Resources Inventory System (ERIS) . The photo interpretation data 
and ground plot measurements were keyed to ASCII text form from the 
data sheets and then converted to a structured ERIS binary file 
format. The Landsat data was input into ERIS files by means of 
several IPIMS functions which create image summaries. 

2. 3. 1.1 Ground Data to Ground-Based Estimates. 


2. 3. 1.1.1 Range Ground Data Reductions . 

There were three basic sets of range ground data 
for each of the 135 ground plots: 

(1) Transect summary data 

(2) Weight estimation and vegetation characterization 
data 

(3) Dry /green weight conversion data 

2. 3. 1.1. 1.1 Transect Summary Data Reductions . 

This data set was a summarization of the 200 points laid 
out in 5 transects of 40 points each for every plot. The cover type 
was identified for each subplot. Nonvegetation types were separated 
into 9 categories. Vegetation types were identified by species where 
possible. Each vegetation subplot had a basal type determination 
and also had type determinations by canopy, if any, up to 3 canopies. 

The data set, consisting of measurements on 135 plots, was 
keypunched into ASCII text format and then converted into ERIS binary 
format. The vegetation basal and canopy type determinations were in 
alphanumeric form. The types occurring were summarized and checked 
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against the project plant and cover type list. There were some cate- 
gories occurring in the data which did not appear on the list. After 
referring to the original data sheets, most of the reasons for these 
discrepancies were apparent, and corrections were made. There still 
remained a few categories not on the list, though, and these were 
added to the list. The final amended list appears in Table 2-28. 

Another screening procedure performed was to verify that 
there were 200 transect "hits" recorded per plot. It was found that 
25 of the 135 plots had either more than 200 or less than 200 hits 
recorded. Examination of the data sheets resolved 7 of these - either 
the dot counts were not summarized correctly on the sheets or the 
nun\bers were not keyed correctly. Of the remaining 18 plots all but 
1 were in the 198-203 range, the other being 190. VIhen these counts 
were used at a later stage for plot vegetation summaries it was not 
assumed that 200 hits existed - proportions based on the recorded total 
were used. 


In order to proceed further with the analysis, it was nec- 
essary to restructure the data file. The alphanumeric codes for 
vegetation/cover types were changed to numeric codes (see Table 2-28) . 
Also, the numbers of hits recorded referred to multiple canopy com- 
binations which often involved more than one plant type. The hits 
in these cases were distributed equally among the vegetation cate- 
gories which occurred. The data could then be formatted in terms of 
number of hits by plant type within each plot. 

2. 3. 1.1. 1.2 Weight Estimation and Vegetation Characterization 

Reductions . 

This data set contained the forage weight measurements in 
grams for subplots within a plot. There could be as many as 20 

{ 
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measured subplots from each plot (out of 200 total) , although the 
number was usually much smaller because nonvegetation types had no 
measurements recorded. The measurements were ocular estimates which 
were made separately by species and by four height categories. Also 
recorded were availability and utilization factors, phenology, average 
height, average crown diameter, age class and form class. The ocular 
estimates, species, availability and utilization were screened and any 
apparen errors were corrected. 

To produce the subplot weight estimates, the values for 
the four height classes were first added together. The results were 
then input into a prediction equation to yield estimates of air dry 
weight. A separate equation was used for each of four individuals 
who did the ocular measurements. These equations were based on a 
linear regression performed on a sample of subplots (section 2. 3. 1.1. 3) 

The availability and utilization factors were then applied. 
The availability factors were as follows: 

Code Percent Factor 


A 

100 

1 

P 

75 

.75 

H 

50 

.5 

L 

25 

.25 

U 

0 

0 
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This factor was multiplied by the predicted air dry weight. The 


utilization factors were 

as follows: 


Code 

Percent Range 

Factor 

0 

0 

1 

1 

1-20 

.9 

2 

21-40 

.7 

3 

41-60 

.5 

4 

61-80 

.3 

5 

81-100 

.1 


Originally, thir. factor was used as a multiplier to the air dry 
weight, but it was later discovered that it should have been used 
as a divisor instead. This error was corrected at a later stage 
by applying average factors by Landsat strata (see Section 
2. 3. 2. 3.1). Also retained were the predicted weights without the 
availability and utilization factors applied. This was also carried 
through the estimation, representing total current air-dry weight, 
whether available or not. The primary value, with the factors 
applied, represents air-dry weight available and projected for 
full utilization. 

The weight measurements for those species which are 
unpalatable were then deleted. The criterion of what was to be 
included as palatable was a proper use factor in the summer of 
at least 10%. Table 2-28 shows what species were considered as 
palatable (indicated by a 1 in the column PALAT) . 

The results were then averaged over any subplots which 
had a given species occurring, yielding average weights by species 
within each plot, as shown in equation 2. 3. 1.1; 
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Table 2-28. Range Ground Data Species Summary. 


PLAMTNUM 

PLANTSYM 

SC I NAME 

TRANHITS 

PALAT = 

WTSAMPS 

AVPRWTAU 


= Numerical species code 
* Species symbol 
= Scientific name 
= Number of transect hits 

1 if palatable species^ = 0 if not palatable or non— veg 
® Number of subplots with ocular v;eight estimates 
= Average weight after ocular to dry regression 

adjustment and availability and utilization factors 
applied (-1 indicates no samples). 


PLAimfDN 

PUUrrBYH 

8GIRAJIE 

TK<.NH1TS 

PAUT 

WTSAMPS 

AVPRWTAU 

1 

ACCR 

ACACIA imEGCIl 

3 

9 

1 

17.2968 

a 

AFD 

ANNUAL FORB DIIAO 

781 

9 

9 

-1 .99099 

a 

AFC 

ANNUAL FORB GREEN 

499 

9 

9 

-1 .90099 

4 

ACAVl 

AGAVE (SPECIES) 

19 

9 

9 

-1 .00009 

S 

AGO 

ANNUAL CRASS DEAD 

4128 

9 

9 

-1 .90009 

4 

ACC 

ANNUAL CRASS GREEK 

8 

9 

9 

-1 .00009 

7 

AHDU 

AMBROSIA DUHOSA 

119 

9 

28 

3.98661 

8 

AMUT 

AMELANCHIER UTABENSlS 

39 

9 

12 

4.99949 

9 

ARI8 

ARISTIDA (SPECIES) 

18 

9 

19 

.701964 

19 

AREN 

ARENARIA (SPECIES) 

9 

9 

9 

-1 .00099 

il 

ARFE 

ARISTIDA FENDLERIANA 

4 

9 

2 

1.89211 

12 

ARLO 

ARISTIDA LONGISETA 

29 

9 

16 

2.08747 

13 

ARLU 

ARTEMISIA LUDOVICIANA 

34 

9 

9 

-1 .90099 

14 

ARPA 

ARGHTOSTAPHYLOS PATULA 

199 

9 

9 

-1 .90009 

18 

ARPR 

ARCHTOSTAPHYLOS PRINGLtI 

28 

9 

9 

-1 .00009 

14 

ARPU 

ARISTIDA PURPUREUN 

81 

9 

1 

6.44363 

17 

ARPU 

ARCHTOSTAPHYLOS PUNCENS 

9 

9 

9 

-1 .99009 

18 

ARTE 

ARISTIDA TERNIPES 

4 

9 

9 

-1 .00099 

19 

ARTR 

ARTEMISIA TRIDENTATA 

455 

9 

9 

-1 .00009 

29 

ATCA 

ATRIPLEX CANESCENS 

24 

1 

19 

7.99029 

21 

ATCO 

ATRIPLEX CONFERTI FOLIA 

7 

9 

6 

19.8262 

22 

BG 

BARE GROUND 

2869 

9 

9 

-1 .90009 

23 

BEDROCK 

BED ROCK 

938 

9 

9 

-1 .99009 

24 

BOER 

BOUTELOUA ERIOPODA 

81 

1 

24 

1.79089 

28 

BOCR 

BOUTELOUA GRACILIS 

42 

1 

32 

1.28279 

26 

CECR 

CEANOTHUS GHECGII 

34 

1 

7 

8.49473 

27 

CHLl 

CHILOPSI8 LINEARIS 

38 

9 

9 

-1 .00099 

28 

CHNO 

CHRYSOTHAMNUS M0HAVENSI8 

8 

9 

9 

-1.00099 

29 

CHNA 

CHRYSOTHAMNUS NAUSE06US 

3 

9 

9 

-1 .00009 

39 

CHOLLA 

CHOLLA CACTUS 

4 

9 

9 

-1 .00009 

31 

CNPA 

CHRYSOTHAMNUS PANICULATUS 

88 

9 

9 

-1.00009 

32 

CIIVA 

CHRYSOTHAMNUS VASEYl 

6 


9 

-1 .00099 

33 

CHVI 

CHRYSOTHAMNUS VISCIDIFLCHUJS 

46 

9 

9 

-1 .90009 

34 

CIRS 

CIRSIUH (SPECIES) 

1 

9 

9 

-1.00009 

38 

COBBLE 

COBBLE 

1484 

9 

9 

-1.00009 

36 

CONE 

COWANIA HEXICANA(SAHE AS COST) 

9 

1 

9 

-1 .00009 

37 

COPA 

COMANDRA PALLIDA 

2 

9 

9 

-1 .00009 

38 

CORA 

COLEOCYNE RAMOSISSIMA 

389 

9 

194 

9.68421 

39 

COST 

COWANIA STANSBURIANA 

137 

1 

32 

13.6943 

49 

DAFR 

DALEA FREMONTII 

8 

9 

9 

-1 .00099 

41 

ELYH 

ELYMUS (SPECIES) 

1 

9 

9 

-1 .00009 

42 

ENCE 

ENCELIA (SPECIES) 

1 

9 

9 

-1.90009 

43 

ENFR 

ENCELIA FRUTESCENS 

16 

9 

9 

-1 .90099 

44 

EPTO 

EPHEDRA TORREYANA 

38 

9 

8 

1 .69624 

48 

EPVI 

EPHEDRA VIRIDIS 

179 

9 

62 

9.41823 

46 

ERIO 

ERIOGONUN (SPECIES) 

49 

9 

9 

-1 .00009 

1 

ERGO 

ERIOCONUK CORYNBOSUH 

3 

9 

9 

-1 .00099 

48 

ERCO 2 

ERICERON CONCINNUS 

9 

9 

9 

-1 .00009 

49 

ERFA 

ERIOGONUN FASCICULATUH 

86 

9 

4 

7.73832 

89 

ERIN 

ERIOGONUN INFLATUN 

46 

9 

9 

- 1 .eoooo 

81 

ERLA 

ERIOPHYLLUM LANOSUM 

1 

9 

9 

-1.99099 
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puumnni 

PUItTWH 

8CINAHE 

nuimiTB 

PALAT 

VTSAMPS 

AVPRirrAU 

B2 

Ejum 

ERIOGONUN imBELLATUN 

8 

4 

18 

2.47942 

83 

CfWR 

ERIOGONUN WRIGHTII 

18 

4 

4 

8.84498 

84 

Et'AL 

EUPHORRIA ALRONARGINATOS 

8 

4 

4 

-1.44444 

88 

F.ULA 

EUROTIA LANATA 

14 

1 

8 

4.14234 

84 

FAPA 

FALLUGIA PARADOXA 

27 

4 

1 

1.84814 

87 

FERO 

FFJIROI BARREL) CACTUS (SPECIES) 

2 

4 

4 

-1 .40044 

88 

FRVE 

FRAXINUS VELirriNUS 

4 

4 

4 

-1 .44044 

84 

CAFL 

GARRYA FUVESCENS 

31 

4 

4 

-1.44044 

«• 

CLNC 

GLOSSOPCTALON NEVLJENSI 

19 

4 

4 

8,48633 

«l 

GRAVEL 

GRAVEL 

4244 

4 

4 

-1 .44044 

43 

CUBA 

GUTIERREZIA SAROTHRAE 

322 

4 

4 

-1 .44444 

43 

HIJA 

HILARIA JANESII 

317 

1 

181 

2.24217 

44 

HIRI 

HILARIA RIGIDA 

4 

1 

1 

.278808 

48 

HY8A 

KYNENOCLEA SALSOLA 

18 

4 

4 

-1 .40004 

44 

juno 

JUNIPERUS NONOSPERNA 

1443 

4 

4 

-1 .44004 

47 

JU08 

JUNIPERUS 06TE0PERHA 

141 

4 

4 

-1.44444 

48 

KOGR 

KOELERIA CRISTATA 

4 

1 

8 

1.66887 

49 

KRPA 

KRAHERIA PARVI FOLIA 

84 

4 

13 

8.46136 

74 

LI CHER 

LICHEN 

2 

4 

4 

-1 .44044 

71 

LATR 

LARREA TRIDENTATA 

228 

4 

4 

1 . 14446 

72 

LEFR 

LEPIDIUH FREHONTIl 

2 

4 

4 

-1 .00004 

73 

LIVR 

LIPPIA WRIGHTII 

14 

4 

4 

-1 .44044 

74 

LOWR 

LOTUS WRIGHTII 

4 

1 

1 

1.23984 

78 

LYAN 

LYCIUH ANDERSONII 

34 

4 

2 

1.46181 

74 

MOSS 

MOSS 

222 

4 

4 

-1 .40004 

77 

HAFR 

NAHONIA FREHONTIl 

11 

4 

4 

-1 .40004 

78 

HESQ 

PROSOPIS JULIFLORA — HESQUITE 

8 

4 

4 

-1 .44044 

79 

NUPO 

nUHLENBERGIA PORTERIA 

18 

1 

4 

3.41784 

84 

NA 

NOTHING TO CLIP 

4 

4 

4 

-1 .44044 

8t 

NG 

NOTHING TO CLIP 

4 

4 

4 

-1 .40044 

82 

HP 

NOTHING TO CLIP 

4 

4 

4 

-1 .44044 

83 

HPL 

NON-PERSISTENT LITTER 

4322 

4 

4 

-1 .44044 

84 

OPUHT 

OPUNTIA (SPECIES) PRIC3CLY PEAR 

34 

4 

4 

-1 .40004 

88 

OPAC 

OPUNTIA ACANTHOCARPA 

1 

4 

4 

-1.00044 

84 

ORilY 

ORYZOPSIS HYHENOIDES 

19 

1 

3 

1 .86836 

87 

PENS 

PENSTEHON (SPECIES) 

4 

4 

4 

-1 .44004 

88 

PlEO 

riNUS EDULIS 

2 

4 

4 

-1 .00044 

89 

PINO 

PINUS HONOPHYLLUN 

437 

4 

4 

-1.44044 

94 

PlPO 

PINUS PONDEROSA 

84 

4 

4 

-1 .44044 

91 

PL 

PERSISTENT LITTER 

942 

4 

4 

-1 .40044 

92 

POA 

POA (SPECIES) 

4 

1 

4 

.879843 

93 

POFE 

POA FENDLERIANA 

18 

1 

9 

1 .21447 

94 

POLO 

POA LONGILIGULA 

92 

1 

49 

1.13211 

98 

PONO 

POTENT I LLA H0DE8TA 

2 

4 

4 

-1 .40044 

94 

PRFA 

PRUNUS FASCICULATA 

9 

4 

3 

2.89864 

97 

P8AR 

PSOROTILINNUS AROORECENS 

19 

4 

3 

9.12767 

98 

OUGA 

OUFRCUS GAHBELII 

8 

1 

4 

-1 .00444 

99 

OUTU 

OUERCUS TURBINELLA 

821 

4 

78 

8.91878 

144 

RHTR 

RHUS TRILOBATA 

38 

4 

3 

2.97877 

141 

SAGA 

SALAVIA CARNOSA 

7 

4 

4 

-1.40004 

142 

8GBR 

SCLEROPOGON BREVEF3LIUS 

7 

4 

8 

1.49447 
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PLAimfim 

PLAirrSYK 

SGIRAHE 

TRANBIT8 

PALAT 

VrSAHPS 

AVFRWTAU 

l«3 

simr 

SITANIOS HYSTRIX 

40 

1 

48 

.837029 

164 

8P0R 

SPOROBOLUS (8PEGIES) 

2 

1 

0 

-1 .00000 

les 

SPAN 

8PHAERALGEA ANBIGUA 

13 

1 

18 

3.21408 

106 

8PC0 

8PHAFJULGEA GOGGIMEA 

3 

0 

1 

.000000 

107 

8PER 

SPOROBOLUS CRYPSTAMDRUS 

0 

1 

0 

-1 .00000 

108 

SPGR 

SPHAERALGEA GROSSULARIAEPOLIA 

16 

0 

8 

2.37006 

109 

8P0R 

SPOROBOLUS (8PEGIES UlflOfOVIl) 

0 

1 

3 

1 .47252 

110 

SPNE 

SPOROBOLUS NEALLEYl 

0 

0 

1 

.084024 

111 

STAn 

STIPA ARIDA 

6 

0 

2 

1.68084 

112 

STGO 

ST I PA GOROKATA 

12 

0 

6 

.721519 

113 

8T8P 

STIPA SPEGI06A 

147 

0 

48 

.768248 

114 

TETR 

TETRADYNIA (SPEGIES) 

9 

0 

0 

-1 .09000 

IIS 

TEAX 

TETRADYHIA AXILLARIS 

12 

0 

0 

-1.00000 

116 

THNO 

THANHOSHA HORTANA 

12 

0 

0 

-1 .00000 

117 

THNE 

URKROVR 

2 

0 

0 

-1.00000 

118 

UP 

UlfKROWR PORB 

8 

0 

0 

-1 .00000 

119 

UG 

UlfKNOVR GRASS 

4 

0 

2 

1 .37232 

120 

US 

UmCROVN SHRUB 

115 

0 

1 

2.49675 

121 

YUCA 

YUCCA (SPECIES) 

3 

0 

0 

-1 .00000 

122 

ARSP 

ARISTIDA (SANE AS ARFS) 

0 

0 

0 

-1.00000 

123 

ARIffTlDlA 

ARISTIDA (SANE AS ARIS) 

0 

0 

0 

-1 .00000 

124 

STSC 

7 

3 

0 

3 

.000000 

12S 

BOGU 

BOUTELOUA CURTIPElfDULA 

8 

1 

1 

2.61046 

126 

STONE 

STONE 

969 

0 

0 

-1.00000 

127 

BANU 

BAILEYA NULTIVAOIATA 

1 

0 

0 

-1 .00000 

128 

ERPO 

ERIOGONUN POLIPOLIUN 

2 

0 

0 

-1.00000 

129 

UG 

UNKNOWN CONPOSITE 

8 

0 

0 

-1 .00000 

130 

YUAN 

7 

36 

0 

0 

-1 .00000 

131 

YUBA 

7 

68 

0 

0 

-1 .00000 

132 

YUBR 

7 

6 

0 

0 

-1 .00000 

133 

YUDU 

7 

2 

0 

0 

-1.00000 
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2 . 3 . 1 . 1 . 1.2 


Continued. 



'^ijkl ^ijkl (2. 3. 1.1) 


where: 


^ijkl “ 

^ijkl “ 
^ijkl “ 
°ijk " 

4jk “ 


predicted air-dry weight in grams for palatable 
species k on subplots 1 from plot j of PSU i 


availability factor 


utilization factor 


set of subplots of palatable species k in plot j 
of PSU i which have weight samples 


cardinality of 


2. 3. 1.1. 1.3 iJry /Green V7eiqht Conversion Data Reduction . 

A subset of the subplots which had ocular weight estimates 
were also weighed after clipping and again after air drying. The 
relationship between the ocular estimates and the air-dry weights 
required estimation to predict the air-dry weight for the larger 
sample of ocular subplots. There were 401 subplots in total which 
had both air-dry weights and ocular estimates, the latter being 
performed by four different individuals. Actually, the coding showed 
three other people did a small number of subplots (about 20) . These 
were combined with the individual whose data structure was most 
similar. A scatter plot of these two variable for all measurements is 
shown in Figure 2-16. The correlation coefficient is 88%, which 
yields a multiple correlation coefficient (r 2) of 77% for a linear 
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2 . 3 . 1 . 1. 3 


Continued. 


regression. Table 2-29 sununarizes the results of separate 
regressions for the different ocular estimators. The R^-values 
are seen to be 55%, 73%, 80% and 95%. The linear regression 
coefficients in Table 2-29 were used for prediction on the full 
set of ocular subplot estimates. This procedure proved to be a 
very effective use of resources for measuring forage weight. 


2. 3. 1.1. 4 Combination of Inputs to Obtain Plot Estimates. 


It was then possible to combine the no'ber of transect 
hits by species/type with the average weights and then sum over 
species/types to get single estimates for each plot. This was 
done as follows: 




K ^ 

^ijk ^ijk 


(2. 3.1.2) 


where 


^i jk 


i jk 


as defined in equation 2. 3. 1.1, 
if weight samples exist for 
that species in that plot 

the overall species k average, 
if weight samples do not exist 
for that species in that plot 


if species is unpalatable 


and 


K 


total number of species/types 


2-117 



Table 2~29 . Dry/Green Weight Regression Summary. 
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2 . 3 . 1 . 1 . 1 .4 — Continued 



number of hits for apeclee k In plob. j of 
PSU 1 


C 


ij 


C 


ij 


total number of transect hits from all 
species/types for plot j of PSU 1 



C 


ijk 


Note that the result is in grams per subplot. Since the subplot size 
was one- tenth square meter, the following cor, versions were made: 


kg/hectare ■ 

100,000 

lO'OD” 

A 

^ij 

Ib/acre ■ 

40,468.6 

A 

^ij 


This yielded the final ground-bases plot estimates of air-dry forage 
weight. 


2. 3. 1.1. 2 Woodland Ground Data Reduction . 

The following data were collected for each stem on the 16 


woodland ground plots: 

Tree number 
Stem number 

Species 

DGH 

DSH 

Height 

I 


Identification number within plot 

Identification number for multiple 
stemmed trees 

Pinyon “ 1 ,Juniper»2 

Diameter at ground height in 
inches 

Diameter at stump height in Inches 
Total height in feet 
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2. 3. 1.1.2 — Continued 


Crown Major 
Crown Minor 
Crown Average 
Age 

Growth 


Major crown diameter to nearest foot 

Minor crown diameter to nearest foot 

Average crown diameter in feet 

Total age from boring; only measured 
for a rample of trees 

Last 10 years growth in inches; only 
measured for a sample of trees. 


This dat<^ set was edited for coding » completeness and consistency. 
A small number of edits were needed for the mu <«urements. A major 
problem was that the ground plot tree numbers did not match the 
tree numbers coded for the corresponding photo interpretation. 
Through additional photo interpretation, the correspondence 
between the trees, photo to ground, was determined for each plot 
and the ground tree numbers were recoded to match the photos. 

Next the records were separated into pinyon and juniper. 
The stem volumes were then computed as a function of DSH, height 
and average crown diameter. For pinyon, the formula used was 
the following: 

Cubic foot volume - Antilog^g Uilog^o (D^H) 

+ ajlo9^„C + ajl 

where a^^ * 1,06779 
a2 « 0.17084 
a^ » -3.10179 
D B dsh in inches 
H « height in feet 
C = average crown diameter in feet 


1 
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2 . 3 . 1 . 1 . 2 — Continued. 


This is the source equation for the pinyon volume table provided by 
the BLM, based on 392 trees in Colorado and New Mexico in 1974 (see 
Table 2-30) . For juniper, the form is the same, but the coefficients 
are: 

a^ - 1.11462 
-*0.02496 
a^ "-3. 08600 

This is the source equation for the juniper volume table provided 
by the BLM, based on 206 trees from the same inventory (see Table 
2-31) . 

The volume for each stem on multiple-stemmed trees were 
added together, which reduced the data from 109 stems to 62 
trees. These were located on a total of 16 ground plots occurring 
in 8 PSU's. 

These tree volumes were then matched agains the photo 
interpretation measurements for the same trees (see section 2. 3. 1.2. 2. 
This was done to get photo-based estimates of tree volumes which 
could then be expanded into plot estimates for each of the 675 photo 
plots. The ground measured tree volumes, however, were not expanded 
into plot estimates directly. 

2. 3. 1.1. 3 Forest Ground Data Reduction. 


2. 3.1.1. 3.1 Tree volume estimates . 

The following data were collected for each stem on the 27 
forest ground plots: 


2-121 


Table 



























original page is 
OF POOR QUALITY 



2-123 





















Identification number within plot 

Identification number for multiple- 
stemmed trees 

Diameter at breast height (4.5*) in 
inches 

Total height in feet 

Total age from boring; only measured 
for a sample of trees 

Last 10 years growth in inches; only 
measured for a sample of trees 

As in the woodland ground data, little editing was required for the 
recorded measurements, but the tree numbers had to be recoded to match 
the photo tree numbers. 

The tree volumes were then computed as a function of DBH and 
height. The formula used was the following: 

2 

Board foot volume ■ aj^ + a2DH, 

if D^H < 20,000 
bi + b2 D^H, 

if D^H > 20,000 

where D » DBH in inches 
H ■ height in feet 
aj^» -1.97333 
a2> .000889682 

bj^» -11,3012 
b2» .00138797 


2. 3. 1.1. 3.1 — Continued . 

Tree number 
Stem number 

DBH 

Height 

Age 

Growth 
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2. 3. 1.1. 3.1 — Continued 


This Is the source equation used In the U. S. Forst Service blackjack 
ponderosa pine table provided shown In Table 2*'32. There were some 
negative estimates (minimum of -1.6) which were changed to zeros. 

Only two trees had more than one stem (two each) . Their 
volumes were added together, resulting In tree volumes for 111 trees 
In total. These were then matched against the corresponding photo 
Interpretation measurements and regression analysis was performed 
(see section 2. 3. 1.2. 3). 

2. 3. 1.1. 3. 2 Plot volume estimates. 


Since the method of sampling trees for ground measurement was 
based on line transects drawn on phots, the ground plot volume esti- 
mates depended on photo measurements. The estimator for board foot 
volume per acre for ground plots Is; 


.. _ 43,560 ^li 

*ij i 1. . E 

A A 

'^Iju / ^Iju (2. 3. 1.3) 

A 

where: 

A 

estimated board foot volume for tree 
u In plot j of PSU 1 as described 
In section 2. 3. 1.1. 3.1. 

W. . 

1]U =» 

estimated crown diameter In feet 
perpendicular to transect 

1 . . “ 
ID 

length of each of the two transects 
on the photo converted to feed on 
the ground 


Number of ground sampled trees on 
plot 

43,560 * 

square feet per acre 
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2. 3. 1.1. 3. 2 — Continued 


The quantities depend on photo measurements And scale 

adjustments and are described in Section 2. 3. 1.2. 2. Using equation 
2. 3. 1.3, resulted in estimates for 27 ground plots occurring in 13 PSU's. 

2. 3. 1.2 Photo Data To Photo-Based Estimates (Including Combinations 
With Ground Data) . 

2. 3. 1.2.1 Range Photo Data . 

2. 3. 1.2. 1.1 Data Reductions . 

This data set consisted of the following for each of 
1620 plots $ 


Scale 

Template 

Percent cover 
Foliar density 
Height class 
Homogeneity index 
The codes and cover types are 


The inverse scale, e.g., scale « 777 
is 1:777. 

Identifier for which template was used 
to overlay the line transects; chis was 
dictated by the scale. 

For each of 48 categories. 

Several categories. 

Several categories. 

Several categories. 

defined in Section 2. 2. 3. 2. 
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2 . 3 . 1 . 2 . 1 . 1 — Continued 


In order to make comparisons between the photo-interpreted 
covcr/density/height and the ground measured forage weights » the 48 
species/cover categories were grouped into these 8 categories s 

Trees 

Desert shrub 
Mountain shrub 
Riparian woodland 
Grasses 
Cactus 

Nonvegatatlon 

Shadow. 

Table 2-33 describes this grouping. The quantities which were added 
together in these groupings were : 

1. Percent crown cover. 

2. Percent crown cover multiplied by foliar density and 
by height (midpoints were used for density and height 
categories) . 

3. Percent crown cover multipled by foliar density. 

Scatter plots and correlations were produced for comparing each of 
these quantities for each group against the forage estimates for the 
135 ground plots. Figures 2-17 and 2-18 are examples of the scatter 
plots for case 1. above. Table 2-34 summarizes the correlations for 
mountain shrub, grasses, desert shrub, and some multiple correlations 
Seve il conclusions can be made from the table. While mountain shrub 
correlated best with the ground measurements, it vras still low. The 
multiple correlations improve the correlations only a moderate amount 
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Table 2-33. Photo-Interpretation Category Groupings 


Group 

Trees 

Desert Shrub 


Mountain Shrub 

Riparian Woodland 

Grasses 

Cactus 

Nonvegetation 


Categories 

Ponderosa pine# Pinyon pine# Juniper# Other 
tree# Aspen# Locust# Dead trees# Fir# Jushua. 

Creosote# Bursage# Blackbrush# Big sagebrush# 
Big rabbitbrush# Cliffrose# Little rabbit- 
brush# Mesquite-acacia# Shadscale# Snakeweed# 
Turpentine bush# Winter fats# Dead brush# 
Arroweed, Apache plume# Four-wing saltbush# 
Other desert shrub# Other low desert sage. 

Gambel's oak# Turbinella oak# Manzanita# 
Serviceberry - ceanothus# Mixed chapparal/ 
mountain shrub. 

Cottonwood# Willow# Riparian shrub. 

Perennial grasses# Annual grasses and forbs. 
Yucca# Cholla# Other cactus. 

Bare ground rocky# Bare ground sandy# Water. 


Shadow 


Shadow. 
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Figure 2-17 
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Figure 2-18 


Table 2-34. Range Photo/Ground Correlations 



Crown 

Cover 

Cover/Density 

Cover/Density/Height 

Photo Category (s) 

A 

B 

A 

B 

A 

. B 

Mountain Shrub 

.185 

.255 

.176 

.248 

. 209, 

.284 

Grasses 

.084 

.051 

.044 

.017 

.006 

-.009 

Desert Shrub 

-.115 

-.153 

-.098 

-.135 

-.063 

-.097 

All 3* 

.226 

.280 

.194 

.259 

.211 

.284 

All 3 and Trees 
and Cactus* 

.328 

.364 

.245 

.302 

.248 

.315 


A s ground estimated forage (amount available and projected for full 
utilization) 

B s ground estimated forage (current amount whether available or not) 
*Multiple correlation coefficient. 

This table is based on 135 plots. 
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2. 3. I. 2. 1.1 — Continued 


The primary estimation parameter is A, with the availability and 
utilization factor applied » but those ground estimates in general 
correlate lower than B. This is because availability and utilization 
were not taken into account for the photo-interpretation (nor could 
they with much reliability) . Another more Important reason for the 
overall low correlations is that unpalatable species could not be 
excluded from the photo interpretation. It is also seen from Table 

2-34 that the height data increases the information more than 
the foliar dersity data alone. 

2. 3. 1.2. 1.2 Bin-Sort Stratification . 

Another photo-interpreted attribute was the bin-sort 
stratification which was used in the sampling of plots for ground 
visits (see Section 2. 2. 1.2.1). Analysis of variance was used to 
see how well this stratification performed in separating differ^snt 
I levels of ground estimated forage. These are shown in Figure 2-19 

for the primary parameter (with availability and utilization factors 
applied) and in Figure 2-20 for the secondary parameter (without 
either factor applied). These show a moderate amount of separability. 
The factor sum of squares divided by the total sum of squarec^ 
represents the proportion of variation explained by the photo strata. 
The square root of this proportion is the analog of the correlation 
or multiple correlation used in linear regression. In Figure 2-19 
this is .385 and in Figure 2-20 it is .395. Since these correlations 
are slightly higher than the beat results achieved from the photo/ 
ground correlations (Table 2-34), the photo stratification is some- 
what more effective than using the regressions corresponding to the 
correlations in that table alone for forage estimation. This by 
itself does not preclude using the regressions as a supplement to 
the stratification. However, analysis of variance between the 
photo stratification and photo crown/density /height measurements 
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ERROR 
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99866.8 
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N 

MEAN ST. DEV. 

1 

18 
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2 

16 
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2. 3.1. 2. 1.2 -> Continued 


yielded analog correlations of well over .8. This Indicates that the 
measurements would add comparatively little new Information. Also, 
predictions under multiple regression models of that type are known 
to be unstable In many circumstances. These considerations led to 
the conclusion that It would be Inadvisable to use the regression 
predictions based on cover, density and heights for the productivity 
estimates. The photo bin-sort stratification was used, however. 

2. 3. 1.2. 2 Woodland Photo Data. 


2. 3. 1.2. 2.1 Data Reductions . 

The woodland photo interpretation data, taken from a 
data set of 675 plots, divided into two sections. The first is the 
same information as recorded for rangeland: species composition, 
foliar density and height. The Second is the result of interpreting 
characteristics of individual trees whose crowns intersect line 
transects overlayed on the slides (see Section 2. 2. 3. 2. 2). The 
atrributes recorded were: 

Species (Juniper = 16, Pinyon pine = 22) 

Crown length (longest dimension) on slide to nearest 

1/100 of an inch 

Crown width (perpendicular to longest dimension) on slide 

to nearest 1/100 of an inch 

Estimated height of tree to nearest 5 feet. 

The scale and template number were also recorded for each plot. This 
latter set of data was used for comparison with the ground data. 
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2, 3. 1,2. 2. 2 Adjustmenta for Scale of Photography . 


Before the photo- Interpreted crown measurements could be 
compared with the ground measurements, they had to be adjusted for 
the scale of the photography. The recorded scales for each plot 
were estimates based on Interpolated altimeter (true altitude) readings 
and ground elevations. These were adjusted using a special paired 
set of ground and photo measurements for 158 plots. Distances between 
specified ground features were measured in inches on either the slides 
or prints (the print measurements were divided by 3.5 to adjust for 
the print/slide scale difference) . The corresponding distances were 
measured on the ground in feet, along with a slope percentage measure- 
ment. The slope percentage was used to convert the ground line-of- 
sight distance to horizontal ground distance as seen from an aerial 
platform: 


g = 

+ ( 100 )^ 

where g * ground distance in feet along line of sight 

a * slope expressed as a percentage (e.g., 54® angle 
from horizontal is 100%) 

g' = horizontal component of of ground distance in feet. 
The scale is then determined by the following: 

S = 

p 

where p = distance in inches on slide. 
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2. 3. 1.2. 2. 2 — Continued. 

f 

The computed scale la then portrayed as 1 t s. These scales are 
plotted against the estimated scales in Figure 2-21, with a .916 
correlatic i. The associated linear regression is shown in Figure 
2-22 » Those coefficients were then used to predict the actual scale 
from the scale as estimated from the aircraft: 

d - (1.146) s' - 43.2 

where s' » aircraft estimated scale. 

This was performed for all the woodland photo plots and all the 
forest photo plots as well. 

These predicted scales were then u?ed to find the approximate 
crovm lengths and widths for the plnyon and juniper trees: 

' °1 “ °1 S / 12 ( 100 ) 

(2. 3. 1.4) 

c^ « c^ S/12 (100) 

where c£ = crown length in 1/100 inches on photo 
c^ = crown width in 1/100 inches on photo 
* crown length in feet on ground 
c^ * crown width in feet on ground 
Also obtained was the geometric mean of the crown length and width: 



2-138 


ORIGINAL page W 
OF POOR QUALITY 


• . • 8 
. - •*«••• • 
«MN S«« N 


» • 


MN ««« N 

«• ' 5 »"< 

*» « 






!I8 


• • **^»M*I» 


s 

i 


I 


i 


s 

i 


tisststssf 

I I i I I i I I i i 


! 

i 


i 

.1 


s 

g 


<N 

I 

<N 

U 

M 

0 

O' 

•H 

Pm 


2-139 



LIHEAR RECRESSKMI FOR PREDICTIIIC EXACT SCAIXTROR EffTIHAlXP 
iiM budatioh is 


I 


ORIGINAL PAGE IS 
OF POOR QUALITY 




X 

0^ 






iisisisjiisiissiisssis 






• •€9*“ • *1# • 

mtm 


* • 




|!3SS8<SSSSS98 


« 

CM 

CM 

I 

CM 

0) 

3 

O' 

•H 

Cm 


I 


2-140 


2. 3. 1.2. 2. 2 — Continued . 

! 

This was used as an estimate of the crown diameter perpendicular to 
the transect, or quantity which would be used in the photo plot 
estimates. 

2. 3. 1.2. 2. 3 Photo/Ground Tree Volume Regressions . 

The photo measurements were matched with the ground measured 

tree volumes (see Section 2. 3. 1.1. 2) and comparisons made separately 

for juniper and pinyon pine. The regressslon and correlation results 

for various photo independent variables are given in Table 2-35. For 

2 

juniper, cases 6 and 12 have the highest R values. However, case 6 

has the double advantage of the lower standard error and having fewer 

terms, so that was the regression selected for predicting volume from 

the juniper photo measurements. For pinyon pine, case 12 had the 
. 2 

highest R . However, it was felt that the crown diameter quadratic 
term (C^) would add instability in a prediction equation. Therefore, 
the equation using only height, crown length and their product, case 6 
was used for pinyon pine as well as juniper. The juniper equation was: 

0 * .260 -.067 h + .022 he -3.29 (2. 3. 1.6) 

The pinyon equation was: 

0 = .567 Cj^ -.098 h -.OOl** hC -3.22 (2. 3. 1.7) 

These were applied to produce photo-based volume estimates for 2004 
juniper trees and 507 pinyon trees. The small number of negative 
estimates obtained were changed to zero. 
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Table 2-35. Woodland Photo/Ground Tree Volume Regressions 


Dependent variable: Cubic foot volume estimated from ground 

measurements 


Pi nyon Pine Juniper 


Photo variables 


si 

R 

S .£. 

Ri 

R 

S . E . 

1 . 



.332 

.576 

2.77 

.198 

.445 

3.67 

2 . 

=1 


.360 

.600 

2.71 

.316 

.562 

3.39 

3 . 

h, Cp 


.337 

.581 

2.85 

.246 

.496 

3.60 

4 . 

h, 


.373 

.611 

2.78 

.336 

.580 

3.38 

5 . 

h, Cp, hCp 


.357 

.597 

2.92 

.247 

.497 

3.64 

6 . 

h , ® 


.373 

.611 

2.88 

.339 

.582 

3.41 

7 . 

^P' V 


.340 

.583 

2.85 

.198 

.445 

3.71 

8 . 

Cl, 


.410 

.640 

2.69 

.317 

.563 

3.43 

9 . 



.342 

.585 

2.95 

.246 

.496 

3.64 

10 . 

°P' °P^' 

hCp 

.390 

.624 

2.96 

.248 

.498 

3.69 

U . 

h, 


.411 

.641 

2.79 

.337 

.581 

3.42 

12 . 

h, ^ 2 .^ ^1 ' 

hCj^ 

.446 

.668 

2.82 

.339 

.582 

3.46 


2 

R = proportion of variation explained by regression 
R ^ multiple correlation coefficient 

S.E. » standard error of the mean estimate 
h - tree height in feet 

= crown length in feet 

c » estimated crown diameter perpendicular to transect. 
P 


2-142 


2. 3. 1.2. 2. 4 Photo Plot Volume Eitlmafa 


The pinyon and juniper volume eatlmates for the photo plots 
reflected the llpe transect sampling method employed. The estimator 
for cubic foot volume per acre was: 



(2.3. 1.8) 


where 



estimated cubic foot volume for tree m in plot 
j of PSU 1 as in equations 2. 3.1.6 or 2. 3. 1.7. 


- estimated crown diameter in feet perpendicular to 
transect (from equation 2. 3. 1.5). 


1 


ij 


length of each of the two transects on the photo 
converted to feet on the ground. 


> number of ground sampled trees on plot. 
43,560 > square feet per acre. 


These estimates were computed for both pinyon and juniper for all 675 
photo plots. 


2 . 3 . 1 . 2 . 3 Forest Photo Data. 


2. 3. 1.2. 3.1 Data Reductions . 

The forest photo interpretation data, taken from a 
data set of 690 photo plots, is similar to the woodland data. 
Again, percent cover, foliar density and height were recorded 
by species and cover types. The tree interpretation was made 
for ponderosa pine, but the same measurements and estimates were 
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2 . 3 . 1 . 2 .3.1 -- Continued 


recorded and the same data sheets used. The measurements were, crovm 
length and width on the slides to the nearest 1/100 inch. The photo 
estimates were for total tree height to the nearest 5 feet. These 
were recorded for all trees whose cro%ms intersected the line tran- 
sects for 690 photo plots from 46 PSU's (although many photo plots 
had no ponderosa pine trees) . 

The crown measurements on the slide were converted to feet 
on the ground by applying the adjusted scale of the photography. This 
was done using exactly the same procedure as for the woodland photo 
measurements and is detailed in Section 2. 3. 1.2. 2. 2. The crown 
diameter perpendicular to the transect was again estimated with the 
geometric mean of the crown width and length (square root of the 
product) . 


2 . 3 . 1 . 2 . 3 . 2 Photo/Ground Tree Volume Regressions . 

The computed volumes , based on ground measurements, for 
111 trees were matched against the corresponding photo interpretation 
data. Two methods for estimating board foot volume through regression 
were investigated. Method 1 was directly based on the crown diameter 
and height. Figures 2-23 and 2-24 are scatter plots of these varia- 
bles against the ground measured board foot volume. The resulting 
prediction equation was: 


bj^ - -5.564 c -.2785 c^ + .0202 c^ 

-3.036 h + .4715 he -.0137 c^h 56.31 


(2. 3. 1.9) 


where c * estimated crown diameter perpendicular to transect 
in feet. 

h » estimated height in feet. 
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-23. Ponderosa Pine Photo/Ground Comparisons. 
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24. Ponderosa Pine Photo/Ground Comparisons 



2 . 3.1. 2.3.2 — Continued 


2 

The R value waa .787 (multiple correlation R ■ .887) with a standard 
errr£ of the estimate 15.41. Method 2 was a two-step procedure In 
which first the photo crown diameter and height were used to estimate 
ground DBH and height. Figures 2-25 and 2-26 are scatter plots 
showing some of the photo/ground relationships Involved. Then a 
linear function of (DBH) x height was used to predict board foot 
volume as was done In producing the ground-based tree volume esti- 
mates. The resulting equations were: 



ys 

‘>2 

»' (.001352) d^h -5.078 

with 

r2 

s .777 and standard error » 15.39, 

where 

A 

d 

= -1.782 c + .0558 c^ + .00087 c^ 



-.4179 h + .0853 he 



-.00244 c^h + 15.33 

with 


s .859 and standard error = 2.81, 

and 

h 

* .7449 h + .8147 c + .02061 c^ t 10.15 

with 

r2 

= .635 and standard error = 10.67. 


(2.3.1.10) 


(2. 3.1.1: 


These were both acceptable in terms of fitting the photo data to the 
ground data. 

However, in applying these coefficients to the full set of 
1168 photo interpreted trees it was found that the range of crown 
diameters, 3.2 to 42.9, was falling well outside the range In the 
regression data set, 4.6 to 32.5. This was a serious difficulty 
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Pine Photo/Ground Comparisons. 


2. 3. 1.2. 3. 2 — Continued 


because of the presence of quadratic and cubic terms for crown 
diameter In both methods above. Some of the estimates were very 
Inflated (over 500 board feet). To counteract this effect; a simpler 
model of height, crown diameter and their product was used for the 23 
trees whose crown diameters exceeded 30 feet: 

£3 » -2.338 c -2.068 h .1710 he + 31.38 (2.3.1.13) 

Scatter plots of volume estimates against crown diameter and height 
are found In Figures 2-27 and 2-28 for Method 1. There It Is seen 
that near the origin equation 2. 3. 1.9 Is unstable. The correspond- 
ing plots for Method 2 are In Figures 2-29 and 2-30. There the 
instability does not exist. On the basis of this analysis, It was 
decided to use Method 2, where equation 2.3.1.10 Is used for crown 
diameters no greater than 30 feet and equation 2.3.1.13 is used for 
those greater than 30 feet. Those estimates which were negative 
(minimum » -4) were changed to zero. 

2. 3. 1.2. 3. 3 Photo Plot Volume Estimates . 

The photo plot board foot volume per acre estimates were 
obtained from the tree volume estimates in exactly the same was as 
the ground plot estimates (equation 2. 3. 1.3) except that here the 
tree volumes were predicted from photo measurements and 

photo/ground regression coefficients. This resulted in estimates 
for 690 photo plots. 

2. 3. 1.2. 3. 4 Photo/Ground Plot Volume Regressions . 

In order to produce volume estimates the photo and ground 
plots were tied together using a regression estimator. The 27 forest 
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27. Ponderosa Pine Photo/Ground Comparisons 
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Pine Photo/Ground Comparisons. 
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-30. Ponderosa Pine Photo/Groxind Comparisons 


2. 3. 1.2. 3. 4 — Contlriued. 


ground plot estimates were matched with the corresponding photo plot 
estimates and linear regressions were performed. One regression was 
performed for all 27 and also the full set was divided into tvto 
groups for separate regressions. The first group was comprised of 
photo strata 20, with 2 plots, and 21, with 15 plots. The second 
group was the 10 plots from photo stratum 22. Figure 2-31 is a 
scatter plot of the photo and ground plot estimates. Table 2-36 
shows the regression results. The two regression groups differ 
widely in their least squares coefficients. These results were 
applied in a regression estimator which is described l.i Section 
2. 3 .2. 3.3. 


2. 3. 1.3 Evmmary of Landsat Classification for Sample Data . 

The 27 summary categories from the Landsat classification 
were used in the sampling procedure at the PSD level but not at the 
level of selecting plots within PSU's for photo and ground data 
collection. The following summarizes the representation of photo and 
ground plots within the Landsat categories. 

2. 3. 1.3.1 Range Summary . 

Table 2-37 sumiaarizes the number of photo and ground plots 
by Landsat stratum. Note that several categories are grouped 
together - this was done at the sampling stage and so was carried 
through into the productivity estimation. The 135 ground plots were 
located on 57 of the 108 flight lines. The Landsat information was 
combined with the photo and ground data to obtain forage estimates 
by the method described in Section 2. 3.2. 3.1. 
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Table 2-37. Summary of Landsat Classification 
for Range Sample Data 


w 
























0 












•M iH 












0 eu 

o 


<n 


in 

VO 

o 

d 

00 

VO 

r* 

• lO 

M 

r4 







fH 















0 












u 












o 









• 

• 


W 
























0 












•M H 












0 0» 

o 


KO 


VO 

in 

o 

00 

<*» 

o 

Cl 

• 0 
0 

r* 

rs 

OV 

in 


Cl 



fH 

*• 

C! 

00 

o 

H 

Z Q 












A 












eu 
























Tl jC 












o o> 












«M (0 4J M 












O - 0 0 

S 0 h c 

in 



o 

o 


o 


15 

00 

00 

• to H 












0 04 4) U iJ 
Z MO 
























(0 (fl 

0> «0 







rH 




in 

• 

M C 

<« e (« 

VO 


fH 

o 

t 

VO 

m 

• 

fH 

• 

• 

o 

• 

Cl 

• 

■P H 0) 

u 0 

O 0 

<0 

ro 

o 

iH 

r* 

o 

Cl 

r« 

o 

Cl 

in 

n 

VP 

fH 

n 

fH 

B Xt 
























tn 












*0 







c- 




ro 

• 

M C 

(U C (0 

VO 


VO 

H 

• 

o 

• 

Cl 

• 

o 

• 

• 

o 

• 

M M 

o 

in 

r* 

o 

VO 

00 

o 

VO 

VO 

o 

o 

0 0 
< 0 

9> 

<N 

fH 



fH 



00 


CO 

£ 












■M 
















(U 








C 

rH 

•H 

r-4 

•H 


■u 

q 









0 

0 


01 






•9 


4J 

a 

M 

>i 

M 

>1 


0 

s 






3 

u 

jc 


M 

U 

(0 

M 

U 

S 

»0 

c 

0 

CO 


u 

'S 

M 


•9 

•0 


4) 

O' 

q 

(0 

01 

iH 

q 

(0 


Q 

(1) 

(U 


Xi 



a 


S 

\ 

■3 

01 

c?» 

<0 

S' 

0) 

u 

(0 


A 01 
M 0 

ns 

H 


•9 

CO 

CO 

P 

M 


(0 

(0 



4J 

•g 


Vt 

(9 

JS 

(0 


M 

3 

14 

0) 


4J U 

q 


xs 

u 

c/1 

H 


0 

a 


M q 

0 

01 

M 

o 

*0 

u 

? 

4) 

GO 

9$ 

J3 

41 U 

z 

n 

\ 

\ 


1 

(U 

% 

4) 

01 

01 \ 


q 

•a 

•0 

q 

l« 

0) 

0) 

d 

0 

St ^ 

0 

IM 

s 

41 

X 

4J 


4J 


M 

Q 

q 

U) 

a 

S 


0 

0 

0 


J3 

0 

■H 

\ 

#H 

9 

41 

z 


(0 

M 

0> 

q 

,!^ 

>0 01 


A 

01 

\ 

G 

0 

0 

0 

>3 

U 

41 0 

q 

0 

n 

M 

V 

3 


(U 

V 


q 

X 4) 

a 

M 

q 

q 

O' 

M 


h 

VI 

a 

fH 

“H M 

•H 

A 

M 

c 

m 


u 

O 

u 

9 

n 

S U 

B 

M 

xs 

M 

(0 










r' 












Cl 

n 


(Q 









\ 

fH 


M • 

(Q 0 

H 

Cl 

ro 


m 

VO 


00 

f>l 

fH 

\ 

o 

11 

H Z 





1 




\ 

f-1 


u 





1 




O' 





Table 2-37 . — Continued 
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Table 2-37. — Continue 


No. of 
Ground Plots 


135 










No. of 
Photo Plots 

113 

1620 










No. of 
PSU's 
Selected 
for Flight 
Lines 

00 

108 










Hectares 

in 

thousands 

• 

00 

215.3 



1 







Acres 

in 

thousands 

00 

• 

o 

fN 

532.0 










Landsat Class Description 

Upland Desert Shrub/Blackbrus 

Total 










Class 

No. 

26 

1 
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2. 3.1. 3.2 Woodlanct Sununary . 


The number of woodland photo and ground plots by Landsat 
class are summarized in Table In the woodland sampling - 

frame, instead of using Landsat to stratify PSU's the sum of PSU 
pixel counts for classes 14, 15, 16 and 16 was used as the weight 
for unequal probability sampling (see Section 2. 2. 1.2). Two plots 
each from 8 PSU's comprised the ground data set. 

2. 3. 1.3. 3 Forest Summary . 

Table 2-39 shows a summary of the number of photo and 
ground plots by Landsat stratum. The 27 ground plots were located 
on 13 of the 46 flight lines. 
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Table 2-38. Summary of Landsat Classification 
for Woodland Sample Data 



Sage/Mixed Shrub 402 163 
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Table 2-39. Summary of Landsat Classification 
for Forest Seunple Data 


(A 




m 








■P 












0 
























0 o< 

o 

O 

o 


o 

o 

o 

O 

o 

o 

O 

. fO 












0 S 












as 9 












0 












u 












o 












n 




IjHI 

■■ 

■■ 



** 



p 












0 












tp H 












0 P4 

o 

o 

o 




PO 

o 

o 

o 

rH 

• 0 












0 P 












Z 0 












£ 












04 












p 












T) X. 












fl) O' 












tp (fl p -H (A 
0 - 0 ip a> 
D 0) b C 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

• cn H -p 
0 04 0) P 4^ 

z yj 0 












tp 












(A 












0) 

ro 

<N 

CM 

O' 

CM 

m 

r«- 

X 

PO 

o 

cn 

fO 

Jj 

H 

in 

iH 

a 

CM 

rH 

O' 

O' 

o 


r^ 


o 

O' 

iH 

4<f 

CM 



00 

fH 

CO 

0 

0 

X 


% 



fH 

M* 




fH 

CM 

o 

o 










fH 

fH 










(A 

0> 

OJ 

o 

O' 

r* 

M* 

fH 

o 

l>* 

X 

H 

o 

r* 


PO 

VO 

fH 

rH 

M* 

PO 

00 

O' 

o 


00 

fH 


in 

M* 

CM 

CM 

X 

X 

iH 

0 

*4* 


« 

% 


<n 

% 



fH 

CM 

r~ 

VO 

Tl* 

CM 



O 







<N 

CM 

rH 



fH 












\ 










<u 


0) 






c 

H 

iH 


p 


p 






0 

•H 

•H 


2 


2 







0 

0 




m 







CO 

CO 


0 


0 






04 

*r4 

>1 

>1 


0) 

p 


0 

M 








TJ 


u 


O 






Q 

u 

s 


\ 








tA 

0) 

s 

(3 

CO 




% 

P 









p 


3 

"O 





U 

0) 

o 


jC 


P 

c 




JG 

(A 

(A 

fl 

O' 

cr 


X 


A 

3 



3 

3 

(Q 

(0 

(U 



X 

fH 


3 

3 

P 

(A 

(A 

p 

p 



*0 


p 

3 

43 

p 

P 

3 

p 


P 

0 


r: 

P 

X 

O 

0 

3 

04 

0) 


P 

o 

W 

X 

U 


X 

X 


(A 

Jd 

0) 

z 

3 

N 

\ 

•d 

jJ 


(U 

(A 

CA 


3 

<0 

•O 

3 

fl) 

(A 

n 

0) 

<U 

<u 

Q 

3 

0) 

c 

P 

c 

3 

X 

p 

P 

p 


P 

a 

3 

O 

3 


•p 

0 

0 

0 

•0 

D 

CA 

•P 

\ 

fH 

» 

S 

(A 

(A 

(A 

c 

M 

-0 3 

P 


n 

3 

\ 


0 

0 

0 

to 

U 

A) P 

3 


3 

X 

3 

0) 

0) 

0) 

iH 

3 

X U 

04 

P 

3 

3 

P 

p 

P 

P 

a 

fH 

•P fl 

•H 

43 

p 

G 

3 


o 

o 

O 

p 

X 

X u 


X 

o 

X 

X 

(A 

(A • 
(0 0 
H Z 

H 

CM 

(*> 

■ 

m 

VO 

1 

00 

O' 

o 

H 

11 

O 




HI 



■ 




L_J 


2-165 
























Table 2-39. — Continued 


HI 

** 

0 P4 
. X 

0 

o 

o 

— 

o 

o 

<*> 


rH 

o 

o 

o 

IT 

11 

(0 










' 


** 












0 
























0 » 

in 

o 

a 

in 

m 


in 

o 

in 

O' 

IN 





<n 

00 


CM 


rH 


00 

• 0 

0 -P 





H 





rH 














0« 






■ 






■p 












*0 jC 
























tp 0) P >H 0) 












0 - 0 H « 

o 

o 

o 

o 

o 

o 

o 

o 

M> 

in 

r- 

D (U Cu C 










(N 


• W iH -H 












0 04 0) M 












7i (no 












M-l 












(0 












0) 

o\ 

« 

m 

or 

in 

in 

m 

IN 

in 

in 

m 

M 

H 

cn 

VO 


ro 

o> 

o 

VO 

VO 

rH 

M* 


o 


in 

ov 


o> 

m 


IP 

in 

a 

U 

% 



% 

% 

■rH 

H 




rH 

0 

o 


fO 

00 

<T\ 







0 

B 

H 


H 


CN 








CO 

m 

o 

p' 

r* 

a> 

in 

M» 

IN 


in 

(0 

in 

a\ 

IN 

o 

r~ 

IN 

rH 

in 

o 

m 

in 




in 

m 

Mr 

o> 

p* 

rH 

O' 

rH 

un 

u 

% 


% 

% 


M* 

ro 



«i 

M* 




m 

lO 

n 





rH 


< 

m 


m 

O' 

r~ 





rH 


a 












0 












•rt 












P 












a 












•H 












p 




,0 






•d 


u 



0) 

3 






0) 


U) 



O' 

p 






X 


0 



<0 

j: 





10 

•p 


Q 



cn 

cn 





o 

X 





\ 

\ 



H 



\ 


(A 



p 

p 

P 

b 

10 


0) 

0) 

0) 

(A 



0) 

« 

0) 

P 

P 


c 

e 

c 




04 

04 

04 

P 

P 


•p 

•p 

•p 




•H 

•H 

•H 

£ 

<0 

0) 

a 

a 

a 

u 



c 

c 

c 

cn 

a 

p 






n 

3 

9 

3 


10 

9 

10 

10 

10 

4J 


3 

*o 

•9 


c 

x: 

P 

0) 

n 

(0 

(9 


P 

\ 

\ 

''V 

•H 

u 

rH 

0 

0 

0 

CA 



c 

C 

0 

10 


9 

p 

p 

p 



W 

0 

0 

0 

P 

•o 

U 

0) 

0) 

0) 

Q 

0) 

P 


>1 

>1 


(1) 


•0 

•0 

•0 

5 

cn 

-.P 

c 

c 

c 

§ 

X 

p 

6 

c 

c 

M 

(0 

lO 

•H 

•p 

•H 

£ 

•p 

O' 

0 

<5 

0 


CO 

cn 

04 

04 

04 

X 

s 

< 

04 

iP 

a 

(0 












(A • 

CM 

m 


in 

VO 

r~ 

00 

a> 

o 

rH 

IN 

(Q 0 

h_sb 

fH 

H 

*H 

rH 

P 

rH 

rH 

rH 

IN 

IN 

IN 

o 













2-166 






































Table 2-39. — Continued 


No. of 
Ground Plots 

o 

o 

o 

o 

o 

27 






No. of 
Photo Plots 


o 

o 


o 

069 






No. of 
PSO's 
Selected 
for Flight 
Lines 

o 

o 

o 

o 

o 

46 






Hectares 

1071 

H 

42 

8119 

4724 

153,114 






Acres 

2646 

CM 

103 

20,063 

11,673 

378,352 






Landsat Class Description 

Shadow 

Hater 

Bare 

Upland Desert Shrub/ 
Blackbrush 

Sage/Grass 

Total 
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CM 
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2.3.2 


Analyses for Productivity Estimation and Map Verifi- 
cation. 


Pour types of analysis were carried out in the produc- 
tivity estimation and map verification task. These analyses were 
computer class descriptions, acreage estimation, productivity 
estimates and a detailed cost analysis. This work provided 
quantitative estimates of species composition and terrain des- 
criptions for each computer class in the Landsat classification. 
Area by vegetation type and class description verification data 
was also produced. Productivity estimates consisted of forage 
production and grazing carrying capacity for range areas. 

Woodland area productivity was estimated as cubic volume of 
Pinyon and Juniper. Forest production consisted of board 

foot volume for Ponderusa pine. 

2. 3. 2.1 Computer Class Descriptions . 

The 1979 large scale photogiaphy was used to describe 
the 117 computer classes. A portion of the photos were randomly 
selected for use as an independent verification data set. This 
data set consisted of 586 photos selected from the 3000 available 
on the basis of 20% of the photos in classes with at least 30 
photos for the computer class. This criteria was used to help 
ensure that an adequate number of samples would exist to des- 
cribe and verify the maximum number of computer classes. This 
data set was selected when the 27 summary classes were being 
considered as the classes with which to work. As the analysis 
continued, significant problems were found with the original 
groupings and the work shifted to consideration of new groupings 
based on the 117 classes. 
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2. 3.2,1, 1 Analysis of Variance . 


The 1979 large scale photo plot point locations were 
digitized from U.S.G.S. topographic quads provided by the B.L.M. 
The digitized points were then Input to GRNDSPOT, along with 
the ARIZT2 transformation (see Section 2, 1.1. 1.4. 2. 2) and the 
classification Images (27 and 117 classes) to produce ERIS flies 
containing the point Identifier and computer class for that 
location. This data was then paired (In ERIS) with the species 
composition photo Interpretation data (sec. 2. 2. 3. 2). This 
paired data set was then Input to MINITAB to produce estimates 
of the mean and standard deviation of percent cover by computer 
class. This was done for both the description and verification 
data sets and the combined full data set. The 11'78 photo- 
interpretation data was also processed in this 

manner to provide backup data for those classes not adequately 
sampled in 1979. Figure 2-32 shows the computer class description 
menu form; Table 2-40 shows the photointerpreted categories 
that go into making up the class description species groups. 

The computer class descriptions for the sampled classes can be 
found in Appendix 2-G. 

The BLM staff, using the 117 class descriptions, 
placed each class into the Arizona strip vegetation framework 
(see Appendix 1-A ) . This grouping resulted in 14 level 3 
groups or 10 level 2 groups. Table 2-41 shows the assignment of 
the 117 classes to the vegetation framework. 

2. 3. 2. 1.2 Contingency Analysis 

Elevation, slope and aspect information was obtained 
for each computer class by performing a contingency table analysis 
(IDIMS function GOUTABLE) between the Landsat classification 
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Figure 2-32. Computer Class Description Menu Form 


NASA/BLM APT PHASE II 

ARIZONA CLASSIFICATION VEGETATION/TERRAIN MENU DESCRIPTION 

Summary Claaas No. of Photo Samples: I of Total: ■ ■ 

Name: No. of Acres: No. of Hectares: 

Spectral Classes: Z of Area: 


A. Trees Mean 

Ponderosa Pine 

Plnyon Pine 

Juniper 

Other Tree 


C. Sh~ubs — Mountain Mean Std.Err. 

Gambel's Oak 

Turblnella Oak 

Other Shrub _____ 

E. Grasses Mean Std.Err. 

Perennials 

Annuals 

G. Non-Vegetation Mean Std.Err. 

Barren (Rocky) 

Barren (Sandy) 

Water 

Shadow 


B. Shrubs — Desert Mean Std.Err 

Creosote 

Bursage 

Blackbrush 

Big Sagebrush 

Other Shrub 

D. Riparian Woodland Mean Std.Err 

Cottonwood 

Willow 

Other Shrub 

F. Cactus Mean Std.Err 

Yucca 

Other Cactus 


I. VEGETATION 
(X Cover by Species) 

Std.Err. 
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Figure 2-32. Continued. 

ELEVATION COMPUTER CLASS GEOBLOCK 



Z HAS NO ASPECT (FLAT) 






Table 2-40. Class Descrlpticrv Species Groupings. 


Class Description 
Group 


Photo- Interpreted Spocies/Species Groups 
(Number from Table 1) 


Other Tree 


Other Desert Shrub 


Other Mountain Shrub 


Other Riparian Shrub 


Fir (12), Locust (18), Other Three (34), 
Dead Tree (45) 

Apache Plume (1), Arroweed (2), Big 
Rabbitbrush (4), Cliffrose (9), Pour- 
Wing Saltbush (13), Joshua (15), Little 
Rabbitbrush (17), Mesquite-Acacia (21), 
Shadscale (26), Snakeweed (27), 
Turpentine Bush- Jimmy Weed (29), Winter- 
fat (31), Other Desert Shrub (33), Sage- 
like Low Desert (43), Deadbrush (44) 

Manzanita (19), Mixed Chapparal/Mixed 
Mountain Shrub (20), Serviceberry- 
Cecnothus (25) 

Aspen (3), Saltcedar (24), Riparian 
Shrub Ml) 


Other Cactus 


Cholla (8), Other Cactus (35) 
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Table 2-41. Assignment of the 117 Computer Classes to 
the Vegetation Framework. 


Vegetation Framework classification Computer Classes 


1) 

1— 

Agriculture 

15-24, 114 

2) 

211 

Ponderosa Pine Forest 

64, 65 

3) 

311 

Pinyon-Juniper Woodland 

33, 35, 41-43, 46, 47, 
52, 53, 55-59, 61-63, 
66, 83-85 

4) 

322 

Riparian Woodland 

25, 26 

5) 

412 

Upland Desert Shrub 

1-10, 12-14, 74, 75, 
81, 86, 87, 89-100 

6) 

421 

Great Basin Sagebrush 

11, 28-30, 32, 34, 38, 
40, 44, 48-50, 54, 
101-103, 105, 106, 108 
110, 112, 115-117 

7) 

423 

Blackbrush 

88 

8) 

424 

Other Tall Shrub 

111 

9) 

425 

Half Shrub 

45 

10) 

432 

Oakbrush 

60, 113 

11) 

433 

Other Mountain Shrub 

27, 31, 39, 51 

12) 

511 

Perennial Grassland 

36, 37, 77, 104, 107, 
109 

13) 

6 — 

Barren Land 

76, 78-80 

14) 

7 — 

Water 

67-73, 82 


2. 3. 2. 1.2 — Continued 


and the terrain images derived from the DMA digital terrain 
data (sec. 2. 1.1. 6). This analysis was performed on the 27, 
117 and 14 class Images as well as on the 10 and 14 clasfi 
■smoothed (10 acre aggregation) images (see Section 2.3.J). 

The results of this analysis are plotted on the class descrip- 
tion menu forms for the smoothed 14 class Images (see Appendix 
2-G) . 


2. 3. 2. 2 Acreage Estimation and Verification of Map Accuracy. 

2. 3 „ 2. 2.1 Acreage Estimation . 

The allotment and pasture boundary digitizing was 
overlaid on the Landsat classification with the function POLYSUM 
to produce pixel counts by pasture and computer class. 

Regression analysis was used to produce adjusted acreage esti- 
mates. The dependent variable was the proportion of large scale 
photo (LSP) plots in each flightline in a stratum based on photo- 
interpretation. The independent variable was the proportion of 
LSP plots in each flightline in a stratum based on the pixel by 
pixel Landsat grouped (27 classes) classification. Figure 2-33 
shows an example of the regression analysis output. These 
regression equations were then applied to the pixel counts by 
class to produce the adjusted acreage estimates (see Appendix 
2-H) . 


2. 3. 2. 2. 2 Verifleabion of Class Description by Anova . 

The verification data set was analyzed the same as 
described above to produce a second set of descriptions. These 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


2. 3. 2, 2. 2 —Continued . 

descriptions were then manually compared to the descriptions 
obtained from the other data set. This comparison turned up no 
discrepancies in assignment of a class to the level 3 categories. 
The verification descriptions are presented in Appendix 2-*I. 


>REGn PIPROP I L8PR0P 
THE REGRESSION BRUATION IS 
Y»< .509 >X1 

ST. DEV. T-RATIO ■ 
COLUMN COEFFICIENT OF COEF. COEF/S.D. 

XI LSPROP .869380 .102674 4.99 

THE ST. DEV. OFiY ABOm' REGRESSION LINE IS 
3 • .179146 

WITH ( 26- 1 ) > 19 DEOIEES OF FREEDOM 

ANALYSIS OF VARIANCE 


DUB TO 

DP 

SS 

MS-SS/'DF 

REGRESSION 

1 

.799166 

.799166 

RESIDUAL 

19 

.669729 

.3269I6E-61 

TOTAL 

26 

1 .46889 



Figure 2-33. Example of Regression Analysis for 

Acreage Estimation. 
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2 . 3 . 2 . 3 Productivity Estlmatea . 

A, 11 the productivity estimation tasks described below were 
performed using the Earth Resources Inventory System (ERIS) with the 
exception of a few hand-held calculator tasks. 

2. 3. 2. 3.1 Range Forage Estimates . 

The primary parameter to be estiamted was air-dry available 
forage biomass from palatable species and estimated for full utiliza- 
tion (e.g., if plant is 20% utilized, the measured weight is multiplied 
by 1/(1- .2)). The secondary parameter was air-dry forage biomass 
at time of inventory, whether usable or not. These estimates were 
required by Landsat strata, and by allotment and pasture. Confidence 
intervals were required by Landsat strata. A description of the 
reduction of the ground and photo sample data is given in Sections 

2. 3. 1.1.1 and 2. 3. 1.2.1. The following section describes how the 
Landsat information was integra.ed with the photo and grouxid data 
through the sample design framework to yield final estimates. 

2. 3.2. 3. 1.1 Forage Estimates P e r Unit Area . 

Both Landsat and photo stratification were used in locating 
the rangeland ground plots, as described in Section 2. 2. 1.2. There 
were 135 ground plots located in 57 PSU's. 

It was originally envisioned to use a straightforward 
stratified two-stage estimator. However, there were two features of 
the design which affected the applicability of this estimator. First, 
the Landsat PSU stratification was overlapping. This resulted in 
stratum estimates containing components which were "shared” with other 
strata and therefore could not be simply added together to form total 
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2 . 3 . 2 . 3 . 1.1 


Continued 


estimates. Second , the selection of ground plots was made across all 
PSU's from within the photo plot strata (and then additional plots 
are selected from the same PSU's represented by the Initial sample). 
The second stage of sampling Is In Itself a double sample of photo 
and ground plots. Therefore the subsampling Is not performed li^de-* 
pendently within PSU's as Is generally assumed with multi-stage 
estimation. Furthermore the second stage Is a double sample for 
only some of the selected PSU's and Is a single sample for the 
remainder. 


These complications In the sampling method suggested 
modifications to the estimator. It was decided that available theory 
and methodology did not sufficiently associate these exact sampling 
frames and selection probabilities with an applicable estimator. 
Therefore « approximations were pursued which still reflected the 
most important featureti cf the design. One simplifying aspect was 
that flight lines were sele ;ted with equal probability within Landsat 
strata. Also, the total number of pixels within each Landsat stratum 
was available as well as the Landsat class value for each photo and 
ground plot. 

These considerations resulted In the following intermediate 
estimators for forage biomass across all Landsat strata based on 
double sampling for stratification: 
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2. 3. 2. 3. 1.1 


Continued. 


“1 

Y <■ Ib/acre forage estimator 
M « 

Yv 

* H gm ^ Lt ^hijro 

^ (ij) es 

■" m 


—1 

Y « total lb forage estimator 

-1 ~1 
Y^ » A Y-^ 

A 

where “ Y^j from eq. 2.3.I.2. 

(ground plot forage for Landsat stratum h, 

PSU i, SSU j and photo stratum m) 

» set of ground plots from photo stratum m 

g^j^ * number of ground plots from photo stratum m 

» number of photo plots from photo stratum m 

Jl = total number of photo plots 

M * miinber of photo strata 

A = total number of acres in sampling frame 

For the kilogram per hectare results, the appropriate valuer were 
substituted for and A. Also the forage estimates without the 

availability and utilization factors applied were obtained by substi- 
tuting the appropriate Y^ijm* 
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2 . 3 . 2 . 3 . 1.1 


Continued 


The Landsat stratum intermediate estimates used the follow- 
ing post-stratified ratio-type estimator: 



Ib/acre forage estimator for Landsat stratum h 


y 



( 2 . 3 . 2 . 3 . 2 .) 


At this point an adjustment was made to correct for the error in 
application of the utilization factor (see Section 2. 3. 1.1.1. 2) . The 
final Landsat stratum estimates were then computed to perform this: 


— —2 

Yh - yVu^ 


where Uj^ = average utilization for all plants in Landsat stratum h 


= total lb forage estimator for Landsat stratum h 

_ A 

■ ^h ^ 


where 


^h 


'h 

H 


average value of stratum h 
number of acres in Landsat stratum h 
number of Landsat strata 
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2. 3. 2. 3. 1.1 


Continued 


The final estimates across all Landsat strata were then computed 
as follows: 


y 


I 


A 



Y - Y/A 

These estimates are recorded in Table 2-42 for the case with avail- 
ability and utilization factors applied and in Table 2-43 without 
the factors applied. 

Standard errors and 80% confidence Intervals were also 
computed for the overall estimates and the estimates by Landsat 
strata. This again required some simplifications but did reflect 
the main features of the design and estimation processes. 


The overall variance estimates were computed following 
the double sampling for stratification framework. The variance 
estimator is adapted from Honijn (1 ) , p. 177. 

VAR (Y) - ^ [Tj^ '^2 " *^3 ’ ’’ 4 ] (2. 3. 2. 3. 3) 

where Tj - ^ J ^ (y„ - Y) 

m 

_ 2 2 

T » §5* 

2 A ^ ^ 2 gm 



RANGE FORAGE ESTIMATES BY LANDSAT STRATA 
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RANGE FORAGE ESTIMATES BY LANDSAT STRATA 
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(1) INSUFFICIENT SAMPLES Paa» 2 of 3 
































RANGE FORAGE ESTIMATES BY LANDSAT STRATA 
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(1) INSUFFICIENT SAMPLES Pw 3 of 3 
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2.3«2.3.I.l 


Continued 


T, 


A-1 1 V Ira Sm^ 

A r Z* ^ 


in 


- 


1 - 


1 I 


sra 


ra 


wherei Sin 


2 




y « average of the y's within photo etratum m 
■* m 

This estimator assumes proportional allocation which was not entirely 
true in this case but was not drastically different than the actual 
allocation. Even though the samples were clustered by being located 
within PSU's, no correction was used because the sample size was very 
small within PSU (usually two ground plots) . This meant that the 
effect of clustering was mild. 

The variance estimates by Landsat strata used the wlthin- 
strata sample variances: 


VAR(Yj^) 




km 


Since is large, this is approximately: 


VAR(Yh) - 


( 2 . 3 . 2 . 3 . 4 ) 


The standard errors (square root of variance estimates) and relative 
standard errors (standard error divided by estimate) were computed 
and are summarized in Tables 2-44 and 2-45. 
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RANGE FORAGE ESTIMATES BY LANDSAT STRATA 
STANDARD ERROR SUM1ARY 
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Table 2-44. — Continued. 

RANGE FORAGE ESTIMATES BY LANDSAT STRATA 
STANDARD ERROR SUMMARY 
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Table 2-44. — Continued. 

RANGE FORAGE ESTIMATES BY LANDSAT STRATA 
STANDARD ERROR SUMMARY 

PALATABLE SPECIES ONLY AMOUNTS AVAILABLE AND PROJECTED FOR FULL U 
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Table 2-45. 

RANGE FORAGE ESTIMATES BY LANDSAT STRATA 
STANDARD ERROR SUMMARY 


Table 2-45. 

RANGE FORAGE ESTIMATES BY LANDSAT STRATA 
STANDARD ERROR SUMMARY 

\BIE SPECIES ONLY CURRENT AMOUNTS. AVAILABLE OR UNAVAILABLE 
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RANGE FORAGE ESTINATES BY LANOSAT STRATA 
STANDARD ERROR SUMMARY 
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2 . 3 . 2 . 3 . 1.1 


Continued 


Confidence intervals at the 80% level were found using 
the normality assumption. The degrees of freedom were 9]^~1 ^or the 
stratum estimates and g-1 ■ 134 for the total estimates. These 
results are found in Tables 2-42 and 2-43. 

The forage biomass estimates by allotment and pasture were 
found by applying Landsat stratum estimates to the pixel counts 
within allotments and pasture. 

A 

y ■ Ib/acre forage estimator for pasture p of allotment a 
ap 

^ap - Z ’h 'Slap / Z "Slap 

h-1 h-1 

A 

Y » Ib/acre forage estimator for allotment a 
& 

^ H ^ H 

"^'a • I I 

h*l h*l 

where ^ap “ acres in pasture p of allotment a with 

Landsat class h 

s number of acres in allotment a with Landsat class h 

When an allotment or pasture had acreage for a Landsat svunmary class 
which had no estimate, both the nximerator and denominator sximmation 
did not include that class. Again, the analogous computations were 
made for the case without the availedsility and utilization factors 
applied and both versions were converted to kg/hectare. The estimates 
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2. 3. 2. 3. 1.1 


Continued 


by allotment are found in Tablet 2-46 and 2-47. The estlmatet by 
allotment and pasture are found In Appendix 23, 

2. 3. 2. 3. 2 Woodland Estimates For Plnyon Pine and Juniper . 

Cubic foot volume estimates for plnyon pine and juniper 
were computed by combining the photo and ground plot data with the 
Landsat Information. The photo plot estimates were computed from 
individual tree volume predictions based on regressions with ground 
measured volumes » as described In Section 2. 3. 1.2. 2. The PSU volumes 
were computed from these In the following manner: 

.15 

Ai Y 

Yi - ^ 2* Yij (2. 3. 2. 3. 5 

j»l 

where « area in acres for PSU 1 (usually ^.22. A) 

Yij » estimated volume for photo plot j of PSU 1 from 
eq . 2 . 3 . 1 . 8 . 


Since the PSU's were selected with unequal probability, the 
selection weight, were used in the estimator. 


Y 


total cubic foot volume estimate 

IT I 

1-1 


( 2 . 3 . 2 . 3 . 6 ) 
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Table 2-46. 

RANGE FORAGE ESTIMATES BY ALLOTMENT 
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Table 2-47. 

RANGE FORAGE ESTIMATES BY ALLOTMENT 
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2 . 3 . 2 . 3.2 


Continued 


7 - cubic foot volume per acre eatlmate 

-s N 

-V/J 

1-1 

N 

where P. - W. / Z W. 

^ ^ 1-1 ^ 

- number of pixels from Landsat eummary categories 
14, 15, 16 or 17 in PSU 1 

N » number of PSU's in population with a positive 

It should be noted that those PSO*s with W^-0 had no chance of being 
selected for photo and ground sampling. These comprise 41 of the 388 
PSU’s and cover 5405 of the total of 65,358 acres. So those 5405 
acres were excluded from the sampling process by the Landsat class!' 
fication because of no apparent woodland cover types. 

The estimates by the three allotments (Wolfhole, Wolfhole 
Mtn. and Blackrock) were made in two steps. First, estimates by 
Landsat summary classes were found as follows: 

A 

*h ■ " (2. 3. 2. 3. 7) 

lyh'S. 

« total cubic foot estimator for Landsat class h 
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2 . 3 . 2 . 3. 2 


Continued. 


*h - V*h 


■ cubic foot per acre estimator for Landsat class h 


where 7^ ■ average photo predicted volume for Landsat class h 

> number of acres from Landsat class h 

Second, the estimates by allotment were computed using the Landsat 
class per acre values and the Landsat class breakdown within allotment. 

A 

» cubic fopt volume per acre estimator for allotment a 
H H 

* I ^h ^ha/ J ^ha (2. 3. 2. 3. 8) 

h“l h*l 

where A. » number of acres in allotment a with Landsat class h 
na 

These estimators were applied separately for pinyon pine and juniper 
and the results were converted to cubic meters per hectare. The 
pinyon pine results are in Tables 2-48 and 2-49. The juniper results 
are found in Tables 2-50 and 2-51. 

The variance for the overall estimates was computed as 
follows (Raj (2), p. 49): 
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Table 2-48. 
PINYON ESHMf^TES 
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Table 2 .9. 

PINYON ESTIMATES IN METRIC UNITS 


ORIGINAL PAiuL IS 
OF POOR QUALITY 
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(32,204:62.762) 
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24.263 

TOTAL 

AREA IN 
HEaARES 
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26,450 
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NUMBERS IN PARENTHESES ARE BOX CONFIDENCE INTERVALS. 
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Table ^ il. 

JUNIPER ESTIMATES IN METRIC UNITS 
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2. 3.2. 3.2 


Continued. 


VAR(Y) 





2 


) 


A 

VAR(Y) - VAR(Y)/A^ 


N 

where A ■ J A^^ 
i-1 


(2. 3. 2. 3. 9) 


The standard errors and relative standard errors were found frost these. 


The variances by allotment were computed in two steps. 
First, the variances for standard post-stratification were used 
(Koni jn (1) , p. 185) . 


A 

VAR^(Y^) 




(2.3.2.3.10) 


where A area in acres for allotment a 
a 

■> number of photo plots in Landsat class h 
» sample variance within Landsat class h 

2 I <yhij - 

^ (i,j)e h 


m 


2 . 3 . 2 . 3. 2 


Continui^d. 


The term used a» an estimate for the number of photo 

plots In allotment a within Landsat class h, which was unknown. 
Equation 2.3.2.3.10 reduces tot 


VAR^(y^) 



h 



(2.3.2.3.11) 


The second step was to reflect the fact that these sample data were 
clustered by means of the two-stage sampling. The effect of this Is 
to somewhat Increase the variance of the estimates. The amount of 
the increase can be estimated by using the following equation 
from Konljn (1) , p. 310 because the subsampling fraction is small 
(15/360 - .042): 

VAR(t ) - VAR^{^ ) 11 «»• (m-l)5I (2.3.2.3.12) 

Cl cl 

Where m » number of photo plots per PSU (15) 

6 " intracluster correlation coefficient 

The quantity (S is a measure of how homogeneous the measured resource 
attributes are within PSU's. If 5 is high, the effect of clustering 
is strong because little additional information is added by sampling 
many plots within the PSU. The value of 6 was estimated as .378 
for pinyon and .453 for juniper from the data. These values are 
fairly high and gave multiplying factors of 2.51 and 2.71 for the 
standard errors. 
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2. 3. 2. 3. 2 — Continued. 


The confidence intervalis at the 80% level were computed 
based on these standard errors and the normality assumption. -The 
degrees of freedom used are summarized in Table 2-52. The allotment 
degrees of freedom were estimated as follows (see Cochran (3), p. 95): 


degrees of freedom 



(2.3.2.3.13) 


where I 


ha 


, ^ha 
h 

a 


u 


ha 



ha 




The resulting confidence intervals are presented along with the 
estimates on Tables 2-48 and 2-49 for pinyon pine and Tables 2-50 
and 2-51 for juniper. 

The average size and age are also reported on the pinyon 
and juniper volume tables. These are the simple averages from those 
trees on the ground plots which had the measurements recorded. 
Confidence intervals are also given for those values. 

2. 3. 2. 3. 3 Forest Board Foot Volume Estimates for Ponderosa Pine . 

The ponderosa pine estimates were obtained by using the 
photo/ground plot regression estimator along with the photo and 
Landsat stratifications. The estimator by Landsat stratum was: 
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Table 2-52. 

Allotment 
Wo If hole 

Wolf hole Mountain 
Blackrock 

Overall 


Degrees of Freedom for Plnyon Pine 
and Juniper 80% Confidence Intervals 


Estimated 
Degrees of 
Freedom 


Plnyon 

Juniper 

86 

77 

106 

121 

182 

231 

600 

600 


Two Sided 
t-dlstrlbutlo‘n 


Value 


Plnyon 

Juniper 

1.292 

1.293 

1.290 

1.289 

1.285 

1.284 

1.282 

1.282 
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2. 3.2. 3. 3 — Continued 


■ board foot volume per acre estimator 

" X ^ ^ ‘’m(^hm ■ X j] 

m»l " 

^ (2.3.2.3.14) 

Jl vt” " 

where * munber of photo plots Interpreted as stratum m from 
Landsat stratum h 

A 

bj^ « linear regression coefficient for on X^ijm 

A 

^m “ constant regression coefficient for on 

s number of photo plots from Landsat stratum h 

Xj^ = average of photo plot estimates by Landsat stratum 

h and photo stratum m 


Xjjj = average of photo plot estimates for those plots 

which also have ground plots and which are in photo 
stratum m 

= average of ground plot estimates from photo stratum m 

A A 

The coefficients a , b_ are taken from Table 2-36. For those photo 

mm 

plots which were not in any of photo strata 20 , 21 or 22, the com- 
bined regression coefficients were used. 
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2 . 3 . 2 . 3 . 3 — Continued 


f 


Just as In the woodland case, the Landsat forest 
stratification was designed to focus the sampling in the areas of 
interest. Only those PSU's which had Landsat summary class 20, 21 
or 22 occurring were possible flight line locations. This reduced 
the sampling frame to 87,015 acres from the population size of 
378,351 acres <441 of the original 2063 PSU's). 

Some Landsat strata which had positive acreage in the 
sampli: ^ frame had no samples. Four of these had insignificant 
acreage (less than 10 acres) and were assigned zero volume. The 
other two had 99 and 138 acres and were assigned the weighted average 
for all the non-ponderosa pine strata. 

The overall estimates were then computed: 

^ H 

^ “ Z ^h ^h (2.3.2.3.15) 

h=l 



where » area in acres of Landsat class h 

The Landsat strata and overall estimates are in Table 2-53. 

The forest area was comprised of two separate, disjoint 
areas, called A and B. Estimates were made for each of these areas 
by Landsat strata in the same way except that coefficients for 
regressions through the origin (no constant term) were used to add 
stability for the smaller sample sizes: 


( 
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PONDEROSA PINE VOLUME ESTIMATES 
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OTHFR NUMBERS IN PARENTHESES ARE 80% CONEIDENCE INTERVALS 



























































PONDEROSA PINE VOLUME ESTIMATES 
AREA: A & B COMBINED 



(1) NO SAMPLES. INSIGNIFICANT ACREAGE - ZERO VOLUME ASSIGNED 

(2) NO SAMPLES. SIGNIFICANT ACREAGE - POOLED ESTIMATE USED FROM ALL STRATA EXCEPT 20. 21. 22 

(3) MOT ESTIMABLE - INSUFFICIENT SAMPLES 

OTHER NUMBERS IN PARENTHESES ARE BOX CONFiOENCE INTERVALS 



2 . 3 . 2 . 3 . 3 — Continued 


• board foot volume estimator for area t, Landsat 
nc 

stratum h 



where ^hmt " photo plots in Landsat stratum h, photo 

stratum m and area t 

» number of photo plots in Landsat stratum h and area t 

^hmt “ average of photo plot estimates by Landsat stratum h, 
photo stratum m and area t 

» least squares regression coefficient for model with- 
out constant term 

Weighted sums were again used to obtain the estimates across all 
Landsat strata. The estimates for Landsat strata 20 , 21 and 22 were 
adjusted to produce additivity over the two areas, that is: 



ht 


? ^t 


ht 


This assures that the area A and B estimates add to the overall 
estimates for those three strata. The estimates by Landsat strata 
are recorded in Table 2-54 for area A and Table 2-55 for area B. 
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PONDEROSA PINE VOLUNE ESTIMATES 
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PONDEROSA PINE VOLUME ESTIMATES 



(1) NO SAMPLES. INSIGNIFICANT ACREAGE - ZERO VOLUME ASSIGNED 

(2) NO SAMPLES. SIGNIFICANT ACREAGE - POOLED ESTIMATE USED FROM ALL STRATA EXCEPT 20. 21 

(3) NOT ESTIMABLE - INSUFFT:IENT SAMPLE: 

nruro miimocdc tn paocMTurcrc aar rnNRnFNr.F INTERVALS 


Table 2-55. 

PONDEROSA PINE VOLUME ESTIMATES 
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PONDEROSA PINE VOLUME ESTIMATES 
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(3) NOT ESTIMABLE - INSUFFICIENT SAMPLES 













































2. 3. 2. 3. 3 — Continued 


The variances were estimated under the model of double 
sampling for stratification applied to domain estimation# as found 
in Konijn (1)# p. 179. 


where 



VAR 


estimated variances for Landsat stratum estimates 



(2.3.2.3.16) 




2 


And the variances across Landsat strata then were: 


VAR 




The same method was used for variances by area. These results 
were adjusted for the effect of clustering In the seune manner as 
described for woodland in equation 2.3.2.3.12. The relevant param- 
eters were : 




Intracluster 

Correlation 

Coefficient 

Standard 

Error 

Multiplying 

Factor 

Area 

A 

.156 

1.784 

Area 

B 

.319 

2.338 

Overall 

.343 

2.409 
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2. 3.2. 3.3 — Continued. 


The degrees of freedom for the confidence intervals were computed in 
the same manner as in the woodland estimation (equation 2.3.2.'3.14) . 
The standard errors, relative standard errors and 80% confidence 
intervals are included with the parameter estimates in Tables 2-53, 

2-54 and 2-55. 

2 . 3 . 3 Costs Analysis . 

The costs associated with this project can be separated 
into recurring and non-recurring components. Only recurring cofts 
will be considered in this analysis. In other words, "false" starts, 
development of procedures and non-special management costs have been 
omitted to the extent possible. The objective of this analysis is to 
associate a cost per acre with both the vegetation mapping and the 
productivity estimation elements of thi project. Table 2-56 sum- 
marizes this analysis and shows that at the average price per ESL 
manhour over the life of the project, the mapping effort for the 
2.2 million acre site cost 3<: per acre. The productivity estimation, 
performed on 550,000 acres, cost 16.3^ per acre given the mapping has 
been performed. These figures are based on the costs associated with 
achieving the objectives and accuracies of this project. If similar pro 
jects were undertaken but with higher accuracy and precision specifi- 
cations, it should be anticipated that the costs per acre could be 
greater. 

2.3.4 LANDSAT Digital Products . 

The LANDSAT digital hardcopy output products consisted of 
4 application maps and 2 sets of color-coded classification images. 

The application maps were generated from the digital classification 
quantitative Inventory and ancillary data in the form of digital 
terrain data and administrative data. The map products were: 
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Table 2-56. Phase II Cost Summary 


Project Element 

Total ESL _ 
Man-Hours 

Mapping ^ 
Component 

Estimation 

Component 

Additional 
> Costs 

Landsat Processing 

570 

450 

120 


Multistage Sampling 

270 

120 

150 

4800^ 

Data Collection 

Photo Acquisition 

50 

50 

- 

11,300^ 

Photo Interpretation 

700 

-CO 

200 


Ground Data Collection 

100 

- 

100 

31,970^ 

Productivity Estimation 

1250 

240 

1010 

700^ 

4 

Output Products 

Map Output 

180 

180 

- 


Tabular Summaries 

120 

40 



Totals 

3240 

1580 



1660 

$48,770 

Vegetation mapping cost 

» [ (1580 hrs 

X $3. 17/hr.) 

+ 11,700) 



f 2,2 million acres « 3<!/acre 

Productivity estimation » (Given the mapping is completed and avail- 
able) 

= ((1660 hrs X $31. 70/hr) 37,470) 

i 550,000 acres ■ 16.3<:/acre 


1. RSRP, U.C. BerJceley, in survey planning model runs and estimation 
woric. 

2. Acquisition of 3000 large scale aerial photography plots. 

3. Collection of ground data on 179 ground plots. 

4. Does not include generation of final report. 
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2.3.4 


— Continued. 


1) Rangeland Suitability defined aet 

1) current production of ueeable forage above 
20 Iba. per acre, 

11) located within a four mile radius of water, 
ill) slope less than 51 percent. 

2) Potential Rangeland Suitability defined ast 

1) current production of useable forage above 
20 lbs. per acre, 

ii) service area of water is greater than four mile 
radius, 

ill) Slope less than 51 percent. 

3) Sagebrush Treatment Area defined as: 

1) Great Basin Sagebrush Vegetation conununity. 
ii) Slope less than 15 percent. 

4) Fire Flash Fuel Areas defined as: 

i) Annual grass and forbs cover 30 percent or 
greater of ground. 

11) Aspect southwest - south • Southeast. 

These products were presented with a mapping minimum of 10 acres 
(4 hectares) at a scale of 1:126,720 as a halftone grey-coded map 
on stable base transparency material. Maps 1 and 2 were limited to 
the 10 allotments sampled for range productivity. Maps 3 and 4 
represent data over the full project area. 
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2.3. 4 


— • Continued. 


The 5~ map output products were for the forest administrative 
area A and consisted of 8 x 10 inch Polaroid hardcopies of the color- 
coded computer class groupings. These products were presented at 10 
acre (4 hectare) minimum mapping units for the Level II framewor)c and 
for the Level III framework groupings developed by the BLM Arizona 
Personnel. An IDIMS computer compatible tape (CCT) was generated 
with this data for the full project area as well. 
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3.0 


EVALUATION OF THE PRE-INVENTORY PLANNING FOR P.J. 


The pre-inventory planning model runs were made using the 
first five parameters shown in column 1 of Table 3-1. The planning 
model estimated that 13 photo transects (PSUs) containing four photo 
plots (SSUs), with two of these visited on the ground » would be 
required to achieve a t20% allowable error eight times out of ten 
over a 300,000 acre planning unit with a 20% crown closure of P-J. 

The budget for this effort was estimated to be $2192.40 (Table 3-2, 
column 1) . Using the actual parameters (Table 3-1, column 2) obtained 
from the inventory, the model estimated that 63 photo PSUs with one 
SSU per PSU with 14 ground visits would be required to achieve the 
same results. See Appendix 3-A for a more complete description of this 
planning model. 

The difference in the estimated parameter and the actual 
parameter (Table 3-1, column 3) accounts for the differences in the 
budgets. The effects of the changes in each of the parameters is 
shown in Table 3-3. 

The largest change in budget comes from the change in intra- 
cluster correlation. This parameter affects the value of additional 
plots taken on a transect. In the original design the maximum number 
of SSUs per PSU was arbitrarily set to 4. However, the high intra- 
cluster correlation (.597) indicates in the model that aided plots 
in the PSU are of very little value to the inventory. Therefore, 
the system design allows only one SSU per PSU. 

The next most significant parameter was the coefficient 
of variation. It was estimated at .75 in the pre-inventory planning — 
by simulation of P-J populations and sampling this population in the 
computer. The error in the estimate of cv probably came from two 
sources: 1) the simulation did not vary the relationship between 

crowa diauneter and tree volume, and 2) the distribution of the trees 
over the simulated stands was random rather than clusters as th^ 
usually are in nature. 
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Table 3-1. 

Pre and Post Inventory Estimation of Parameters 
Used in the Inventory Planning Model 



Pre- Inventory 
Estimate or 
Assumption 

Actual 

■1 

Landsat vs Photoground 
Correlation 

.50<2) 

.58 

+ .08 

P-J Photo Plot Volume vs Ground 
Estimated Volume 

.B5<2> 

.748 

-.168 

3 

P-J Volume per Acre (ft ) 
Coefficient of Variation 


.93 

+ .18 

Average Number of P-J Trees/ 
Transect 

3.74 t 2.3^^^ 

5.6 

+1.9 

Intra-Cluster Correlation for 
P-J Volume 

oil) 

.597 

+ .597 

3 

P-J Volume per Acre (ft ) 

293.6 ± 48^^^ 

241 

-52 

Relative Standard Error of the 
Estimates 

8.2^^^ 




Values based on 150 foot transect with an average spacing of 
1368 meters. 


(1) No data available 

(2) Based on other related studies 

(3) Estimated using Resource Inventory Services planning model 
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Table 3-2. 

Comparison of Pre and Post Inventory Planning Model 
Runs for P-J Volume Estimation 



Pre-Inventory 
150 ft. Transect 
with 483 ft. 
Spacing 

Post- Inventory 
150 ft. Transect 
with 446 ft. 
Spacing 

Number of Plot PSUs to be 
Interpreted 

13 

63 

Number of Ground Visited PSU 
from the Photo PSUs 

2 

14 

Number of SSUs per PSU (Photo 
and Ground) 

4 

1 

Photo Acquisition Budget ($) 

391.85 

898.00 

Photo Interpretation Budget ($) 

125.93 

271.00 

Ground Data Collection Budget ($) 

1674.72 

3984.00 

Total Budget ($) 

2192.49 

5154.00 


Parameters set In the model are shown In Table 3-1. 























Table 3-3. 

Estimated Sampling Budget for Achieving a t20% AE 
at the .8 Probability Level beginning with the Pre- 
Inventory parameter from Table 1 and Progressing to 
the Actual Parameter. 



Budget ($) 


Original Budget Estimate 

2192.00 


Landsat to Photo Ground .58 not .50 

2018.67 

-174 

Photo Volume to Ground Volume Cor- 
relates .748 not .85 

2160.93 

+ 142 

Intra-Cluster Correlation of .597 
not 0 

3846.99 

+ 1686 

Coefficient of Variation .93 not .75 

4818.87 

+ 972 

Number of Trees per Transect 5.6 
not 3.7 

5153.64 

+ 335 
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3.0 


— Continued. 


The other errors were within the expected errors in the 
model and simulated runs. 

The P-J inventory could have been conducted in several ways. 
Pour possible methods are 

1) a ground based inventory 

2) Landsat classification followed by ground data 
collection 

3) Large scale photographs followed by ground data 
collection 

4) Landsat classification, large scale photography 
and ground data collection. 

Using the parameters from the results of the inventory, the planning 
model was used to estimate the sampling cost for the four methods. 

To compute the cost effectiveness ratio for each method 
the cost of establishing the sampling frame (Landsat classification 
cost or manual interpretation) must be added to the cost shown in 
Table 3-4. If these costs are assumed to be equal, the cost effectiveness 
ratios relative to the least expensive (Landsat, photo, ground) are 


X) 

Landsat) Photo, Ground 

1 

2) 

Landsat Ground Only 

1.1 

3) 

Photo Ground Only 

1.32 


4) 


Ground Only 


1.48 



Table 3-4. 

Sample Requirements and Expected Cost for P-J 
to Achieve a ±20 9 .8 Probability Using 
the Post Inventory Estimates of the 
Planning Model Shown in Table 2 - 2 , Color. n 2. 


COST EFFECTIVENESS OF 
THE DATA SOURCE 

GROUND 

ONLY 

LANDSAT 

GROUND 

ONLY 

PHOTO 

GROUND 

ONLY 

LANDSAT 

PHOTO 

GROUND 

Number of Photo PSU 


0 

0 

92 

63 

Number of SSU/PSU 


1 

1 

1 

1 

Number of Ground PSUs 


38 

25 

22 

14 

Photo Acquisition 

($) 

0 

0 

1171.95 

898.09 

Photo Interpretation 

($) 

0 

0 

395.67 

270.95 

Ground Data Collection 

($) 

7618.66 

5677.41 

5244.34 

3984.60 

TOTAL COST 

($) 

7618.66 

5677.41 

6811.91 

5153.64 
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— Continued. 


If these results hold over the entire study site (2,500,000 
acres) and the average planning unit is 300,000 a?res, total total savings 
in sampling would be $20,541.00. 

It is beyond the scope of this phase of the APT to conduct a 
complete cost effectiveness analysis. However, if the cost of Landsat 
classification and sample frame development were compared to the cost 
of conventional photo interpretation and sample allocation methods, 
the cost effectiveness ratios might be even greater in favor of the 
Landsat , large scale photo and ground method. 
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APPENDIX 1-A 


VEGETATION FRAMEWORK 


The following pages detail the hierarchical structure 
of the vegetation framework classification as developed by 
BLM personnel from the Arizona Test Site and the Denver Federal 
Center. Also included are the narrative descriptions of the 
framework categories. 
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I. Water 

II. Silty Water, Clear Water 

III. River, Stock Reservoirs 


I. Barren Land - Less than aX Vegetative Cover 

II. Natural Features 

III. Exposed Surface Rock - gravel size to large 
rock and cobble lying on the surface. 

IV. Rock - Five CM or greater In diameter 

IV. Gravel - Less than five CM In diameter 

IV. Gravelly soli - mixture of surface soil and 

gravel. Gravel generally contiguous 

with small Interspaces of base soil. 


III. Other Exposed Rock - Exposed bedrock, lava 
flows, cinder cones, talas slopes, alluvial 
sand and gravel. 


III. Bare Soil - Soil with no rock or gravel present - 
sand dunes or flats, salt flats, other bare 
soil areas. 


II. Cultural Features - Man-made barren areas. 

III. Roads & Road Cuts ,. Mines Gravel Pits, 
dry stock ponds. 
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I. 


ORIGINAL PAGE IS 
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Agrlculturil Land 

II. Cropland & Pasture 

III. Cropland - field crops, row crops, 
orchards, etc. 

III. Pasture - irrigated permanent pastures 


I. Forest - tree crown closure is lOZ or more of total area. 

Under story species may be pinyon and/or juniper 
and/or shrub communities or non-exlstant . 

Generally found at elevations above 6500 feet. 

II. Coniferous Forest - Majority of trees present 
are conifers. 

.III. Ponderosa Pine Forest - Majority of trees 
present are Plnu s ponderosa species. 

IV. Open Ponderosa Pine - crown closure of 

mature trees is 66Z or less of total area. 
Understory species may Include any of the 
following: Juniper ( Juniperous 

osteospemia) ; pinyon ( Plnus e dulis 
and P. monophylla) ; big sagebrush 
( Artemisia tridentata) or other sagebrush 
species; gambel's oak ( Quercus gambelii) ; 
shrub live oak ( Q. turbinella) ; chaparral 
( Artostaphylos and Garyig spp.); cool season 
perennial grasses; and secondary growth 
ponderosa pine. 


1-A-ll 


0RIG\^^^ 
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IV. Closed Pondeross Pine - Crown closure of 

mature trees is greater than 66!L. Understos 
species generally non-existant excispt secom 
i * growth ponderosa pine. 


111. Mixed Conifer Forest - Approxima' uly equal 
amounts of ponderosa pine and other conifer 
species, generally true fir ( Abis spp.) and 
Douglas fir ( Psuedotsuga roensy-ii) . Very limited 
in distribution, mostly at highest elevations 
and/or north facing slopes. 


IV. Ponderosa Pine - Fir - Mixture of ponderosa 
pine and true fir and/or Douglas fir. 


I-. Grassland - ground cover is primarily herbaceous and is 

five X or greater of total vegetation cover 
but usually less than 507*. May be characterized 
by varying amounts of half-shrubs and tall shrubs 
Annual forbs usually present; amounts depending 
on current growing conditions. 

II. Plains Grassland - grasses may be annual or perennial, 
short or mid grasses. Occurs generally above AOOO 
feet elevation. 


l-A-12 


1X1. Warn Season Perennial Grassland - perennial 
grasses less than 30 cn tall. Comniunlty co> 
dominated by warm season species of the genera 
Hilaria end Bouteloua . Occurs at elevations of 
AOOO to 5500 feet. 


Calleta - Crama Steppe - broad valleys and 
rolling hills with at least 5% perennial 
_ grass cover. Association dominated by 

ORiQlNAL PAGE IS 

OF POOR QUALITY galleta ( Hilaria iamesil) , and/or black 

grama ( BouteloUa eriopoda) . Blue grama 
(B. gracillisl may also be a dominant or 
co-dominant bn some sites. Subdominant 
species include burro grass ( Scleropogon . 

e 

brepffolius) and sand dropseed ( Sparobolus 
cryptandrus) . Cool season grasses, Indian 
ricegrass ( Oryzopis hymenoides) . Needle 
grasses ( Stipa spp.). and squirreltail 
( Sitanion hystrix) . may also occur. 

Varying amounts of shrubs are scattered 
throughout the type. Half- shrubs, 
snakeweed ( Xanthacephelum spp.), little 
rabbitbrush ( Chrysothamnus viscidflorus) 
and wlnterfat (C eratiodes lanata) are quite 
randomly distributed throughout the associ- 
ation. Tall shrubs, yucca ( Yucca baccata) and 
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Brlghan tea ( Ephedre spp.) occur mainly on 
hilleidea and rldgetopa. Shrubi comprise 
lass than total cover. Annual forbs 

are scattered throughout. Russian thistle 
( Sal sola spp.) can occur in extensive stands 
In overgrazed bottomlands In this association. 

IV. Galleta - Crama - Shrub Steppe • 

* Same as preceding class except that 

shrubs, half-shrubs and tall shrubs 
. may form up to of the. total cover 
of the area. 

Cool Season Perennial Grassland - 
Perennial grasses up to 100 cm tall- 
community dominated by cool season grasses of 


the genera Sltanlon , 


L» St li5a and 


can be found at various elevations 


but Is most common above 4300 feet. This 
community Is widely scattered but limited 
in size. Thought to be the remnants of a 
much larger and diverse cool season grass 
community that has been reduced by over-grazing 
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IV. Squirreltail - Cool season grassland 

dominated by squirreltail ( Sitanion hystri> 

IV. Squirreltail - Shrub - Squirreltail 
association w.fth up to 5X of ground 
cover provided by Shrubs. 


3V. Keedlegrass -> Cool season grass associa- 
tion dominated by Stipa species. 

IV. Western Wheatgrass - Cool season grass 

f 

association dominated by Agropbron s mithii . 

3V. Indian Ricegrass - cool season grass 
. association dominated by Aryzopsi4 
hymenoides . 

IV. Mixed Cool Season Grasses - Cool season 
grass association but with no one 
particular species dominant. May also 
include up to 5X cover by shrubs. 

III. Annual Grassland - Annual grasses usually less 
than 30 cm tall. Twenty percent or greater 
of total cover is provided by annuals, primarily 
Bromus species. Community has usually become 
established in a shrub or woodland community 
after a fire. Occurs generally at 3000 to 
SOOPfeet elevation. 
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IV. Cheatgrasi - 20% to 1007. of cover of 

area la provided by B. tectorlum and/or 
B. rubeora . Six veeka feacue ( Festuca 
octoflora) may alao be present. 


IV. Cheatgrass - Shrub - same as preceding 
except that shrubs provide up to 5% of 
total vegetation cover of area. Shrub 
species may be any of those found in 
the test site. 

III. Introduced Seedings - Vegetation treatment 
areas dominated by introduced cool season 
mid grasses, 30^100 cm tall, usually vheat- 
grasses ( Agropyson spp.). Russian wildrye 
( Elymus cinereus) may also be present. 
Community occurs above 4500 feet elevation. 

IV. Mixed Wheatgrass - Mixture of seeded 
wheatgrass species, usually Agropyron 
desertorum , A. spicatum , A. pubescence . 
May also Include russian wildrye, Elymus 
jdn etis and mountain brome, Bromus 
ii- i^: -nis . 

■>. . II , 1 111— 1» 


IV. Mixed Wheatgrass - Shrub - As above 
but with tall and half shrub species, 
gambel oak, live oak, sagebrush, little 
rabbitbrush and snakeweed reinvading. 
Shrubs less than 5% of total cover. 
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IV. 


Annual Rye - pure community of annual 
Secala cereale. Usually small 


areas, up to AO acres in size. 

I,* Woodland - Crown closure of trees (excluding forest species) 

is 10% or greater. Tree species primarily ever- 
green, piriyon and Juniper; also included decidu- 
ous species: aspen (Populus tremulides) ; tree 
farms of gambel oak; and riparian species, 
cottonwood ( Populus fremontii) . and salt cedar 
( Tain’arix spp . ) . 

11. Evergreen Woodland - evergreen species, pinyon and/or 
Juniper, comprise majority of woodland cover. 

111. Pinyon-Juniper Woodland - Majority of woodland 
cover is comprised of pinyon. Juniper or 
mixture of both. Pure pinyon stands are jammCaa 
rare.. Pure Juniper stands and mixtures of 
pinyon and J\iniper are equally common. Tree 
heights up to 6 m for Juniper and 10-12 m for 
pinyon. The community covers a large area 
of the test site. Occurs at elevations above 
AOOO feet. Under story species generally shrubs: 
sagebrush, cliff rose, oak, blackbrush and 
man^anita. Sagebrush is most common. Grass 
species also occur, usually blue grama, galleta 
or squirreltail. 
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IV. 


Juniper - Crown cover 10% or greater 
Juniper species. Other tree species 

origu^m. i'-; ■ 

OF POOH less than 10% cover. Onderstory species, 

shrub or grass, less than 5%. 


IV. Juniper-Shrub - Crown cover primarily 
Juniper, 10% or greater. Shzrb species 
comprise at least 5% of total cover. Most 
common shrubs are sagebrush, blackbrush, 
cliffrose, gambel oak, rabbitbrush, 
snakeweed. 

IV. Juniper -Pinyon - Crown’ cover comprised of 
both Juniper and pinyon with Juniper 
. making up at least 50% of tree total cover. 
Pinyon crown cover is 10%-50% of total 
tree cover. 

IV. Juniper- Pinyon- Shrub - Same as preceding 
with at least 5% of total cover, comprised 
of understory shrubs. Shrubs are same 
as those listed under Juniper-shrub 
association. 

IV. Pinyon - Crown cover primarily pinyon 
species, 10% or greater. Other tree 
species less than 10% cover. Under story 
shrubs or grass less than' 5% cover each. 
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Pinyon-Shrub - Crown cover of pinyon 
lOX or greater. Shrub species comprise 
at least 5X of total cover. 

PlnyoXi- Juniper - Crbwn cover comprised 
of both pinyon and Juniper with pinyon 
making up at least SOX of total tree cover 
Juniper cover is 10X-5QX of total tree 
cover. 

Pinyon-Juniper-Shrub - Same as preceding 
with at least 5X of total cover comprised 
of understory Shrubs. 

II. Deciduous Woodland - Majority of woodland cover made 
up of deciduous tree species. Some species, 
i.e., gambel oak,LSalt cedar and locust occur 
in both tree and shrub forms. Trees are 
defined as single stems or multiple stems 
arising from same basal point and attaining 
height of 3 meters or greater. 

III. Upland Woodland - Woodland types not associated 
with water or wet areas. Occurs primarily 
in mountainous areas above 5,000 feet in 

t 

small, scattered, rather isolated patches. 

Under story shrubs such as sagebrush, shrub 
oaks, cliff rose may occur. 


IV. 




IV. 


OF POO^ 


IV. 
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IV. 


Cambel Oak - Majority of crown cover 
made up of tree form of gambel oak. 

Trees usually slngli: stemmed and may 
occur in clones or small groves. Other 
tree species may or may not be associated. 

IV. Gambel Oak-Locust - Crown cover made of 
varying percentages of tree forms of 
oak and New Mexico locust. Only known 
occurrence is on Mt. Turnbull area. May 
be associated with ponderosa pine 
and/or aspen. 

IV. Aspen - Majority of crown cover made up 
of aspen. Very limited distribution . 
Occurs only as small, widely scattered, 
isolated patches above 6000 feet. 

Riparian Woodland - Tree species associated 
with perennial or intermitten streams or ponds. 
May be found at all elevations within test site. 

IV. Cottonwood - Crown cover 107. or greater 
cottonwood trees or saplings. Other tree 
species, if present, comprise less than 
10% crown cover. Sh rubs comprise less 
than 5% ground cover. 


IV. Cottonwood • Shrub - Crovm cover 10% or 
greater cottonwood species, but with at 
5% shrub under story. Shrubs may include 
arroweed ( Pulcea spp.) and shrub forms 
of salt cedar ( Tamar lx pentandra) and 

page t ’ 

OF POOR QUAUT’t willow ( Sallx spp.) 


IV. Salt Cedar - Crown cover 10% or greater 
of tree forms of salt cedar. Other tree 
species. If present, comprise less than 
10% crown cover. Shrubs comprise less 
than 5% ground cover . 

IV. Salt Cedar-Shrub - Crown cover 10% o^ great 
of salt cedar trees with at least 5% shrub 
under story. Shrubs Include shrub form of 
sal t cedar . 


I. Shrubland - Shrubs comprise at least 5% ground cover. 

Includes tall shrubs up to 3 meters and half 
shrubs up to .5 meters. Includes both hot and 
cold desert types. Some species, notably 
snakeweed, lyclum, yucca, and little rabbitbrush 
occur at all elevations in many vegetation types. 
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Mohave Desert Shrub - "hot desert" shrub types. 
Occur generally at elevations up to about 3500 
feet although some species range up to 6000 feet 
depending on soil, vegetation community and 
rainfall. 


111 . 




Riparian Desezc Shrub - cover dominated by 
those species associated with intermittently 
flowing desert washes. Taller species such 
as cottonwood and tree forms of salt cedar 
may be present as scattered individuals. 

In most cases total vegetation cover in this 
type will not exceed 10-20^. 


IV. Salt Cedar - Arrowweed - Cover dominated 

by shrub forms of salt cedar and arrowweed 
( Pluchea sericea) . 

IV. Willow - Cover dominated by willows, either 
Salix s pecies, true willow, or Chilopsis 
species, desert willow. or both. 

IV. Mesquite-Acacia - Cover dominated by 

species of Mesquite ( Prosopis) , Acacia 
( Acacia) or both. Generally less than 3 M 
in height although some species may exceed 
3 m. 
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IV. Rabbitbrush - cover dominated by the 

tall, greater than 1 m in height, species 
of rabbitbrush. 


III. 




Upland Desert Shrub - cover dominated by 
those shrub species not associated with 
riparian areas. This type is found 
generally on mesas, bajadas, hills and • 
other desert uplands. Elevatlonal limit 
limit generally 3500 feet, although there 
is much variation, depending on the site. 


IV. Creosote - cover is dominated by creosote 
( Larrea tridentata) . Other vegetation 
species may or may not be associated. 

IV. Mixed Creosote-Bursage - Creosote and 
Bursage ( Ambrosia dumosa) are co- 
dominants with approximately equal 
amounts of each in the vegetation cover. 

IV. Bursage - cover dominated by bursage. 

Can from pure stands to diverse 

mixtures with other species. Perennial 
grasses may also be present In the 
under story. 
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IV. Yucca - cover dominated by yucca species. 
This would include Spanish bayonet (Yuccf 
baccate) , Joshua tree ' (T. bretffolia) , 
and Mohave yucca (Y. Shidigera) , cover 
can be comprised of any one or a combina* 

a 

tion of these species. Other desert shru 
species are usually associated. 


IV 


)WG»NAL PAGE iS 
)F POOR QUALITY 


Mixed Desert Shrub • cover comprised of 
several different species in approximate!; 
equal amounts. No one species is domin- 
ant. Also can be applied to those areas 
that don't fit the above associations of 
desert shrub. Species likely found in 
this class would be lycium, snakeweed, ^ 
rabbitbrush, range ratany, ( Kramerla 
parirfolia) , spiny hop sage ( Graya 
spinosa) , cholla and other cactus type 
species. 


II, Great Basin Desert Shrub - "cold desert" shrub types. 
The major component of this type is Great Barin 
Sagebrush. Occurs at all elevations above 3500 feet 
in the test site. May contain scattered patches or 
individuals of Mohave Desert Shrub. 
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PftGE IS 
QUAUTf 


Great Basin Sagebrush - cover is dominated 
b> one or more of the various species of 
sagebrush ( Artemisia) , usually big sagebru 
tridentata) . Consounity can occur as 
pure stands, with grass landerstory, with 
other shrubs, with pinyon and/or Juniper 
(less than 10% crown cover) or as any 
combination of these. 

IV. Big Sagebrush - occurring as a pure 
stand or with other understory or 
shrub species forming less than 
5% cover. 

• 

IV. Big Sagebrush > Perennial Grass: 

Big sagebrush comprises the dominant 
cover with perennial comprising 5% 
or greater of the cover. 

IV. Big Sagebrush - Mixed Shrub: 

Big sagebrush comprises the dominant 
cover with other shrub species forroi 
5% or greater of the cover, as a sin 
species or in aggregate with other 
species. 


IV. Big Sagebrush - Tree: Big sagebrush 

comprises the dominant cover with tree 
species (pinyon* juniper, ponderosa pine) 
forming 5X - 10% of the cover. 


III. Saltshrub cover dominated by great basin 
saltshrub species, almost exclusively four- 
wing saltbush ( Atriplex canescens ) and shad- 
scale ( A. confertifolla) . Large expanses of 
saltshrub, primarily shadscale, occur only in 
the clayhole area on the east side of the test 
site. Occurs elsewhere as intermitten scattered 
patches, from .^OOOfeat to 6000feet elevation. 


ORIGJNAL 
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IV. Shadscale: cover dominated by shadscale. 

Cover usually does not exceed 10-157.. 

Can occur as pure stands but often with 
russian thistle (Salsola kali). 


IV. Fouxwing Saltbush: cover dominated by 

fourwing saltbush. More scattered than 
shadscale. Occurs at greater elevational 
range, 3000 feet to 6000 feet and with 
a variety of other vegetation. species, 
including grass and shrubs. 
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111. Blackbrush - cover dominated by blackbrush 


( Coleofcyne ramosslslma ) . This species forms an 
ecotone between hot desert and cold desert 
species in the test site. Occurs at 3000 to 5000 
feet elevation. May be in pure stands or in 
association with perennial grass, other shrubs 
or pinyon and Juniper. 


IV. 
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Blackbrush: cover dominated by blackbrush. 

May be a pure stand or up to 5% grass or 
other shrub iipecies, or pinyon and Juniper. 
May- b o a - pur e o t - cwd or up t o 5 % grao s er 

shrub— a p e-e ico i — o r p i nyon a nd junipc g^. 


IV. Blackbrush - Tree: cover dominated by 

blackbrush but with STL to 10% pinyon or 
Juniper. 


IV. Blackbrush - other desert shrv'iib: cover 

dominated by blackbrush but with other 

shrub species comprising more than 5%' 

cover, either as single species or in 

aggregate. 


III. Other Tall Shrub: cover dominated by cold 

desert species other than those listed above. 
These types are generally rather small and widely 
scattered. 
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IV. Cliffrots: cover dominated by cliffrose. 

Uaually aaeoclated with other shrub species, 
grasses or pinyon- juniper which comprise 
less than 5X cover as Individual species 
or in aggregate. 

IV. Great Basin Riparian Shrub - consists 

almost exclusively of apache plume 

( Falleegia paradoxal and bi^ rabbitbrush 

ORIGINAL PAGE IS species ( Chryothamnous spp.). Occurs 
OF POOR QUALITY , ^ 

along Intermittently flowing washes at 

elevations generally above AOOO feet. 

IV. Winterfat - cover dominated by winterfat 
( Cefatoides lanata) . Occurs from lowest 
elevations up to 6000 feet, but as an 
extensive/ie ml y in the Hurricane Vclley 
grassland area. Associated with peren nial 
grass and/or snakeweed and little rabbitbrush 

III. Halfahrub - shrub species less than .5 meter in 
height. Applies primarily to snakeweed and 
little rabbitbrush. Occurs from lowest eleva- 
tions up to 6000 feet. Associated usually with 
other shrubs, hot and cold desert, and/or grass- 
land. 

0 

f 

i 
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IV. 
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IV. 


II. Mountain Shrub - Tall shrub occurring in mountainous 
areas above 5500 feet elevation. Primary species 
include turbinella oak and gambell# oak ( Quercus 
turbinella a nd ^ ga mbelii) ; Manzanita ( Arctostaphylus 
SPP*)* serviceberry ( Amelanchler spp.); and ceanothus 
( Ceanothus spp.). 

III. Mountain Chapparal - Cover comprised primarily 

of manzanita and silk tassel { Garry gl fluYescens) . 
Turbinella oak often associated. Sagebrush in 
small amounts may also be present as « ell as 
scattered pinyon and Juniper. 


Snakeweed - cover dominated by snakeweed 
species ( Xanthocephalum spp.) with little 
rabbitbrush often occurring as a codominate 
or subdominate. Extensive stands occur 
in the grassland type, but also occurs in 
all shrub types. Probably most widespread 
species in the test site. 

Little Rabbitbrush - cover dominated by 
little rabbitbrusb ( Chrysothamnus spp . ) . 
Snakeweed often assoicated as co-dominant 
or sub-dominant. Occurs in much the same 
area as snakeweed. 
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IV. Manzanita - cover dominated by manzanlta sped 
Other Bpeciee, if present, comprise less 
than STL each. 
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IV. Mixed Chapparal - cover dominated by 
chapparal species but with other shrub 
species comprising greater than 5Z cover. 
Pinyon or Juniper can form up to lOX crown 
cover. 


III. Oakbxrush cover consists primarily of shrub 

oak, gambel's oak and/or turbinella oak. Other 
shrub species, sagebrush, scrvicebefry , manzanita 
and cliffrose may also be associated. ' 




IV. Gambel's Oak - cover dominated by gambel 

oak with other species, if present, forming 
lesr than 5% cover either as a single 
species or in aggregate. 

.IV. Mixed oak - cover dominated by one or both 
oak species. Associated shrubs comprise 
5% or greater cover either as single 
species or in aggregate. 


III. Other Mountain Shrub - consists primarily of 

ceanothus and serviceberry , usually associated 
with oakbrush or manzanita. 
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IV. Servlceberry • Ceanothus - cover dominated 
by servlceberry, ceanothua or both. Other 
species forming less than 5%, either singly 
or in aggregate. 

IV. Mixed Mountain Shrub - cover consists 

primarily of serviceberiy , ceanothus, with 
• • 

at least 57. cover formed by other species. 


ORIGINAL PAGE IS 
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APPENDIX 1-B 


PROJECT PLOW DIAGRAM 


The following flow diagram illustrates the order of 
the retailed tasks and their interrelationships with other tasks 
that were necessary to completely define the approach to the 
Phase II technology demonstration. 
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APPENDIX 2-A 


PHOTO/GROUND PLOT SELECTION PROCEDURES 
(Resource Inventory Service) 
January 1979 


INTRODUCTION 

Over the past year, there has been a significant amount of new 
research in the areas of range, pinion juniper, and forest photo 
interpretation and ground data collection. This work has been completed 
by the U.S. Forest Service, the Bureau of Land Management, and a number 
of universities. To ensure that the BLM ASVT program utilizes the most 
cost-effective procedures available to date in Arizona, a literature review 
and rework of the photo interpretation and ground data collection 
procedures was done. 
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OBJECTIVES 


There were four major objectives considered when the literature 
review and rework of the photo Interpretation and ground data collection 
design was done. 

(1) Complete an Inventory for the demonstration area that will 
meet the needs of the resource manager. 

(2) Demonstrate the utility of Landsat data In the Inventory 
procedure. 

(3) Evaluate the potential of large-scale aerial photography as a 
cost-effective element In the Inventory procedure. 

(4) Provide adequate data to allow the design and optimization 
of an operational remote sensing-based Inventory. 

During the literature review, a number of facts were uncovered: 

(1) The volume of pinion and juniper can be predicted rather 
precisely from measunjments of crown area.^*^®^' 

(2) The error In the estimation (on the ground or on ptotos) of 
basal area, cords, and timber volume for individual trees Is 
proportional to the size of the estimate.. 

(3) The error In estimation of forage production (In the range 
environment) Is proportional to the size of the estimate of 
forage production 

(4) The correlation between estimates made on the ground and estlmi^tes 

made from the photo range from 0.8 to 0.95 for most resource 
parameters pertinent to this study 2» 4, 13, 14, 15, 16) 
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(5) The correlation between dot count of % cover and ocular 
estimates of percent cover Is extrerjiely hlqh (0.95).^^®^’ 

(6) The ratio of cost of photo samples to the cost of ground samples 
Is a significant factor In the optimization function 

Small changes In the cost of photo Interpretation are very 
significant when the cost of photo Interpretation Is low. 

RECOMMENDED PROCEDURES 

Given the above set of references and the general references 
associated with the project In the past, the following recommendations are 
made: 

(1) For range resources: the primary photo Interpretation will be 
the ocular estimation of percent cover by species or cover 
type using a % cover comparison key. 

(2) In the forest and woodland strata the trees for photo Interpretation 
and ground data collection should be selected proportional to 
crown width after the plots have been selected from the flight 
line. 

(3) The photo Interpretation and ground data collection should be 

done In several phases: a) photo Interpretation for stratification 
over the entire study area; b) the ocular photo Interpretation 
of the forest, range and woodland; c) a subset of those plots then 
visited on the ground. 
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SAMPLING LEVEL 


There are three sampling levels proposed: 

(1) Using the UCB planning model, primary sample units will be 
selected. This would be a stratified random sample optimized 
for variance, strata size, and expected resource value. This 
sample selection would incorporate administrative boundaries, 
digital terrain data, Landsat computer classification results, 
1978 large-scale photo interpretation results, and the results 
of interpreting existing resource photography. The data would 
be used in the simulation model, allowing the estimation of 
resource values, resource type, the variability of the resource 
and the cost of data collection. After the simulation and 
stratification and computation of sample size, the primary 
sample unit for large scale aerial photo acquisition would be 
selected at random by strata with replacement. 

(2) Using the results of he manual photo interpretation for 
stratification and the detailed information from the planning 
model run, the PSU's and SSU's for ground visit will be selected. 
This selection is completed at this point -n the inventory to 
allow the ground visits to be completed during the upcoming field 
season. If time permitted, this selection would be more 
efficiently accomplished after a detailed interpretation of the 
large scale photography followed by an objective selection of 
plots for further PI and ground data collection. The plots 
selected for ground visit should be based on cover type, 
physiographic position and administrative unit. These samples 
will be restricted to forest samples in the administrative 
forest area, range samples in the range administrative area, 

and woodland samples in the woodland administrative area. 

Forest, range and woodland plots encountered outside the pre- 
scribed administrative boundaries will not be ground visited. 


2-A-4 



Using the same Information as used In the selection of 
ground plots, a portion of the photo plots will be selected 
for ocular PI. All plots selected for ground visit will be 
included in the ocular estimate set. 


PHOTO INTERPRETATION METHODOLOGY 
Photo Interpretation for Stratification. 

(1) Without regard for details of scale, place each of the large 
scale photos into broad strata based on the percent cover 

of trees, shrubs, grasses, litter and nonveg. 

(2) Identify the dominant species or cover type category. 

(3) Place the plants in the dominant species into one of three 
categories: 1) small, 2) medium or 3) large (optional). 

(4) Determine the spatial distribution of the plants in the plot 
into three categories: uniform, mottled or clumped (optional). 

(5) Identify the homogeneity index for each of the plots as 
described in Table 2-26, page 2-98. 


INTERPRETATION FOR PRECISION PLOT SCALE ESTIMATION 
ON ALL PLOTS THAT HAVE BEEN GROUND VISITED 


1. Locate the photo scale points marked by the ground crew on the 
stereo pair. 

2. Measure and record the distance between the points on both photos 
In the stereo pair. 

3. Calculate the average of the two measurements and use this average to 
calculate the difference between actual photo scale and estimated 
photo scale. 

This photo scale measurement should be done Independently of all 

other Interpretation for the plots that were ground visited. 

WOODLAND PHOTO INTERPRETATION 

1. Compute photo scale. 

2. Compute the photo length of the transect. 

3. Locate the transect/plot template on the photo with the center of 
the stereo coverage as the transect center, and the transects 
oriented parallel to the long axis of the photo. 

4. Beginning at the north or west end of the transect, number the trees 
whose crowns are Intersected by the transect. Assume the upper edge of the 

transect line is "razor thin" to minimize the bias Introduced by a 
"wide" transect. 

5. Record the species and condition of each of the numbered trees. 
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6. Measure and record the length of the crown (along the long axis of 
the crown) and the width of the crown (perpendicular to the long axis) 
for all of the numbered trees. 

7. Ocular estimate (using percent cover key) the percentage cover for 
each of the required species or cover types in the plot. For 
computer class description only, (optional) 

8. Estimate the average stand height class for the pinion and juniper 
(optional ). 

RANGE PHOTO INTERPRETATION (OCULAR) 

1. Compute the photo scale. 

2. Overlay the appropriate plot template (rectangle covering the same 
area as the five transects). 

3. Estimate and record the percent cover of ths required cover type 
and species using the percent cover key and type identification key. 

4. For each of the identified species, identify and record the height 
category. 

5. For each of the identified species identify and record the foliar 
density category. 
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FOREST PHOTO INTERPRETATION (OCULAR) 


1. Compute the photo scale. 

2. Overlay appropriate transect plot template and height estimation key. 

3. Beginning at the north or west end of the transect, number and record 
the species of each tree Intersected by the transect. 

4. Measure and record the diameter of the long axis and the diameter 
perpendicular to the long axis for each of the numbered trees. 

5. Estimate and record the height of each of the numbered trees on the 
transect (five height categories). 

6. Count and record the number of standing dead trees 
Intersected by the transect. 

7. Measure and record the width of each dead tree. 

8. Estimate and record the percent cover of each required species or 

cover type (10/^ Increments) for computer class description only, (optional) 
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Transect Length for Forest and Woodland Photo Interpretation and Ground 
Data Collection. 


Using the results of the photo interpretation for stratification, 
each of the forest and woodland plots will be placed into one of three strata 
based on their ranges of percentage crown closure oy species. A transect 
length will be selected for each percentage crown closure strata to provide 
an optimum number of trees per transect. The length will be selected 
based on an expected coefficient of variation by line length and cost of 
extracting data from each line length for each percent cover strata. 


Gro und Data Collection for Range Resources. 

The ground data collection will follow as closely as practical the 
SVIM ground data collection procedures including the ratio of ocular estimates 
to clipped plots and characterization plots to other plots. On each of the 
SSU's (photo plots) five transects will be laid out in the format described 
by U.e. Berkeley. Each of the transects will contain 40 points and four 
weight estimate plots, systematically located along the transect. This will 
provide a total of 200 points and 20 one-tenth meter square plots. On each 
of these points and plots the SVIM procedure will be followed in detail. 

The point and plot data will be summarized using the algorithms currently 
in use by the Bureau of Land Management. The weight and characterization 
plots will be every tenth point on the transect, beginning with the fifth 
point on the first transect. 

In addition to the procedures in the SVIM manual the following 
steps will be taken on the ground: 

(1) The name of the estimator and recorder will be entered on 
the data collection form to allow the development of a ratio 
estimator for each individual providing ocular estimates of weight. 
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(2) The Individual makins the ocular estimates will not know the 
plots selected for clipping prior to making the ocular 
estimates on all plots on the SSU (ground plot). 

(3) The plots to be clipped will be selected at random from the 
ocular estimate plots with one plot taken from each transect. 

(4) Two points, at least 20 meters apart, appearing in stereo, will 
be located on the ground and on the photo and pinpricked on the 
large-scale photography. The level distance between the two 
points will be measured to 1 centimeter and recorded on the form. 
Using an Abney or similar device, the slope to 1 degree will be 
measured between the two points with the direction and degree 

of slope recorded on the form. 

Forest and Woodland Ground Data Collection. 

On the photo plots selected for ground visit, the threes whose 
crowns were Intersected by the transects will be ground visited and measured 
using existing BLM tree measurement, recording and data reduction procedures. 
If It Is determined on the ground that a crown that appeared to be one tree 
from the photo was a clump of trees, or contained obscured trees of adequate 
size to be measured, these trees will be added to the list of trees for 
ground measurement. The clump of trees that appeared as one crown on the 
photograph will be Identified on the field form to allow the data to be 
summarized Into one unit of volume for regressing between photo estimate and 
ground estimate by clump of trees. 

Calculations. 

In addition to the procedures outlined by UCB, five major elements 
should be Incorporated In the estimation procedure: 

(1) The parameter for selected prediction equations should be 
estimated using the Information obtained from the photo and 
ground data collection effort from this project. Existing 
regression coefficients will only be used as a guide. 

I 
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(2) The cost of extracting data, processing the data, and the 
correlation between each of the attributes, imist be 
thoroughly documented to meet the objective of providing 
information for the design of future remote sensing based 
Inventories. 

(3) The error In scale estimation between the large scale 
photography and ground must be Incorporated In the forest 
fnd woodland data reduction procedure because of the need 
to estimate area rather than percent cover. This should be 
done by relating the ground measurement between points to 
the photo measurement between the same points. 

(4) To estimate the error In area estimation from dot counts, 
the procedure described In BLM's Resource Inventory Notes 
should be used rather than the normal binomial error 
estimators. 

(5) Regression equations should be developed that relate 
Individual tree measurements taken from the photo and volume 
estimates taken on the ground. These should Include: 

a) crown width to ground volume for pinyon 

b) crown width to ground volume for Juniper 

c) crown width and tree height to ground volume for pinyon 

d) crown width and tree height to ground volume for juniper 

e) crown width and height class to ground volume for 
Ponderosa pine 

f) crown width and tree height to ground volume for 
Ponderosa pine. 
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APPENDIX 2-B 


DESCRIPTION OF THE U.C. BERKELEY SURVEY PLANNING MODEL 

PURPOSE AND METHOD 

Extracted from the 1979 Forest Resource Inventories 
Workshop Proceedings, Colorado State University, 

Ft. Collina, pp. 983-992. 


BACKGROUND 


THE SURVEY DECISION PROBLEM 


origiwal page JS 
^OOR QUAl.in. 


Figure 1 lllustretes the steps that must be followed In selecting 
the "best" sampling design to be used In e particular sample survey 
problem. However, the basic Cecislon problem focused upon In the frl low- 
ing discussion has two basic components: .(I) given a specific feasible 
design, what Is the optimum allocation of sampling resources and (2) 
given the optimum allocations for two or more designs, which one should 
be chosen? 

Optimum Sample Allocation 

In a multiresource Inventory, optimization Involves dividing up of 
the survey effort to achieve the "best possible results," I.e. , when 
(1) the Information desired Is obtained and (2) the valu^ of some ob- 
jective function it minimized. 

Two principal approaches have been used In defining the Information 
requirements of the mul tlresource survey and each approach has Its own 
objective function to be minimized. These are (1) minimize the cost of 
the survey subject to the constraints on the minimum precision that Is 
required for each parameter estimated and (2) minimize the cost plus 
tome loss function that is related to the precision that is obtained for 
each parameter estimated. 

These approaches have been discussed by Hazard (197A, 1977) • 

Hamilton (1976), Burkhart et al . (1978), and others. In the development 
that follows^ the first approach (minimize cost subject to constraints) 

Is assumed. (However, the first approach could also use the "cost plus 
loss" criteria If appropriate loss function^ could be formulated.) 

• 

For each alternative sampling design, then, we wish to choose the 
combination of sample sizes that minimizes the total variable cost 
subject to the following constraints: (I) The variance, of the estimator 
for each parameter is not more than the maximum specified, and (2) the 
sample sizes are not more thun the population sizes for each stratum or 
less than some minimum size for degrees-of-freedom considerations. Since 
the variance functions are usually nonlinear functions of the sample 
sizes, this problem becomes a nonlinear programming problem that can be 
solved for the sample sizes, using procedures discussed below. 

Choosing the Optimum Design 

After the optimum combinations of sample sizes are determined for 
each design, the designs are then ranked by cost. If all factors have 
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been considered in formulating the cost function and the constraints, 
then the design with the lowest cost should be selected as optimum. 
However, many factors are difficult to include in the design problem 
formulation. For example, Aldred (1371) found that fixed costs in terms 
of equipment availability, staff expertise and experience, and possible 
measurement errors had to be considered in the final selection of the 
"best" sampling desilgn. 

Finally, remember that this ranking of alternatives is based upon 
approximate co^ts of a survey that has_ not yet been ruiv^ _ Those w^^ 

famll lar with unexpected’cost bverr uns w i 1 1 oot for the 

of the designs that are ranked highly. It is often on the question of 
risk that more com(.’ex designs are bypassed for designs that use ratio 
or regression estimav^ors. 

THE SURVEY PLANNING MODEL 

With the decision problem formulated above, three problems remain: 

(1) to develop a representation of the mui ti resource population to be 
sampled, (2) to use this representation and other available data to de** 
velop the necessary cost and variance relationships, and (3) to obtain 
the optimum allocation of survey effort over the various parts (strata, 
stages, etc.) of the design. 

To solve this problem, a series of computer programs were developed 
by the staff of the Remote Sensing Research Program (RSRP) , University 
of California. These programs, referred to collectively as The Survey 
Planning Model (SPM) > follow the concept initiated in the univariate 
case by Aldred (I97i). That is, a spatially-referenced population "rep- 
resentation" Is constructed and the sampling process is simulated to 
develop the necessary variance and cost information required for a 
minimum-cost allocation of effort for alternative sampling designs. 

Detailed description of the theory and applications of the SPM are 
given by Titus ^.1378), Titus et al. (1376, 1377), and Thomas and DeGloria 
(1973). Figure 2 presents the SPM program structure. The individual 
programs developed as part of this model are described below. 


Program SIMULATE 


The first program, SIMULATE, produces a spatially referenced popu- 
lation of resource parameter values appropriate to the inventory problem 
at hand. Two inputs lire required by SIMULATE. First, a cell-by-cell 
spatial class assignmr*,nt file is required representing a cover type map 
for the survey area registered to some (X, Y) coordinate system. 

This file can be developed from digitizing existing maps or strata 
on aerial photographs, or from classification of digital Landsat, 
aircraft, ar.d/or terrain data. The second input to SIMULATE is a set of 
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rtsource porameter estimates (mean vector and covariance matrix) for 
each cover type. The parameter set chosen should represent parameters 
of ultimate interest on the ground (e.g., forest volume, growth, useble 
forage production) or parameters strongly correlated with these when 
estlMtes of ground parameters are not available. These parameter esti- 
mates will normally be obtained from previous surveys, pilot surveys, or, 
at lest resort, from educated guesses. 

SIMULATE assigns to each cell in the spatial class file a vector of 
resource parameter values. These values are chosen at random from the 
multiviiiate normal distribution (Naylor et al., 1966) defined by the 
mean vector and the covariance matrix supplied on input for each cover 
type. The result Is a simulated, spatially-referenced population of 
resource parameter vectors. A key assumption (untested) of the SPM is 
that this population representation reflects the actual spatial distri- 
bution of resource parameter values in the inventory area. The validity 
of this assumption depends upon the accuracy of the resource parameter 
estlaiates supplied to the program. 

Program CLUSTER 

Elementary sampling units are g rouped" togethe r "to ~form sa nf~un i ts 
with program CLUSTER. In the current SPM implementation, the elementary 
units are defined as the cells In the spatial class assignment file. 

The dimensions of the clusters and the sample grid origin are specified 
by the user. Cluster sample units created in this manner are assigned 
to strata by CLUSTER in one of two ways. The first method, known as the 
"pHurality" assignment rule, ass'?ns a cluster unit to the stratum (con- 
sisting of some combination of cover types) having the largest number of 
cells within the cluster. A second method uses the simulated resource 
parameter values thems/slves to compute an average parameter vector for 
the cluster. The pregram then assigns the cluster to the stratum having 
a mean vector Riost similar to the average ve(^tor computed for the cluster. 
The population size of first siage sample units classified into each 
stratum by either method is summarized for input to SUMT. 

CLUSTER also produces dispersion statistics and probabilities of 
sample unit selection required by SUMT. Both first sample stage and 
second sample stage ccvariance sv^trices are computed directly from the 
cluster structure and the probability selection method used. First 
sample stage selection probabilities are provided for each first stage 
sample unit In each stratum. Variable selection probsbilities, if used, 
are defined as the area occupied by cover types of interest in a given 
first stage unit relative to the total area in that unit. 
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The SUMT nonlinear programming procedure, developed by Mylander et 
el* (1971 end based upon the theory presented by Ftacco and McCormick 
(1968), is used to solve for the optimum allocation of sampling units 
among sample stages and strata. In addition to the CLUSTER output file, 
ground- to- aerial photograph correlation estimates must be specified by 
the user, as well as cost coefficients appropriate to each sample design 
alternative to be evaluated. These data include the costs of ground 
measurement, photograph measurement, satellite imagery measurement, and 
travel within and between sample units, all exp'^*ssfd on a per sample 
unit basis. Previous surveys or pilot surveys can be used as a guide In 
specifying these correlations and cost coefficients. 

SUHT minimizes the total variable cost (TVC) for a given sample 
design subject to variance constraints on the estimators and to con- 
straints on the sample sizes. The decision variables solved for in this 
process are the sample sizes by stage and stratum required to simul- 
taneous I y meet the precision requirements on each resource parameter 
evaluated by the SPM. 

The solution is optimal in the sense that, for the given design, 
total variable cost is minimized subject to variance constraints. SUMT 
provides a global minimum solution to the sample size allocation problem 
provided that the objective function (total variable cost In this case) 

Is convex, and that the constraining functions (precision levels, maximum 
and minimum sample sizes) are concave. However, Bracken and McCoriiiick 
(1968) have shown that a local minimum may be a global minimum even in 
some cases where the convexity assumptions are not met. 

Ranking Alternative Designs 

To aid In the selection of the sample design to be Implemented on 
the Inventory area, all alternative designs Should be evaluated by SPM. 
These designs should be ranked by total variable cost based on the number 
of sample units at each stage. However, variable cost ranking alone is 
not sufficient to select a final design. The true "optimal" design 
alternative will be selected from the set of designs having a low total 
variable cost, and which also satisfies fixed budget constraints. Is 
easiest to understand and Implement, can be implemented and completed 
with time constraints, minimizes opportunities for measurement error, 
and provides the resource manager with the most flexibility in satisfying 
ongoing information needs. 


The program SUMT was used by Hazard and Promnitz (I 97 M In a similar. 
*” context where they were determining .the optimal replacement fractions 
for continuous forest inventory with partial replacement. 
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APPENDIX 2-C 

WOODLAND/ FOREST TRANSECT LENGTH 
OPTIMIZATION PROCEDURES AND CALCULATIONS 

The objective of this task was to estimate the expected coefficient 
of variation by percent cover of timber volume In the forested area, wood volume 
In the woodland area, and estimate the number of trees per plot expected In each 
type by percent cover class. This Information would then be used In the sample 
size and plot size calculation In the optimization model. 

The simulation was done using the system described In Appendix 3-A 
(Population and Sampling Simulation). 

The data used In the Pinyon-Juniper (P-J) population simulation was 
obtained from Howell.^ A composite stand table was entered. Table 1, Into the 
computer and four populations that covered the expected range of conditions In 
the study were simulated: 

(1) 15X of normal 

(2) 33% of normal 

(3) 66% of normal 

(4) 100% of normal . 

These percentages of a "normal" stand created stands with crown closure percents 
of; 

(1) 10.5% CC 

(2) 23.1% CC 

(3) 46.2% CC 

(4) 70% CC. 

The Ponderosa Pine stands. Table 8, were simulated using volume Information by 

2 

crown diameter from Dllworth, and Information on trees per acre and volume using 
the BLM's preliminary data from Black Rock Inventory. 

It appeared from the results that the P-J simulated population did 
represent the range of conditions that would be encountered In the area. However, 
the Ponderosa Pine simulation produced excessively high volumes per acre for the 
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percent cover and the number of trees per acre used. The coefficient of variation 
for volume and number of trees per transect obtained from the simulation provided 
a reasonable estimate for the purpose of sample design In both the Pine type and 
the P-J type, Tables 2, 3, 4, 9. 

The coefficient of variation (cv) followed the expected shape when it 
was calculated as a function of line length, clearly showing the gain In precision 
obtained by Increasing line length In a given cover percent. The data obtained 
from the simulation runs for P-d was run through a least square regression model 
to obtain the following prediction equation for cv, given percent crown closure 
and line length, 

(2.98 - .391 • L - .00706 • - .419 Vp) 

cv « E 

E = 2.71828 
L * Transect Length 
P - % Crown Closure 

2 

The R for this equation Is 91.8 with 28 degrees of freedom. The standard deviation 
of Y about the regression line Is .15. 

The general form of the model was obtained for O' Regan & Arvanitls' 
and modified to meet the needs of this project. 

The number of trees per transect given transect length and % crown 
closure was obtained from the simulation data. Table 4 and Table 9. This data 
was then used to produce the following least square regression estimate for 
expected number of trees per transect: 

-3 

No. of trees/transect = 1.26x10 (Line Length) * {% Cover) 
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SELECTION OF LINE LENGTH AND ESTIMATION 
OF PARAMETERS FOR THE PLANNING MODEL 


The optimization of line length and number of lines requires that two 

functions be estimated when optimizing simple random sampling and a constraint 

20 

on either budget or precision be established (Arvinetas ). When more complex 
techniques (multi-stage and multi -phase are used, the optimization becomes much 
more complex requiring cost, correlation, auto correlation and variance being 
estimated for many combinations of sampling procedures, plot sizes and clusters of 
plots. The approach taken here was to optimize line length, given the conditions 
expected and the budget constraints. 

A number of summary statistics were calculated from the simulated data to 
aid in the selection of optimum line length (plot size) and to determine the sensi- 
tivity of the model to changes in population parameters or plot size. Table 5 
summarizes three of these parameters: 

1) Coefficient of variation for simple random sampling 

2) Coefficient of variation after correction for double 
sampling (photo ground) 

3) Relative Standard Error of the Estimate given double sampling. 

In addition to the cv and number of trees/ transect prediction equation, 
a number of parameters were set to allow the optimization model to operate. These 
parameters consisted of 1) firm constants and 2) soft constants. The firm 
constants were established for the project and rarely changed. The soft constants 
were set for each run to allow experimentation with various values. 

The constants used for this project included: 

1) Size of Area to be Sampled = 300,000 acres 

2) Allowable Error * 1205$ 

3) Probability Level ® .8 
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4) Range of Line Length ■ 50 to 300 feet 

5) Photo to Ground Correlation at the Plot Level » .85 

6) Intra-Cluster Correlation ■ ,0 

7) Landsat to Photo/Ground Correlation ■ .5 

Cost elements considered and values assumed In the optimization model 

Included: 

Ground Costs 


1) 

Travel between clusters 

2.5 mlles/hr 

2) 

Travel between plots In clusters 

1 .25 mlles/hr 

3) 

Establish a cluster 

15 min/ transect 

4) 

Establish a plot 

3 mln/plot 

5) 

Measure each tree 

8 mln/tree 

6) 

Crew cost 

$346.65/day 

Photo Costs 

A. Acquisition Cost 

1) Telephoto photo cost (Including processing & labeling) $1. 50/plot 

■ 2) 

Wide angle photo cost (Including processing & labeling) .50/each 


3) Aircraft Cost 33,ld/m1n 
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4) Flight crew cost/day $365/day 

5) Miles/minute on line 1.73 miles/minute 

6) Miles/minute between lines 1.83 miles /minute 

7) Minutes to establish each line 3 min/line 


B. Interpretation Cost 


1) Set up flight line 6 min/line 

2) Set up plot 1.5 min/plot 

3) Measure tree .5 min/tree 

4) $/min for interpreter 41.75i/min 

There were 20 optimization model runs made to determine the best combin- 
ation of line length, number of plots per flight line (PSU) to be flown, and the 
number of ground samples to be taken (PSU and SSU) for each allotment, The number 
of PSU to be flown and the number of SSUs per PSU dropped as the length of the transect 
was increased. However, given the constraints placed on photo acquisition by 
the contract, the following combination was considered "best": 


1) 23 

2 ) 2 

3) 150 foot 

4) 5 


flight lines (PSUs) 
plots per PSU (SSUs) 
transect length 
ground PSUs visited 


The expected number of Ponderosa Pine (PP) trees per transect and the 
coefficient of variation (Table 5 and Table 9) for PP board foot volume were within 
the error bounds of the PJ prediction equation. Given the limited amount of data 
for PP it was assumed that the optimization would represent the PP as well. 
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Due to the lack of Information in the range environment, it, was not 
possible to provide the same level of planning for the range inventory. Therefore, 
the recomnendations for sampling in A-2 were made to ensure that adequate data 
would be obtained from this inventory to allow an optimum planning effort to be 
conducted for future inventory work in this environment. This recommendation 
included both ocular and dot count methods for photo interpretation, and sampling 
procedures to ensure sampling over the range of conditions present in an unbiased 
fashion. 
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GROUND DATA COLLECTION INSTRUCTIONS 
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3. Ground Data Collection 


a . Range 

1) Plot location 

Use topo maps and LSP prints to locate plot center. Plot center 
is defined to be the center of the area of stereo coverage of 
the two LSP photos (i.e. a stereo pair) per SSU (ground plot). 

2) Transect and subplot location 

Layout transects and subplots as shown is Figure 1. Proceed 
by temporarily staking plot center point first and then locating 
end points of center transect (transect #3) and temporarily staking. 
Locate end points of transects II and 15 and temporarily stake 
each. Then locate end points of transects 12 and 14 and tempo- 
rarily stake each. Remove stakes after measurements are taken. 

Transect end points and plot center point should be pin pricked 
with annotations on the back of one of the LSP prints for the 
plot. The annotation should include plot number and transect 
numbers, as in Figure 1. Stretch a chain with markings at 1.25 
foot intervals between end points of a given transect before lo- 
cating points. Each marking point on the chain will define a 
sample point. After measuremeiits on a given transect, one of 
the designated SVIH weight measurement subplots (circled plots 
in Figure 1.) should be selected for clipping at random using 
a random number table. Record the number of this subplot on the 
recording form. 

3 ) Transect and subplot measurement 

On all points note the single ground cover component and 
canopy components as in SVIM procedures. 

- On 20 systematically located subplots take SVIM shrub and 
weight characterization plot measurements. These subplots 
(.1 M ) are marked as circled hashmarks in Figure 1. 
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Weight and characteriza- 
tion plots are located at 
every 10th point; the 
starting points for these 
plots are: 

Transect Starting Point 


1 

2 

3 

4 

5 


2 

4 

6 

8 

10 


The starting points are 
determined by counting 
from the first point 
sampled on each line 
(starting points are 
counted from alternating 
ends of the transects 
as shown). 



Figure 1. Range Sample Plot Layout. 
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Clippings should be placed In separate bags by plot and by 
species with complete labels on each bag Including date, 
plot no., transect no., name of clipper, species. 

- Take supplementary plot data (see SVIM forms). 

Two points, at least 20 meters apart, appearing In stereo, 
will be located on the ground and on the photo and pinpricked 
on the large-scale photography. The level distance between 
the two points will be measured to 1 centimeter and recorded 
on the form. Using an Abney or similar device, the slope 
to 1 degree will be measured between the two points with 
the direction and degree of slope recorded on the form. 

4) Data Coding Forms 

SVIM forms with minor changes. See Tables 1 through 4. 

5) Field crew-organization, training, scheduling 

a) 2 crews of 2 people each. In addition, one field crew 
supervi sor/coordi nator . 

b) Initial 3 day training period with follow-up training 
based on supervisor plot rechecks. 

c) Crew scheduling will be arranged so that PSU's selected 

In a given cover type In a given allotment will be visited 
sequentially over the growing season. 

b. Woodland/Forest 

1) Plot location: as In range example. 

2) Plot layout: stake transect ends as annotated on LSP. 

3) Plot measurement 

a) Woodland 

Record on coding forms the following Information for each 
pinyon and juniper tree touching the transect 
- species and tree number; trees are numbered sequentially 
from the ^ieglnnlng of the transect. Tree location 
and number will be pinpricked and annotated on corre- 
sponding LSP print. A clump of trees that appears as 
one crown on the photo should be annotated as such on 
the field forms as In multiple stem trees discussed 
below; 

2-F-4 



diameter ground height; 

diameter stump height (1 foot above ground level) If 

more than one main stem occurs below 1 foot, then record 

diameter for each and label stems a,b,c,etc.; 

height as estimated ocularly after using clinometer 

calibration exercise on 3 or 4 trees; 

number of commercial stems 

crown width, major ani minor axes; 

bore 2 pinyon and 2 juniper trees per plot. Bore on 
side facing plot center when possible, otherwise on 
uphill side only. Count number of rings to center of 
tree and record. Measure length of core from outer 
ring to tree center and record. Measure length of 
core for last 10 years growth and record. Refer to 
Instructions In BLM Extensive Forest Inventory Manual 
for more Information on boring and ring counting tech- 
niques (see Table 6). 
mortality data; 

supplementary plot data (see coding forms) 
b) Forest 

Record the following data for each ponderosa pine tree ^4" 

diameter breast height (dbh) touching the transect, 
species and tree number; tree mimbering should be 
sequential from the beginning of the transect. 

Annotate tree location and number on LSP print. A 
clump that appears as one tree on the photo should be 
Identified as such on the field form, 
stratify trees measured Into 2 inch dbh classes (e.g. 

3.0 - 4.99 = 4", 5.0 - 6.99 = 6", etc.) and into 
open versus closed/partially closed canopy situations. 
Using a random number table select one tree from each 
dbh stratum, canopy class stratum combination and bore 
at breast hight. Perform core measurements (age and 
10 year radial growth) defined In the woodland case. 


2-P-5 



Also measure and record bark thickness. In addition, 
height to commercial top and to tree top. 

Mortality as defined In BLM Extensive Forest Inventory 
Manual : 


GUIDE FOR ESTIMATING TIME SINCE MORTALITY - INTHRMOUNTAIN 


Died within past 5 years 

[ Species ; 

Died more than 5 years ago 

Some foliage left 

Ponderosa 

No foliage 

50+% twigs left 

pine 

Big limbs gone 

Most branches left 

(122) 

Much bark sloughing 

Most bark left 


50/^ or less twigs & 50% or 
less branches left 


4) Data Coding Forms 


5) Field Crew Organization 
Similar to range. 
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Table 5. 


BO RING AND RING COUNTING INSTRUCTIONS 
(From 6LM Ext. For. Inv. Man.) 

Item 20. Ten-Year Radial Growth 20. 10-Yr Radial Growth 

Measure the past 10- year radial growth from an increment core bored 
iinnediately below the point of d.b.h. measurement and at right angles to the 
bole. Count 10 rings in from the outer end of the wood core and mark the 
point. Measure the length of the core from the outer edge of the last complete 
summer wood ring to the outer edge of the summer wood rings 10 years ago, 
using the ruler graduated to 1/20 inches. Record the reading to nearest 
l/20th inch as a 3-digit code. Example: 6/20 code 006, 21/20 code 021. 

Bore on side facing plot center when convenient, otherwise on uphill 

side. 


Measure only for growth sample and site trees. 

When not using the first tree of any 2-inch d.b.h. class as a growth 
sample. tree because of abnormal d.b.h., make a boring for growth on the next 
tree encountered in that size class. Enter the radial growth in item 20 tor 
the tree actually bored. 

If the only tree encountered on a location in a particular 2-1 nch 
class has an abnormal d.b.h., enter a dash in item 20 and the words "abnormal 
d.b.h." in the remarks column. This recording in remarks column also applies 
to all other abnormal d.b.h. trees. Also enter a dash in item 20 for all 
rough and rotten trees and all dead trees. 

Item 37. Tree Age at DBH 37. Tree Age 

Record as a 3-digit code the age taken at breast height from increment 
borings for each of the Growth Sample or site trees only. For all other trees, 
enter a dash (*) in item 37. 
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Table 5. --Continued. 


When determining the age of a tree having Determining age of 
a radius greater than the length of an Increment .Ir*^ 

borer, use the following procedure: 


Bore Into tree as far as possible, extract 
core, and count the rings. Measure the diameter of 

the tree and divide by two, then subtract the bark OHiGiNAu PAGE IS 

thickness. This gives the radius of the wood part POCS QUALITY 

of the tree. Measure the length of the core and 

subtract from the radius of wood to determine how 

much longer the core would have to be to reach the 

pith. Count the number of rings In the last 2 Inches 

and extrapolate to the center. Add this to the 

ring count on the extracted core. See example. 


Example: Determine the age of a Douglas-fir 

58.6 Inches d.b.h. having bark 

thickness of 2.0 Inches when a core 

14.0 Inches long has 110 rings and 

the last 2 Inches have 10 rings. 

58.6/2 • 29.3 (radius of 

wood and bark) 

29.3 - 2.0 ■ 27.3 (radius of wood) 

27.3 - 14.0 ■ 13.3 Inches (short of 

hitting center) 


10/2 


• 5 rings per Inch 


5 X 13.3 « 66.5 rings 66 

110 +66 ■ 176 years old, or 

age group code 19. 


Table 6 
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APPENDIX 2-G 


COMPUTER MENU DESCRIPTIONS 
LEVEL 3 CLASSIFICATION - 14 CLASSES 
BASED ON 10 ACRE AGGREGATION ANOVA RESULTS 


This appendix contains menu descriptions of the '1 level 3 vegetation 
framework categories that were originally identified through ANOVA runs on 
the 117 spectral classes of the Landsat classification. After identifying 
the categories, the classification was aggregated into a 10 acre mapping 
resolution, i.e., all areas less than 10 acres in size were reassigned to 
categories represented by surrounding areas greater than 10 acres. The 
descriptions for the categories represented by the 10 acre mapping were 
created by a new set of ANOVA runs of photo samples against the classification 
and reported here. 

The second ANOVA step was performed because the reassignment of a 
computer class could result in a photo sample being associated with a new 
class. Note that Summary Class 7, Blackbrush, had no photo samples at the 
10 acre aggregation, but 12 photo samples before the aggregation step was 
performed. Consequently, the vegetation descriptions for those categories 
with no photo samples at 10 acre resolution have been provided based on 
their original ANOVA results. 
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NASA/BLM APT PHASE II 

ARIZONA CLASSIFICATION VEGETATION/TERRAIN MENU DESCRIPTION 


Sumary Class; 1 No. of Photo Samples: ^ X of Total: 0 

Nam* Agriculture No. of Acres: 1735 No. of Hectares: 694 

Spectral Classes: 15-24, 114 % of Area: <1 

I. VEGETATION 


(% Cover by 

A. Trees Mean Std.Err. 

Ponderosa Pine 

Pinyon Pine 

Juniper 

Other Tree 


C. Shrubs — Mountain Mean Std.Err. 

Gambel’s Oak 

Turbinella Oak 

Other Shrub 

E. Grasses Mean Std.Err. 

Perennials 

Annuals 

G. Non-Vegetation Mean Std.Err. 

Barren (Rocky) 

Barren (Sandy) 

Water 

Shadow 


Species) 

B. Shrubs — Desert Mean Std.Err. 

Creosote 

Bursage 

Blackbrush 

Big Sagebrush 

Other Shrub 

D. Riparian Woodland Mean S td . Err . 

Cottonwood 

Willow 

Other Shrub 

F. Cactus Mean Std.Err. 

Yucca 

Other Cactus 
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COMPUTER CLASS 1 (Agricu lture) GEOBLOCK Arizona 



Z BAS NO ASPECT (FLAT) 





NASA/BLM APT PHASE It 

ARIZONA CLASSIFICATION VEGETATIOW/TERRAIN MENU DESCRIPTION 


Summarv Class : 2 

_ No. 

of Photo Samples: 183 X of Total: 

9 


Name: Pond erosa Pine 

No. 

of Acres: 

20.364 No. of Hectares; 

8146 


Forest 

Spectral Classes: 

64-65 

X of Area: 1 





I. VEGETATION 





(Z Cover by Species) 



A. Trees 

Mean 

Std.Err. 

B. Shrubs — Desert 

Mean 

Std.Err. 

Ponderosa Pine 

26.1 

1.8 

Creosote 



Piny on Pine 

4.3 

.7 

Bursage 



Juniper 

9.0 

.9 

Blackbrush 



Other Tree 

2.0 

.3 

Big Sagebrush 

3.1 

• 6 




Other Shru5 

1.8 

.3 

C. Shrubs — Mountain 

Mean 

Std.Err. 

D. Riparian Woodland 

Mean 

Std.Err. 

Gambel's Gsk 

12.4 

1.2 

Cottonwood 



Turblnelia Oak 

.8 

.2 

Willow 



Other Shrub 

3.3 

.7 

Other Shrub 

.1 

.1 

E. Grasses 

Mean 

Std.Err. 

F. Cactus 

Mean 

Std.Err. 

Perennials 



Yucca 



Annuals 

2.4 

._8 

Other Cactus 



G. Non-Vegetation 

Mean 

Std.Err. 




Barren (Rocky) 

6.4 

1.1 




Barren (Sandy) 

27.2 

1.3 




Water 

.1 

a 




Shadow 

1.1 

._6 
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NASA/BLM APT PHASE II 

ARIZONA CLASSIFICATION VEGETATION/TERRAIN MENU DESCRIPTION 

Summary Class; 3 No. of Photo Samples: 1.222 !K of Total: 58 

Name: Plnvon-Juniper No. of Acres: 1,056,976 No. of Hectares; 

viooaiahd 

Spectral Classes: 33,35,41-43,46-47 , % of Area: 41 

52-53,55-59,61-63,66,83-85 

I. VEGETATION 
(% Cover by Species) 


A. Trees 

Mean 

Std.Err . 

B. Shrubs — Desert 

Mean 

Std.Err 

Ponderosa Pine 

1.0 

.1 

Creosote 



Piny on Pine 

6.5 

.3 

Bursage 



Juniper 

17.8 

.4 

Blackbrush 

.1 

.0 

Other Tree 

1.5 

.1 

Big Sagebrush 

7.0 

.3 




Other Shrub 

2.1 

.1 

C. Shrubs — Mountain 

Mean 

Std.Err. 

D. Riparian Woodland 

Mean 

Std.Err 

Gambel's Oak 

3.2 

.3 

Cottonwood 



Turbinella Oak 

1.6 

.2 

Willow 



Other Shrub 

7.0 

.4 

Other Shrub 


— 

E. Grasses 

Mean 

Std.Err. 

F. Cactus 

Mean 

Std.Err 

Perennials 

.3 

.1 

Yucca 

.1 

0 

Annuals 

2.7 

.3 

Other Cactus 


— 

G. Non-Vegetation 

Mean 

Std.Err. 




Barren (Rocky) 

19.7 

.7 




Barren (Sandy) 

29.1 

.8 




Water 


- 




Shadow 

.4 

.1 
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NASA/BLM APT PRASE II 

ARIZONA CLASSIFICATION VEGETATION/TERRAIN MENU DESCRIPTION 


Summary Clasa; 4 No, of Photo Samples: 0 % of Total: Q_ 

Name; Riparian Woodland No. of Acres: 368 No. of Hectares: 147 

Spectral Classes: 25-26 % of Area: <1 


I . VEGETATION 


(X Cover by 

A. Trees Mean Std.Err. 

Ponderosa Pine 

Pinyon Pine 

Juniper 

Other Tree 

} 

C. Shrubs — Mountain Mean Std.Err. 

Gambel a Oak 

Turblnella Oak 

Other Shrub 

E. Grasses Mean Std.Err. 

Perennials 

Annuals 

G. Non-Vegetation Mean Std.Err. 

Barren (Rocky) 

Barren (Sandy) 

Water 

Shadow 


Species) 

B. Shrubs — Desert Mean Std.Err. 

Creosote 

Bursage 

Blackbrush 

Big Sagebrush . 

Other Shrub 

D. Riparian Woodland Mean Std.Err. 

Cottonwood 

Willow 

Other Shrub 

F. Cactus Mean Std.Err. 

Yucca 

Other Cactus 
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ELEVATION COMPUTER CLASS 4 (Ripari an Woodland) GEOBLOCE 



,Z HAS HO ASPECT (FLAT) 







NASA/BLM APT PHASE II 

ARIZONA CLASSIFICATION VEGETATION/TERRAIN MENU DESCRIPTION 


Summary Claaa; 5 No. of Photo Samples: 410 X of Total: 19 

Name: Upland Desert Shrub No. of Acres: 565,623 No. of Hectares: 226,249 

Spectral Classes: 1-10.12-14, 74-75, X of Area: 

81,86-87,89-100 

1. VEGETATION 
(X Cover by Spu.'^ies) 


A. Trees 

Mean 

Std.Err. 

B. Shrubs — Desert 

Mean 

Std.Err 

Ponderosa Pine 



Creosote 

5.8 

.3 

Pinyon Pine 



Bursage 

2.6 

.2 

Juniper 

.1 

.1 

Blackbrush 

.5 

.2 

Other Tree 



Big Sagebrush 

.8 

.2 




Other Shrub 

7.2 

.3 

C. Shrubs — Mountain 

Mean 

Std.Err. 

D. Riparian Woodland 

Mean 

Std.Err 

Cambel’s Oak 



Cottonwood 

.1 

.1 

Turbine 11a Oak 



Willow 



Other Shrub 



Other Shrub 



E. Grasses 

Mean 

Std.Err. 

F. Cactus 

Mean 

Std.Err 

Perennials 

.1 

0 

Yucca 

.4 

.1 

Annuals 

12.1 

.8 

Other Cactus 

.9 

.1 

G. Non-VeRetatlon 

Mean 

Std.Err. 




Barren (Rocky) 

iixl 





Barren (Sandy) 


_i^ 




Water 






Shadow 
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ELEVATION COMPUTER CLASS 5 (Do land Desert Shrub) GEOBLOCE Arizona 



Z HAS NO ASPECT (FLAT) 






NASA/BLM APT PHASE II 

ARIZONA CLASSIFICATION VEGETATION/TERRAXN :4ENU DESCRIPTION 

Suamary Claast 6 No. of Photo Samples: 535 X of Total: 25 

Name: Great Basin Sage- No. of Acres: 90A.389 No. of Hectares: 361,756 

brush Spectral classes a. 1. 28-30, 32. 34. 38. X of Area: 

40,44 ,48-50 , 54 , 101-103 , 105-106 , 108 ,110,112, 115-117 

I. VEGETATION 
(X Cover by Species) 


A,. Trees 

Mean 

Std.Err. 

B. Shrubs— Desert 

Mean 

Std.Err 

Ponderosa Pine 



Creosote 

.1 

0 

Plnyon Pine 

.2 

0 

Bursage 

.3 

.1 

Juniper 

1.9 

.2 

Blackbrush 

.4 

.2 

Other Tree 



Big Sagebrush 

12.3 

.7 




Other Shrub 

8.7 

.4 

C. Shrubs — Mountain 

Mean 

Std.Err. 

D. Riparian Woodland 

Mean 

Std.Err 

Gambel's Oak 



Cottonwood 



Turblnella Oak 

.2 

.1 

Willow 



Other Shrub 

.8 

.2 

Other Shrub 

— 


E. Grasses 

Mean 

Std.Err. 

F. Cactus 

Mean 

Std.Err 

Perennials 

3.2 

.5 

Yucca 

1.5 

.1 

Annuals 

14.9 

1.0 

Other Cactus 

.3 

.1 

G. Non-Vegetation 

Mean 

Std .Err. 




Barren (Rocky) 

27.2 

1.4 




Barren (Sandy) 

26.5 

1.2 




Water 






Shadow 

1.4 

.5 
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jl HAS NO ASPECT (FLAT) 







NASA/BLM APT PHASE II 

ARIZONA CLASSIFICATION VEGETATION/TERRAIN MENU DESCRIPTION 


Sunmary Class; 7 No. of Photo Samples: n % of Total: 

Name: Blackbrush No. of Acres: 7468 No. of Hectares: 2987 

Spectral Classes: 88 % of Area: <1 

I. VEGETATION* 

(% Cover by Species) 


A. Trees 

Mean 

Std.Err. 

B. Shrubs— Desert 

Mean 

Std.Err. 

Ponderosa Pine 



Creosote 

2.3 

1.3 

Pinyon Pine 



Bursage 

.8 

.8 

Juniper 

2.4 

2.3 

Blackbrush 

9.2 

8.2 

Other Tree 

.3 

.2 

Big Sagebrush 

1.2 

.8 




Other.' Shrub 

_^6 

_L_2 

C. Shrubs — Mountain 

Mean 

Std.Err. 

D. Riparian Woodland 

Mean 

Std.Err. 

Gambcl's Oak 



Cottonwood 



Turbivtella Oak 



Willow 



Other Shrub 



Other Shrub 



E. Grasses 

Mean 

Std.Err. 

F. Cactus 

Mean 

Std.Err. 

Perennials 



Yucca 

_l8. 


Annuals 


2.0 

Other Cactus 

1,5 

.8 

G. Non-Vegetation 

Mean 

Std.Err. 




Barren (Rocky) 

37 = 9 

11.1 




Barren (Sandy) 

24.0 

9,2 




Vjiater 






Qihadow 

5.0 

5.0 





* % ._o*rsr b«!.sed on full resolution data and 12 photo samples. 
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ORIGINAL PAGE 
OF POOR QUALITY 




,Z HAS NO ASPECT (FLAT) 




NASA/BLM APT PHASE II 

ARIZONA CLASSIFICATION VEGETATION/TERRAIN MENU DESCRIPTION 


I 


Summary Class; 8 No. of Photo Samples: 13 Z of Total: 1 

Name; Other Tall Shrub No. of Acres: 7421 No. of Hectares: 2969 

Spectral Classes: HI % of Area: <1 


I. VEGETATION 
<% Cover by Species) 


A. Trees 

Mean 

Std.Err. 

B. Shrubs — Desert 

Mean 

Std.Err 

Ponderosa Pine 



Creosote 



Ploy on Pine 



Bursage 



Juniper 

.5 

.5 

Blackbrush 



Other Tree 



Big Sagebrush 






Other Shrub 

4.5 

1.3 

C. Shrubs — Mountain 

Mean 

Std.Err. 

D. Riparian Woodland 

Mean 

Std.Err 

Gambel's Oak 



Cottonwood 



Turbinella Oak 



Willow 



Other Shrub 



Other Shrub 



E. Grasses 

Mean 

Std.Err. 

F. Cactus 

Mean 

Std.Err 

Perennials 

13.7 

4.9 

Yucca 

2.8 

2.0 

Annuals 

33.1 

7.2 

Other Cactus 



G. Non-Vegetation 

Mean 

Std.Err. 




Barren (Rocky) 

6.5 

6.5 




Barren (Sandy) 

38.8 

5.5 




Water 






Shadow 
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NASA/BLM APT PHASE II 

ARIZONA CLASSIFICATION VEGETATION/TERRAIN MENU DESCRIPTION 


Summary Glass: 9 No. 

of 

Photo Samples: 

12 X of Total: 

1 

Name: Half Shrub No. 

of 

Acres: 8788 

No. of Hectares: 

3515 

Spectral Classes: 


45 

X of Area: <1 



I . VEGETATION 




(% Cover by Species) 



A. Trees 

Mean 

Std.Err. 

B. Shrubs — Desert 

Mean 

Std.Err 

Ponderosa Pine 



Creosote 

2 2 

1.2 

Pinyon Pine 



Bur sage 

2.1 

2.1 

Juniper 



Blackbrush 



Other Tree 



Big Sagebrush 

3.6 

2.2 




Other Shrub 

10.9 

1.7 

C. Shrubs — Mountain 

Mean 

Std.Err. 

D. Riparian Woodland 

Mean 

Std.Err 

Gambel's Oak 



Cottonwood 



Turblnella Oak 



Willow 



Other Shrub 



Other Shrub 



E. Grasses 

Mean 

Std.Err. 

F. Cactus 

Mean 

Std.Err 

Perennials 

1.3 

.7 

Yucca 

.8 

.4 

Annuals 

12.0 

3.5 

Other Cactus 

.4 

.2 

G. Non-Vesetatlon 

Mean 

Std.Err. 




Barren (Rocky) 







Barren (Sandy) 

11., .5 





Water 






Shadow 
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ELEVATION COMPUTER CLASS 9 (Half Shrub) GEOBLOCK Arizona 



X HAS HO ASPECT (FLAT) 





NASA/BLM APT PHASE II 

ARIZONA CLASSIFICATION VEGETATION/TERRAIN MENU DESCRIPTION 


Suimnary Class: 10 

No. 

of Photo Samples: 24 X of Total: 

1 


Name: Oakbrush 

No. 

of Acres: 

3638 No. of Hectares: 

1455 


Spectral Classes: 

60,113 

X of Area: <1 





I . VEGETATION 
(X Cover by Species) 



A. Trees 

Mean 

Std.Err. 

B. Shrubs — Desert 

Mean 

Std.Err 

Ponderosa Pine 

.7 

1.6 

Creosote 



Pinyon Pine 



Bur sage 



Juniper 

.3 

.5 

Blackbrush 



Other Tree 

.3 

.5 

Big Sagebrush 

.5 

1.2 




Other Shrub 

3.1 

2.4 

C. Shrubs — Mountain 

Mean 

Std.Err. 

D. Riparian Woodland 

Mean 

Std.Err 

Gambel's Oak 

12.8 

9.4 

Cottonwood 



Turbinella Oak 



Willow 



Other Shrub 

50.0 

11.3 

Other Shrub 



E. Grasses 

Mean 

Std.Err. 

F. Cactus 

Mean 

Std.Err 

Perennials 



Yucca 



Annuals 



Other Cactus 



G. Non-Vegetation 

Mean 

Std.Err. 




Barren (Rocky) 
Barren (Sandy) 

11.5 

9.1 




Water 






Shadow 







2-G-20 



COMPUTER CLASS 10 (Oakbrush) GEOBLOOC 



Z HAS NO ASPECT (FLAT) 







NASA/BLM APT PHASE II 

ARIZONA CLASSIFICATION VEGETATION/TERRAIN MENU DESCRIPTION 


Suinnary Class; 11 No. of Photo Samples: 10 X of Total: <1 

Name: Other Mti\. Shrub No. of Acres: 6258 No. of Hectares: 2503 

Spectral Classes: 27,31,39,51 % of Area: <1 


I. VEGETATION 
(% Cover by Species) 


A. Trees 

Mean 

Std.Err . 

B. Shrubs— Desert 

Mean 

Std.Err 

Ponderosa Pine 



Creosote 



Piny on Pine 



Bursage 



Juniper 

4.2 

2.7 

Blackbrush 



Other Tree 

.7 

.5 

Big Sagebrush 

.4 

.4 




Other Shrub 

2.8 

1.3 

C. Shrubs — Mountain 

Mean 

Std.Err. 

D. Riparian Woodland 

Mean 

Std.Err 

Gambel's Oak 

^ 

.2 

Cottonwood 



Turblnella Oak 

3.4 

2.2 

Willow 



Other Shrub 

14.1 

3.7 

Other Shrub 



E. Grasses 

Mean 

Std.Err. 

F. Cactus 

Mean 

Std.Err 

Perennials 



Yucca 

.4 

.4 

Annuals 

5.9 

2.6 

Other Cactus 



G. Non-Vegetation 

Mean 

Std.Err. 




Barren (Rocky) 

54.8 

9.9 




Barren (Sandy) 

13.1 

8.9 




Water 







Shadow 


2-G-22 








NASA/BLM APT PHASE II 

ARIZONA CLASSIFICATION VEGETATION/TERRAIN MENU DESCRIPTION 


Suosaary Class: 12 

No. 

of Photo Samples: 0 Z of Total: 

0 


Name: Perennial Grassland No. 

of Acres: 

16,368 No. of Hectares: 

6547 


Spectral Classes: 

36-37.77.104.107 % of Area: 1 





109 






I. VEGETATION* 
(X Cover by Species) 



A. Trees 

Mean 

Std.Err. 

B. Shrubs — Desert 

Mean 

Std.Err. 

Ponderosa Pine 



Creosote 



Plnyon Pine 



Bursage 



Juniper 



BlackbrUish 



Other Tree 


.2 

Big Sagebrush 






Other Shrub 

_^4 


C. Shrubs — Mountain 

Mean 

Std.Err. 

D. Riparian Woodland 

Mean 

Std.Err. 

Gambel's Oak 



Cottonwood 



Turblnella Oak 



Willow 



Other Shrub 



Other Shrub 



E. Grasses 

Mean 

Std.Err. 

F. Cactus 

Mean 

Std.Err. 

Perennials 

13.0 

5.8 

Yucca 



Annuals 

36.5 

5.7 

Other Cactus 



G. Non-VeRetatlon 

Mean 

Std.Err, 




Barren (Rocky) 






Barren (Sandy) 

42.6 

6.1 





Water 

Shadow 


* % cover based on ANOVA on full resolution data with 10 photo samples. 


2-G-24 



15 I HAS NO ASPECT (FLAT) 





NASA/BLM APT PHASE II 

ARIZONA CLASSIFICATION VEGETATION/TERRAIN MENU DESCRIPTION 


Sumnary Claea; 13 No. of Photo Samplaa: 0 Z of Total: 0 

Name: Barren Land No. of Acres: 0* No. of Hectares: 0 

Spectral Classes: 76,76-80 X of Area: 0 


I. VEGETATION 
(Z Cover by Species) 


A. Trees 

Mean 

Std.Irr. 

B. Shrubs— Desert 

Mean 

Std.Err 

Ponderosa Pine 



Creosote 



Piny on Pine 



Bursage 



Juniper 



Blackbrush 



Other Tree 



Big Sagebrush 






Other Shrub 



C. Shrubs— Mountain 

Mean 

Std.Err. 

D. Riparian Woodland 

Mean 

Std.Err 

Ge^bel's Oak 



Cottonwood 



T'jrblnella Oak 



Willow 



Other Shrub 



Other Shrub 



E. Grasses 

Mean 

Std.Err. 

F. Cactus 

Mean 

Std.Err 

Perennials 



Yucca 



Annuals 



Other Cactus 



G. Non-Vegetation 

Mean 

Std .Err. 




Barren (Rocky) 
Barren (Sandy) 

— 





Water 






Shadow 







* No area greater than 10 acres exists. 


2-G-26 


ELEVATION 



Z HAS NO ASPECT (FLAT) 






( 


NASA/BLM APT FHA&E II 

ARIZONA CLASSIFICATION VEGETATION/TERRAIN MENU DESCRIPTION 


Summary Clasa; 14 No. of Photo Samplaa: 5 X of Total: <1 

Name: Water No. of Acres: 6569 No. of Hectares: 2628 

Spectral Classes: 67-7J,82 % of Area: <1 


1. VEGETATION 
(% Cover by Species) 


i 


A. Trees 

Mean 

Std.Err. 

B. Shrubs — Desert 

Mean 

Std.Err. 

Ponder osa Pine 



Creosote 



Plnyon Pine 



Bursage 



Juniper 



Blackbrush 



Other Tree 



Big Sagebrush 






ether Shrub 



C. Shrubs — Mountain 

Mean 

Std,Err . 

D. Riparian Woodland 

Mean 

Std.Err. 

Gambel ' s Oak 



Cottonwood 



Turbinella Oak 



Willow 



Other Shrub 



Other Shrub 



E. Grasses 

Mean 

Std.Err. 

F. Cactus 

Mean 

Std.Err. 

Perennials 



Yucca 



Annuals 



Other Cactus 



G. Non-Vegetation 

Mean 

Std .Err. 




Barren (Rocky) 






Barren (Sandy) 






Water 

ion. 





Shadow 







: i 

j C 

i. 

I 

i 


2-G-28 



ELEVATION COMPUTER GLASS lAfVater) GEOBIOCK Arizona 



Z HAS NO ASPECT (FLAT) 








APPENDIX 2-H 


( 

AREA ESTIMATES BY PASTURE AND ALLOTMENT 
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APPENDIX 2-1 


COMPUTER CLASS MENU DESCRIPTIONS: VERIFICATION DATA 
LEVEL 3 CLASSIFICATION - 14 CLASSES 
BASED ON 10 ACRE AGGREGATIONS ANOVA 


This appendix contains the menu descriptions developed to verify 
those descriptions presented In Appendix 2-G. The 586 photo samples used 
In the ANOVAs for these descriptions were a selected subset from the total 
2985 photo plots for the area and were not used In the original descriptions. 
In some cases, there were Insufficient photo plots available to perform both 
the vegetation framework category description and the separate verification 
step as well. Ai a result, only seven of the 14 summary classes could be 
verified. 
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NASA/BLM APT PHASE II 

ARIZONA CLASSIFICATION VEGETATION/TERRAIN MENU DESCRIPTION 


Sumary Class 


No. 

No. 


Name : P onderosa Pine 
Forest 

Spectral Classes: 


of <?hoto Samples: 4S X of Total: 

of Acres: , No. of Hectares: 

X of Area: 


1. VEGETATION 
(X Cover by Species) 


A. Trees 

Mean 

Std.Err . 

B. Shrubs — Desert 

Mean 

Std.Err. 

Ponderosa Pine 

20.6 

2.7 

Creosote 



Piny on Pine 

5.3 

1.6 

Bur sage 



Juniper 

10.3 

1.8 

Blaekbrush 



Other Tree 

1.8 

.3 

Big Sagebrush 
Other Shrub 

2.9 

2.3 

J.:l 

.6 

C. Shrubs — Mountain 

Mean 

Std.Err. 

D. Riparian Woodland 

Mean 

Std.Err. 

Gambel's Oak 

13.1 

2.2 

Cottonwood 



Turblnella Oak 

1.2 

.7 

Willow 



Other Shrub 

5.2 

1.9 

Other Shrub 



E. Grasses 

Mean 

Std.Err. 

F. Cactus 

Mean 

Std.Err. 

Perennials 



Yucca 



Annuals 

1.6 

.7 

Other Cactus 



G. Non-Vegetation 

Barren (Rocky) 
Barren (Sandy) 
Water 

Mean 

3.6 

31.4 

Std.Err. 

1.4 

2.4 




Shadow 
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NASA/BLM APT PHASE II 

ARIZONA CLASSIFICATION VEGETATION/TERRAIN MENU DESCRIPTION 

Sunmary Class! 3 No. of Photo Samplaa: 3T0 X of Total: _______ 

Name: Pinyon-Juniper No. of Acres: No. of Hectares: 

Spectral Classes: X of Area: 


I. VEGETATION 
(X Cover by Species) 


A. Trees 

Mean 

Std.Err. 

B. Shrubs— Desert 

Mean 

Std.Err 

Ponderosa Pine 



Creosote 



Pinyon Pine 

5.5 

.5 

Bur sage 

.1 

.1 

Juniper 

17.3 

.7 

Blackbrush 

.1 

,1 

Other Tree 

1.5 

.1 

Big Sagebrush 

9.3 

.8 




Other Shrub 

2.5 

.3 

C. Shrubs— Mountain 

Mean 

Std.Err. 

D. Riparian Woodland 

Mean 

Std.Err 

Gambel's Oak 

2.8 

.6 

Cottonwood 



Turbinella Oak 

2.2 

.3 

Willow 



Other Shrub 

6.2 

.7 

Other Shrub 



E. Grasses 

Mean 

Std.Err. 

F. Cactus 

Mean 

Std.Err 

Perennials 

.3 

.1 

Yucca 

.2 

.0 

Annuals 

2.9 

.5 

Other Cactus 

.0 

.0 

G. Non-Vegetation 

Mean 

Std.Err. 




Barren (Rocky) 

16.4 

1.4 




Barren (Sandy) 

32.0 

1.5 




Water 






Shadow 

.3 

.3 
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NASA/BLM APT PHASE II 

ARIZONA CLASSIFICATION VEGETATION/TERRAIN MENU DESCRIPTION 


Sunnary Glaaa! 5 No. of Photo Samplaa: 101 * of Total: 

Nama: Upland Desert Shru b No. of Acres: No. of Hectares: 

Spectral Classes: X of Area: 


I. VEGETATION 
(% Cover by Species) 


A. Trees 

Mean 

Std.Err. 

B. Shrubs — Desert 

Mean 

Std.Err 

Ponderosa Pine 



Creosote 

5.6 

.6 

Piny on Pine 



Bur sage 

■ 

,4 

Juniper 



Blackbrush 

__ii 


Other Tree 



Big Sagebrush 

.7 





Other Shrub 

-Aa 


C. Shrubs— 'Mountain 

Mean 

Std.Err. 

D. Riparian Woodland 

Mean 

Std.Err 

Gambel's Oak 



Cottopvood 



Turblnella Oak 



Willow 

.6 

.6 

Other Shrub 



Other Shrub 



E. Grasses 

Mean 

Std.Err. 

F. Cactus 

Mean 

Std.Err 

Perennials 



Yucca 

.4 

.1 

Annuals 

13.8 

2.0 

Other Cactus 

.9 

.2 

G. Non-Vegetation 

Mean 

Std.Err. 




Barren (Rocky) 

33.4 

3.6 




Barren (Sandy) 

29.3 

3.5 




Water 

1.8 

1.3 




Shadow 

3.0 

1.7 
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NASVBLM apt phase II 

ARIZONA CLASSIFICATION VEGETATION/TERRAIN MENU DESCRIPTION 


Sunnary Claaa: 6 

No. of Photo Samplea: 

118 X of Total: 

Name: Great Basin Sage- 

No. of Acres: 

No. of Hectares: 


brush Clagses: X of Area: 

1. VEGETATION 




(X Cover by Species) 



A. Trees 

Mean 

Std.Err. 

B. Shrubs — Desert 

Mean 

Std.Err 

Ponder osa Pine 



Creosote 

.1 

.0 

Piny on Pine 



Bur sage 

.1 

.0 

Juniper 

1.2 

.4 

Blackbrush 

.6 

.4 

Other Tree 

.1 

.1 

Big Sagebrush 

12.6 

1.4 




Other Shrub 

8.2 

.8 

C. Shrubs — Mountain 

Mean 

Std.Err. 

D. Riparian Woodland 

Mean 

Std.Err 

Gambel's Oak 



Cottonwood 



Turbinella Oak 

.1 

.1 

Willow 



Other Shrub 


.J 

Other Shrub 



E. Grasses 



Mean 

Std.Err. 

F. Cactus 

Mean 

Std.Err 

Perennials 

2.9 

1.0 

Yucca 

1.6 

.3 

Annuals 

13.0 

2.0 

Other Cactus 

.3 

.1 

G. Non-VeRetation 

Mean 

Std.Err. 




Barren (Rocky) 

29.2 

3.2 




Barren (Sandy) 

25.2 

2.5 




Water 






Shadow 

4.0 

Ilz 
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NASA/BLM APT PHASE II 

ARIZONA CLASSIFICATION VEGETATION/TERRAIN MENU DESCRIPTION 


Suonary data; 7 No. of Photo Sanplaa: 1 Z of Total: 

Nana: Blackbrush No. of Acrea: No. of Hectares: 

Spectral Classes: X of Area: 


1. VEGETATION 
(X Cover by Species) 


A. Trees 

Mean 

Std.Err. 

Ponderosa Pine 
Piny on Pine 
Juniper 
Other Tree 

— 

— 

C. Shrubs— -Mountain 

Mean 

Std.Err. 

Gambel's Oak 
Turbine 11a Oak 



Ocher Shrub 



E. Grasses 

Mean 

Std.Err. 

Perennials 

Annuals 

12.0 

- 

G. Non-VeRetatlon 

Kean 

Std .Err. 

Barren (Rocky) 
Barren (Sandy) 
Water 
Shadow 

ZM 

1 .11 ■■ 


B. Shrubs — Desert 

Mean 

Std.Err. 

Creosote 

6.0 


Bur sage 

4.0 

t* 

Blackbrush 



Big Sagebrush 



Other Shrub 

3.0 

- 

D. Riparian Woodland 

Mean 

Std.Err. 

Cottonwood 



Willow 



Other Shrub 



F. Cactus 

Mean 

Std.Err. 

Yucca 



Other Cactus 


- 


2 - 1-6 
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I 

1 NASA/BLN APT PHASE! II 

^ ARIZONA CLASSIFICATION VEGETATION/TERRAIN MENU DESCRIPTION 

Sunnary Claaa; 8 No. of Photo Samplan: 8 Z of Total: ____ 

Name: Other Tall Shrub No. of Acres: _________ No. of Hectarea: , 

Spectral Claaaes: X of Area: _________ 


1. VEGETATION 
(X Cover by Species) 


A. Trees 

Mean 

Std.Err. 

B. Shrubs — Desert 

Mean 

Std.Err 

Ponder osa Pine 



Creosote 



Piny on Pine 



Bursage 



Juniper 



Blackbrush 


1 

Othei Tree 



Big Sagebrush 






Other Shrub 

7.1 

1.7 

C. Shrubs — Mountain 

Mean 

Std.Err. 

D. Riparian Woodland 

Mean 

Std.Err 

Gambel’s Oak 



Cottonwood 



Turbinella Oak 



Willow 



Other Shrub 

— 

— 

Other Shrub 

— 


E. Grasses 

Mean 

Std.Err. 

F. Cactus 

Mean 

Std.Err 

Perennials 

ZiLL 


Yucca 

1.0 

1.0 

Annuals 

19.4 

.111 

Other Cactus 



G. Non-Veaetatlon 

Mean 

Su' Jrc. 




Barren (Rocky) 

11.4 

11.4 




Barren (Sandy) 

40.6 

8.1 




Water 






Shadow 

— 

— 





( 
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NASA/BLM APT PHASE II 

ARIZONA CLASSIFICATION VEGETATION/TERRAIN MENU DESCRIPTION 


Sunnary Claaat 10 Wo. of Photo Sat:>p.les! 3 Z of Total: 

Name: OakbrUSh No. of Acres: No. of Hectares: 

Spectral Classes: X of Area: 


I. VEGETATION 
(X Cover by Species) 


A. Trees 

Mean 

Std.Err. 

B. Shrubs—Desert 

Mean 

Ponderosa Pine 



Creosote 


Plnyon Pine 



Bursage 


Juniper 

.3 

.3 

Blackbrush 


Other Tree 

1.7 

1.7 

Big Sagebrush 

- - 




Other Shrub 

1.7 

C. Shrubs — Mountain 

Mean 

Std.Err. 

D. Riparian Woodland 

Kean 

Gambel's Oak 

13.3 

13.3 

Cottonwood 


Turbinella Oak 



Willow 


Other Shrub 

41.7 

10.9. 

Other Shrub 


E. Grasses 

Mean 

Std.Err. 

F. Cactus 

Mean 

Perennials 



Yucca 


Annuals 



Other Cactus 


G. Non-VeRetatlon 

Mean 

Std.Err. 



Barren (Rocky) 

13.3 

13.3 



Barren (Sandy) 

28.0 

14.5 




Water 

Shadow 


< 

2 - 1-9 


Std.Err. 


1.7 


Std.Err. 


Std.Err. 
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RANGE PRODUCTIVITY ESTIMATES 
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RANGE FORAGE ESTIMATES BY PASTURE 
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RANGE FORAGE ESTIMATES BY PASTURE 
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RANGE FORAGE ESTIMATES BY PASTURE 



9.435 3.818 1.8 18.0 20.1 









































































RANGE FORAGE ESTIMATES BY PASTURE 


























RANGE FORAGE ESTIMATES BY PASTURE 



TOTAL 45,876 18,566 8.6 15.3 17.2 








































































RANGE FORAGE ESTIMATES BY PASTURE 
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RANGE FORAGE ESTIMATES BY RASTURE 
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MUSE FORAGE ESTIWiTES BY PASTURE 
PAUTABLE SPECIES ONLY - CURRENT AMOUNTS AVAILABLE OR UNAVAILABLE 
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Appendix 3-A 


Wildland Resources Inventory Planning Model 

(Resource Inventory Services) 

Design of an optiirum inventory is an extremely complex 
process which requires trade offs or corn. '»mises among 
numerous alternatives. 

This involves the distribution of the inventory effort 
in Loth a horizontal and vertical direction. 

Decisions in the horizontal direction include the deter- 
mination of the size of the sample plot, the shape of a 
sample plot, the clustering size and arrangement of sample 
plots and the number of clusters required at each level of 
data (ground, large scale photo and Landsat). 

Decisions in the vertical direction include the selec- 
tion among alternative sources of data (Landsat, small 
scale photo, large scale photo and ground) and the portion- 
ing of effort between the selected lines of data to achieve 
an optimum mix of data from the selected sources. 

The model utilizes cost, time and population factors to 
determine the optimum allocation of effort in both of these 
directions. It will handle inventories based on ground 
only, photo stratification followed by ground sampling, 
photo stratification followed by large scale photo sampling 
with ground subsampling, and stratification followed by 
complete enumeration with photo sampling and ground 
subsampling. It will optimize plot size, number of plots 
per cluster, number of clusters, and the ratio of large 
scale photo effort to ground effort given the required inputs. 

The optimization can be done by specifying either 
desired precision and probability level which will estimate 
the effort required or by specifying budget limitations 
which will then estimate the expected optimum effort and 
predict the expected precision and probability level given 
the fixed budget. 

Given a number of inputs associated with the cost of 
collecting data and the characteristics of the resource to 
be inventoried, the model outputs estimates of optimum: 

li plot size, 

2. cluster size (number of plots/cluster), 

3. sample size (number of clusters) , 

4. distance between plots in cluster, 

5. ratio of photo to ground in double sampling, 

I 6. expected budget for each element, 

^ 7. hours required for each element. 
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Injiut rccjuiroincnts; 

Jost elements 
. for aircraft 

a. speed between cransects, 

b. speed on transact, 

c. average time to sot up transect, 

d. average ferry tii'.e fro;., operating bane to study 
site, 

e. maximum safe f iying tirAO, 

f. cost o£ photo product, 

g. format of image.ry, 

h. desired stereo coverage, 

i. ground crew ccsi, 

j. flight crow cost, 

k. aircraft hourly rate, 

l. expected overhesd rate, 

cxj)octed error rate in acquisition (fi'Mire into) 
. for ground data collection 

a. road nctv;ork classification, 

b, expected ground speed 

i. walking to flight lino, 

ii. walking belv.-een flig)it lines, 

iii. driving be iv.’een flight lii.n.s, 

c. time from camp to general area, 

d, expected failure rate, 

c. time to measure a sample unit or itcmi within the 
unit in forests, 

f. time to set up for 

i. a transect, 

ii. a plot, 

g. data entry ti:,.c; per plot, 
i. per sample unit, 

h. hourly rate in dollars for crev;, 

i. per diem, 

ii. per expected day length, 

iii. per expected overhead rate, 

. for photo interpretation costs 

a, time to set up for flight line transe;;t.y (I\ ) 

b, time to set up for plot within transect (GSU/ 

c, time to measure attributes on plot 

d, expected useful hours/day for interpreter. 


Population characteristics 

1. expected stratified variance as a fimction of p'^ot 
size, 

. the physical size of the unit to be inver Aoric'’ , 

. the expected auto correlation as a function of ; .ot 
size and spacing, 

4. the expected correlation or coefficient of determina- 

tion between the enumeration estimate derived fx'orA 
Landsat or other imagery and the photo/greund esti- 
mates and the first phase estimate from photo/ground 

5. the correlation between the estimates made ^rom the 

first phase sample photography and the second pbvise 
estimates from tlio ground. 


SUMMARY OF FUNCTIONAL CAPABIJjlTIES 


In addition to providing a scientifically and statis- 
tically sound framework for selecting the best sampling and 
measurement scheme ^ the system interactively guides the designer 
through the inventory planning process. The system asks for 
pertinent information needed to optimize the design, therefore 
minimizing the chance of overlooking important constraints or 
assumptions that affect the best design. 

The system consists of six major elements: 

1. population generator 

2. population sampler 

3. sample summarizer 

4. prediction estimator 

5. sample optimizer 

6. sample selector. 

Population Generator 

The population generator takes information on the 
expected conditions of the areas to be inventoried and generates 
a forest stand in the computer to the specifications of the 
designer. 

The population can be generated in two basic ways: 

1. the entry of x,y coordinates and description of 

each tree 

2. the simulation of the population from a stand table. 

In some experimental forest areas, detailed information is 
known for each tree in the forest, including location and 
description. In this case, the detailed information can be 
entered in the proper format and used as the population in 
the sampling and summary routines to experiment with sampling 
procedure optimization. 




More commonly, however, there is very limited infor- 
mation available on an area and the objective is to select 
the best sampling procedures to be used in an inventory. In 
this case, the population can be simulated by obtaining a 
candidate stand table for the area and estimate the range of 
conditions that may be encountered in the area. First, the 
candidate (normal) stand table is entered in the system. The 
system then asks for information to allow the designer to 
create a range of conditions that may be encountered in the 
study. 

Population Sampler 

The purpose of this set of functions is to simulate 
population sampling procedures commonly used in forest 
inventory to aid in the selection of optimum sampling procedures 
and parameters. The types of plots currently available for 
investigation are: 

1. fixed radius plots with selectable radius 

2. horizontal point sampling, with selectable basal 
area factors, 

3. horizontal line sampling with selectable basal area 
factors and line length 

4. line intersect sampling with selectable transect 
lengths . 

Sample Summarizer 

This function provides two basic capabilities: 

1. the summary of the data obtained from the 
sampling runs by plot 

2. the summary of the plot data into sample estimates 
for the population. 
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The plot sununaries are generated by taking the 
individual cree measurement and running them through the data 
reduction algorithm for th»» sampling method used. The algorithm 
uses either table lookup pi.ocedures or prediction equ^\tions 
to generate the tree estimates and the appropriate mathematical 
model to generate the plot summaries. The population sununaries 
are generated by taking the results from the plot summaries and 
calculating the appropriate population estimates based on the 
sample data. 

At each stage in the summary process, the investigator 
has the option of creating permanent disk files for future use, 
temporary files for use in the next step, or the creation of 
computer listing of summary data. 

Prediction Estimator 

This functional portion of the system takes informa- 
tion from the Sample Summarizer and generates estimators for 
parameters needed by the Sample Optimizer. These estimators 
include : 

1. expected value of inventory items of interest, 

2. expected value of characteristics of the 

population that influence cost, and 

3. variance of the expected values. 

These estimators are derived from characteristics 
easily obtained from small pilot studies or preliminary strati- 
fication of the area to be inventoried. They include; 

1. percent crown closure 

2. estimated tree height 

3. estimated crown diameter 

4. spatial distribution of the trees. 
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The output of the Prediction Estimator can either be 
a functional or tabular representation of the relationship 
between the observed characteristics and the parameters needed 
for sample optimization. 

Sample Optimizer 

This is the portion of the system that allows the 
inventory planner to select from among the alternatives the 
"best" data collection methods and sampling scheme. There 
are two basic optimization alternatives: one is to allocate 
the resources between alternative sources of data to give the 
most precise answer for a fixed budget, and two, to minimize 
the budget to achieve a desired level of precision by 
determining the optimum level of effort to be allocated to 
each of the sources of data. 

In this program the designer can specify: 

o varying levels of budget 

o varying levels of precision 

o varying levels of probability. 

Alternatives that can be evaluated include simple random 

sampling, stratified random sampling, systematic sampling, 
and complete enumeration followed by multistage sampling, 
multiphase sampling, and combined multistage multiphase sampling. 
The evaluation can include sampling with equal and varying 
probability at each phase or stage. 

The optimizer requires the output from the prediction 
estimator and expected cost information for the area to be 
inventoried. The cost data includes: 

1. cost to travel between clusters 

2. cost of data gathering per plot 

3. cost of paperwork per plot 

4. cost of data analysis per cluster of plots 

for ground travel, photo acquisition and photo interpretation. 



f 


Additionally, there are correlation characteristics 
that must be estimated and provided to the optimizer. The 
correlations are: 

1. the correlation between estimates (photo and 
ground, photo and photo, ground and ground) 

2. the autocorrelation between sample units by 
unit size and distance between units. 

The input requirements for optimization are dependent of the 
level of design sophistication required. 

When only simple random sampling is being considered, 
plot size and number of plots are optimized given the precision 
and probability statement or maximum budget constraint. The 
inputs required include: 


1. 

travel cost 


2. 

setup cost 


3. 

measurement cost 


4. 

data analysis cost 


5. 

size of area 


6. 

desired precision 


V. 

probability level 


8. 

anticipated average 


9. 

anticipated variance versus 

plot size function. 

When 

simple random sampling with 

clusters of plots 


(multistage sampling) is being evaluated the following 
additional information is required: 

1. autocorrelation between plots (SSU) in the cluster 
given plot size and distance between plots in the 
cluster 

2. cost of establishing cluster (PSU) 

3. cost of travel between plots in the cluster. 

f 
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I When multiphase sampling (double sampling) is being 

evaluated the following additional information is required: 

1. correlation between estimates (measurements) 
made on attributes of interest from each type 
of data considered 

2. cost of making the estimates from each of the 
data types. 

These costs must be broken down by acquisition cost, fixed cost 
and variable cost by plot size or by cluster and plot size if 
multistage sampling is being evaluated. 

When multiphase sampling with complete enumeration 
at the first phase is being considered, the relationship between 
the cost of the enumeration and the expected correlaUion between 
the first and :irecond phase estimates is required. This 
information allows the selection of the appropriate level of 
effort at the first phase. 

When a multiple resource inventory is being designed 
the correlation between the various resource characteristics 
of the population that are to be inventoried and the cost of 
data collection allows the optimization of the inventory by 
taking into account the natural relationships that exist 
between these resource characteristics, This optimization 
reduces the number of measurements to be taken on highly 
correlated characteristics and allocates the effort to the 
measurement of the most cost effective set of characteristics 
in the seunple. 
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Sample Selector 

This portion of the aystem aids in the selection of 
samples given the sampling scheme desired. It allows the 
implementation of simple random sampling, systematic sampling, 
stratified random sampling and probability proportional to 
size sampling. The function samples fzom lists generated 
by the system based on descriptions of the inventory areas 
to be sampled. The output from the function includes the x,y 
coordinates of the sample unit, the strata identifier and the 
selection probability. 



MAJOR STATISTICAL FORMULAS USED IN THE OPTIMIZATION 


number of plots per cluster (SSU's) 

total number of flight lines 

number of flight lines selected 

index for flight lines 

index for plots within a flight line 

volume per acre for plot j of flight line i 


(volume for flight line i) = S 

^ n 

(average volume per flight line) =* Z Y./n 

i ^ 

(average per plot) «= Ey^/nm = Y/M 
Z Z (Yi.-Y) 

i i i 


n'M-; 


= variance among the elements 


Z (y.-Y) 
i ^ 


/ “ variance among the cluster totals 

\n“x)M 




(M-l)S' 


= intraclusv^er correlation 


SJ = 


1 2 

S tl + (M-l)Ppl » variance of the sample mean 


p *= (simple correlation coefficient) = 
xy 


E(Yj^-y) (X^-X) 
(Sx) (Sy) 
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variance arrtong the element 


4 - - 

given a complete enumeration followed by a random 
sample 

2 2 

S (1 - (1“^)) ■ double sample variance 

among the elements of y given 

1 A 2 ^ 

“ — i_ _ B ratio of ground to photo effort 




t 


cv 


cost of obtaining photo data (per plot) 
cost of obtaining ground data (per plot) 
student’s "t” value 

2 _ 

(coefficient of variation) “ S / y 


AR 


allowable error 


n 


^2 2 
t cv 

^ 

AE 


Reference for formulation: 

Cochran; Sampling Techniques 2nd Edition; 
John Wiley & Sons, Inc; 1963. 



